(12) INTERNATIONAL APPLICATION PUBLISHED UNDER THE PATENT COOPERATION TREATY (PCT) 



(19) World Intdlectuai Property Organizatioa 
Internationa] Bureau 




liimmiifiiiiiiiiiiiiin 



(43) InterDational Publication Date (10) International Publication Number 

28 February 2002 (28.02.2002) pCT WO 02/16956 Al 



<SI) loteroatioaai Pateot Classificatioo'': 
33/56 



(}01R 33/48, 



(21) loteraatiooal Application Nnmben PCT/USO 1/25954 

(22) Interaational FiUng Date: 21 August 2001 (21.08.2001) 



CZ, DE. DK. DM, DZ, EQ EE, ES, R. GB, GD, GE, GH, 
GM, HR, HU. ID. IL. IN, IS, JP, KE. KG, KP, KR, KZ, LC, 
LK, LR, LS, LT, LU. LV, MA, MD, MG, MK, MN, MW, 
MX. MZ, NO, NZ. PL, FT, RO. RU. SD. SE, SO, SI, SK, 
SU TJ, TM, TR. TT, TZ, UA. UG. US. UZ, VN. YU. ZA, 
ZW. 



(25) Filing Language: 

(26) Publication Language: 



English 
English 



(30) Priority Data: 
60/226,667 



22 August 2000 (22.08.2000) US 



(71) Applicant and 

(72) Inventor: AflLLS, Randett, L. [USAJS]; 7 Weatherfield 
Drive, Newtown, PA 18940 (US). 

(74) Agent: MELCHER, Jeffrey, S.;Manelli. Denison & Sel- 
ler. PLLC. 7th Floor, 2000 M Street N.W.. Washington, DC 
20036-3307 (US). 

(81) Designated States (national): AE, AG, AL. AM. AT, AU, 
AZ. BA, BB. BG, BR. BY, BZ, CA, CH, CN. CO. CR, CU, 



(84) Designated States (regional): ABIPO patent (GH, OM, 
KE, LS, MW, MZ, SD, SL, SZ. TZ, UG. ZW), Eurasian 
patent (AM, AZ, BY, KG, KZ, MD. RU, TJ, TM), European 
patent (AT. BE, CH, CY, DE, DK. ES. H. FR, GB, GR, IE, 
IT, LU, MC, NL, PT, SE. TR). OAPI patent (BP, BJ, CF, 
CG. CI, CM, GA. GN, GQ, GW. ML, MR, NE, SN, TD, 
TG). 

Published: 

— with iniernational search report 

— before the expiration of the time limit for amending the 
claims and to be republished in the event cf receipt of 
amendments 

For two-letter codes and other abbreviations, refer to the "Guid- 
ance Notes on Codes and Abbreviations" appearing at the begin- 
ning of each regular issue of the PCT Gazette 



(54) Title: 4 DIMENSIONAL MAGNETIC RES0NANC:E IMAGING 



mmimrm 

us a -J m CfSSffJt 



COPY 



n 



snmt 

IQCAm 



DfU¥E 



mcessett 



r" 



' o aaaaooo 




MAmro 

MfffM 



tr 



IT) 

vo 

^ (57) Abstract: Provided are apparatus and methods for forming a multidimensional image of inanimate or animate objects which 
utilize a magnetization source (1 12) to magnetize a volume of an object to be imaged, a radiation source (118) for applying a radiation 
field to the object to be imaged, an output signal detector (120) for producing output signals in response to the secondary radiation 

Q al a plurality of spatial locations outside of the object as a function of time, processor (126, 126a and 126b) for determining a 
plurality of Fourier components, for associating the Fourier components due to each voxel (14) by phase, and for converting each 

^ set of components into a voxel location, and an image processor (128) for producing an image. 
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4 DIMENSIONAL MAGNETIC RESONANCE IMAGING 

This application claims priority to U.S. provisional application serial No. 
60/226,667, filed on August 22, 2000, the complete disclosure of which is 
5 incorporated h^ein by reference. 

1. Field of the Invention : 

10 This invention relates to method and apparatus for imaging an object, such as a 

body. More particularly, the invention provides high resolution, real-time images in 4 
dimensions with Uttle or no deterioration fix>m motion artifact. 

15 2. Background of the Invention : 

Nuclear Magnetic Imaging (NMR) which is conunonly called magnetic 
resonance intiaging (MRI) entails 1 .) magnetizing a volume widi a constant primary 
magnetic field in a z-direction, 2.) providing a gradient along the axis of the z-directed 

20 field to select a slice in the xy-plane, the plane perpendicular to the direction of the 
primary field, 3.) providing electromagnetic radiation resonant with the Larmor 
firequency of protons in the slice, 4.) providing a pulse of resonant electromagnetic 
radiation to flip the magnetization vector into the transverse plane or plane of the slice, 
and S.) applying a magnetic field gradirat along an axis in the xy-plane of the z-directed 

25 field with excitation at the Larmor firequency to provide phase dispersion of the NMR 
signal along the axis to encode spatial information, and 6.) recording the firee induction 
decay (FID) radio emission signals following excitation, 7.) recording a plurality of such 
FIDs, each recorded following an excitation with a rotated direction of the gradient in the 
xy-plane, and 8.) reconstructing the image bom the plurality of the FIDs. An integer n 

30 of FIDs each having a phase gradient that corresponds to the magnetic field gradient that 
was rotated to n unique directions in the xy-plane comprise a set along two orthogonal 
axes in phase or k-space. A two dimensional Fourier transform of the data set is used to 
reconstruct an n by n pixel image. 

MRI is of primary utility in assessing brain anatomy and pathology. But long 

35 NMR relaxation times, a parameter based on how rapidly excited nuclei relax, have 

prevented NMR from being of utility as a high resolution body imager. The most severe 
limitation of NMR technology is that for spin echo imaging n, the number of free 
induction decays ("FIDs"), a nuclear radio frequency energy emitting process, must 
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equal the number of lines in the image. A single FID occurs over q>pn)ximately 0. 1 
seconds. Not considering the q)in/Iattice relaxation time, the time for the nuclei to 
reestablish equilibrium following an RF pulse, which may be seconds, requires an 
irreducible imaging time of n times 0.1 seconds, which for 512 x 512 resolution requires 

5 £q)proximately one minute per each two dimensional slice. This represents a multiple of 
1 500 tunes longer that the time that would freeze organ movements and avoid image 
deterioration by motion artifact. For example, to avoid deterioration of cardiac images, 
the imaging time must not exceed 30 msec. A method for speeding NMR imaging flips 
the magnetization vector of the nuclei by less than 90 degrees onto the xy-plane, and 

0 records less FIDs, Such a method, known as the flash method, can obtain a 128 x 128 
resolution in ^proximately 40 seconds. Another technique used to decrease imaging 
time is to use a field gradient and dynamic phase dispersion, corresponding to rotation of 
the field gradient, during a single FID to produce imaging times typically of 50 msec. 
Both methods produce a decreased signal-to-noise ratio ("SNR") relative to spin echo 

5 methods. The magnitude of the magnetization vector which links the coil is less for the 
flash case because the vector is flipped only a few degrees into the xy-plane. The echo- 
planar technique requires shorter recording times with a concomitant increase in 
bandwidth and noise. Both methods compensate for decreased SNR by increasing the 
voxel size with a concomitant decrease in image quality. Physical limitations of these 

0 techniques render obtaining high resolution, high contrast vascular images impractical. 

SUMMARY OF THE INVENTION 



It is an object of the invention to provide high resolution multi-dimensional 
images of an object, such as a body, tissue, or working cardiopulmonary system. 

It is a furth^ object of the invention to rapidly acquire the data to provide 
magnetic resonance images of a body with reduced motion artifacts. 

These and other objects of the invention are attained by providing an ^paratus 
for obtaining a magnetic resonance image of a body using data acquired over the three 
spatial dimensions plus time, rather than acquiring data only in time at one receiving 
anterma. A preferred embodiment of the apparatus of the invention includes a radiation 
source for applying a first radiation field having a magnetic component to the body, to 
magnetize the body. The apparatus further includes a source for applying a second 
radiation field to the body, to elicit a radiation field from the body. A detector senses 
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fliis radiation field, and produces a signal that a reconstruction processor employs to 
create ttie magnetic resonance image of the body. 

A NMR image is obtained of a magnetized body fit>m a three-dimensional map 
of the intensity variation of flie NMR signal produced by each voxel of flie magnetized 

5 body, and detected over a three-dimensional volume of space external to the body, herein 
refOTcd to as the "sample ispace-" The data is acquired over three spatial directions plus 
time. Iq an embodim^t, the NMR signals are detected over a three dimensional detector 
array as a function of time. The NMR signals may be sampled at least at the Nyquist 
rate, i.e., at a rate that is twice the highest temporally frequency of the NMR signal and 

10 twice the bluest spatial frequency the Fourier transform of flie NMR image of the 
phantom. Sampling at the Nyquist rate or higher allows the spatial variations of the 
external NMR signal to be acquired. 

In an embodiment, the NMR signal at each detector as a function of time is 
processed by a method such as a Fourier transform operation to give a plurality of 

15 Fourier components each having the same frequency, an mtensity and a phase angle. 
The NMR signal of each voxel at any given detector gives rise to a Fourier component 
with a unique phase angle relative to the Fourier component of any other voxel of the 
phantom at that detector. The set of Fourier components that correspond to the NMR 
signal of a given voxel over the detectors is determined. This may be achieved by using 

20 a first component having a phase angle and calculating the phase angle as a function of 
• spatial position of tiie first detector relative to any other detector and identifying the 
component at each detector having the calculated phase angle. The sets are determined 
for all of the voxels. A NMR image is obtained from the sets wherein each set of Fourier 
components comprises a three-dimensional m^ of the intensity variation of the NMR 

25 signal produced by each voxel of the magnetized body. This practice of flie invention 
preferably employs a Fourier transform algorithm, described in Fourier Transform 
Reconstmction Algoriflun Section, to determine tiie spatial location of each voxel from 
the corresponding set of Fourier components comprises a three-dimensional map of the 
intensity variation of flie NMR signal produced by each voxel of the magnetized body 

30 over tiie sample space. This is repeated for each set to form the NMR image of the 
object. 

NMR images produced according to prior art methods and systems rely on 
applying an additional magnetic field in the direction of the primary field having a 
gradient along an axis in the transverse plane to cause a phase variation of the NMR 
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signal along the axis in the transverse plane. The direction axis of the gradient is varied 
a plurality of times to gives rise to an equivalent number of lines in the reconstructed 
image. In the present invention, the unique phase variation of the NMR signal is 
provided by the combination of I .) the angle 6 suspended between the direction of the 

5 detector and the radial vector, the vector from the dipole to the detector, and 2.) the angle 
0 due to a separation distance r between a voxel and a detector given by the 
wavenumber of the RF field k times r . 

My prior inventions disclosed in U.S. Patent No. 5,073,858 and U.S. Patent 
AppUcation No. 09/191,454, the complete disclosures of both which are incorporated 

1 0 herein by reference, are m part based on the realization that matter having a permeability 
different from that of free space distorts a magnetic flux applied thereto. This property is 
called magnetic susceptibility. An object, herein called a phantom, can be considered as 
a collection of small volume elements, herein referred to as voxels. When a magnetic 
field is applied to the phantom, each voxel generates a secondary magnetic field at the 

15 position of the voxel as well as external to the phantom. The strength of the secondary 
magnetic field varies according to the straigth of the applied field, the magnetic 
susceptibility of the material within the voxel, and the distance of the extanal location 
relative to the voxel. For example, my U.S. Patent No. 5,073,858 and U.S. application 
No. 09/191,454 teache that the net magnetic flux at a point extrinsic to a phantom to 

20 which a magnetic field is applied, is a sum of the applied field and the external 

contributions from each of the voxels. The '858 patent fiirther teaches sampling the 
external flux point by point and employing a reconstruction algorithm, to obtain the 
magnetic susceptibiUty of each voxel from the sampled external flux. 

Unlike the '858 patent that relies on a static response &om a magnetized body to 

25 determme the magnetic susceptibility of the body, my invention disclosed in U.S. Patent 
Application No. 09/191,454 elicits a radiative response fit>m a magnetized body by 
subjecting the body to a resonant radiation field. One embodiment of the '454 invention 
generates a three-dimensional magnetic susceptibiUty image of an object including a 
patient placed in a magnetic field from a three-dimensional map of a radio frequency 

30 (RF) magnetic field extemal to the patient, induced by subjecting selected nuclei of the 
body to a resonant RF field. Application of an RF pulse to the body causes the body to 
emit the RF magnetic flux extemal to the body. A Fourier transform of this extemal flux 
produces its firequency components ("Lannor firequencies"). Each Larmor frequency is 
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used to determine the magnetic susceptibility of the voxels of the body producing that 
Lannor jfrequency. Further, the intensity variation of the Vernal RF field over a three- 
dimensional volume of space is used to determine the coordinate location of each voxel. 
One practice of the inventions disclosed in my U.S. Patent No. 5,073,858 and 

5 U.S. Patent Application No. 09/ 1 9 1 ,454 obtains a three-dimensional magnetic 

susceptibility map of a magnetized body ftom a three-dimensional map of a secondary 
magnetic flux produced by the magnetized body, and detected over a three-dimensional 
volimie of space external to the body, herein referred to as the "sample space." The 
extrinsic magnetic flux is sampled at least at the Nyquist rate, i.e., at twice the spatial 

0 frequency of the highest frequency of the Fourier transform of the magnetic 

susceptibility m^ of the phantom, to allow adequate sampling of spatial variations of the 
external magnetic flux. This practice of the inventions preferably employs a Fourier 
transform algorithm, described in Fouria: Transform Reconstruction Algorithm Section, 
to form the magnetic susceptibility map of the object. 

5 The present invention relates to systems for providing images of distributions of a 

quantity, in a chosen region of the body, by gyromagnetic resonance, particularly nuclear 
magnetic resonance (NMR) techniques. Such techniques may be used for examining 
bodies of different kinds. A particularly beneficial ^plication is the examination of 
patients for medical purposes. Unlike my U.S. Patent No. 5,073,858 and U.S. Patrat 

0 Application No. 09/191,454, the present invention employs nuclear magnetic resonance 
(NMR) to induce a magnetized phantom of essentially constant magnetic susceptibility 
to emit an external radiation having a magnetic field component In particular, 
application of an RF pulse, resonant with selected nuclei of a magnetized body, can 
polarize the nuclei throu^ rotation of their magnetic moments. The polarized nuclei 

5 withm a voxel process about the local magnetic field in the voxel at a Larmor firequency 
determined by the applied magnetic field at position of the voxel. The superposition of 
external RF fields produced by all the voxels of the body creates the total external RF 
field at each detector that is time dependent. The external RF field recorded at the 
detectors as a function of time contains components each having a unique phase angle 

0 relative to other components. Each component corresponds to an emittmg voxel of the 
phantom. The time dependent signal at each detector may be transformed into a series of 
components havmg intensity and phase data. Each set of components of the NMR signal 
over the sample space due to a given voxel is determined finom the phase data and the 
detector positions. The spatial variation of the NMR signal over the sample space is 
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used to determine the location of the voxel in the phantom. This is repeated for all sets 
of components, each corresponding to a' voxel to reconstract the NMR image. 

The radiation source for magnetizing a body to be imaged can be a direct current 
("DC") magnet, including a superconducting magnet. The radiation sources and 
5 amplifiers for applying an RF pulse to the magnetized body are well known in the art, 
and include, but are not limited to, klystrons, backward wave oscillators, Gunn diodes. 
Traveling Wave Tube amplifiers. A preferred embodiment of the invention employs a 
three dimensional array of antennas as detectors for sensing the external RF field. 

One practice of the invention detects the external RF field in the near field region 

1 0 where the distance of a detector sensing radiation from a voxel at a distance r from the 
detector is much smaller than the wavelength A of the radiation emitted by the voxel, 
i.e., r«X (or kr The near fields are quasi-stationary, that is they oscillate 
harmonically as e'^'^ , but are otherwise static in character. Thus, the transverse RF 
magnetic field of each voxel is tiiat of a dipole. In one embodiment, an array of 

1 5 miniature RF antennas sample the external RF field over a three-dimensional voliune of 
space that can be either above or below the object to be imaged. The distance r of the 
detector bom the voxel gives rise to the phase term e'^ of the component of the 
detected RF signal where k is the wavenumber of the NMR signal. The harmonic 
oscillation of each RF dipole is equivalent to the dipole rotating in the transverse plane. 

20 The detector is responsive to a component in this plane. At a point in time, each RF 
dipole is directed at an angle 6 relative to the direction of detection of the detector. The 
phase angle 6 of the RF dipole relative to the direction of detection axis of the detector 
gives rise to a phase angle term e"^^ . In a preferred embodiment, flie simi of the phase 
angles, kr and 9 , are unique for each voxel at each detector. The position of each 

25 detector relative to a different detector may be used to calculated the phase angle of the 
second relative to the first. TTiis may be repeated over all of the detectors to give the set 
of intensities of ttie NMR signal over the sample space due to a voxel. The location of 
each voxel is determined through the spatial variations of the intofisity of the NMR field 
of the set of componrats associated by phase. Thus, the phase of the conq)onents of the 

30 external RF radiation, and the intensity variations of the external RF radiation provide 
the necessary information for providing a NMR image of the magnetized phantom, such 
as a human body. Such a NMR unage can be employed to obtain anatomical images of a 
human body based on selected physiological parameters. 
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In an emboduneat of the present invention, the NMR image of an object 
including a patient placed in a magnetic field is generated fiom a three-dimensional map 
of the transverse resonant radio fiequency (RF) magnetic flux external to the pa&exxt 
The external RF field recorded at the detectors as a fimction of time contains components 

S each having a unique phase angle relative to other components. Each component 
corresponds to an emitting voxel of the phantom. The time dependent signal at each 
detector may be transformed into a series of components having intensity and phase data. 
Each set of components of the NMR signal over the sample space due to a given voxel is 
determined bom the phase data and the detectof positions. The intensity variation of the 

0 transverse RF field over the sample space is used to determine the coordinate location of 
each voxel. The RF field is the near field which is a dipole that serves as a basis element 
to fomi a unique reconstructioa The geometric systena function corresponding to a 
dipole which determines the spatial intensity variations of the RF field is a band-pass for 
k^-k^* Preferably, each volume element is reconstructed independently in parallel with 

5 all other volimie elements such tiiat the scan time is no greater than the nuclear free 
induction decay (^ID) time. 

Secondary Magnetic Field 

} The magnetic moment of each voxel is a magnetic dipole. And the phantom 

can be considered to be a three-dimensional array of magnetic dipoles. At any point 
extrinsic to the phantom, the z-component of the secondary flux, B* , from any single 
voxel is 

) where x, y, and z are flie distances fiom the center of the voxel to the sampling point It 
is shown in APPENDICES I-IV that no geometric distribution of magnetic dipoles can 
give rise to Eq. (1). Therefore, the flux of each magnetic dipole (voxel contribution) 
forms a basis set for the flux of the array of dipoles which comprise the NMR image of 
the phantom. 

) Eq. (1) is a system fimction which gives the magnetic flux output in response to a 

magnetic dipole input at &e origin. The phantom is an array of spatially advanced and 
delayed dipoles weighted according to the magnetic moment of each voxel; this is the 
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bsput function. The secondary flux is the siqpaposition of spatially advanced and 
delayed flux, according to Eq. (1); this is the output fimction. Thus, ttie response of 
space to a magnetized phantom is given by the convolution of Eq. (1) with tiie series of 
wei^ted, spatially advanced and delayed dipoles representing the NMR image of the 
5 phantom. 

In Fourier space, the output fimction is the product of the Fourier transform (FT) 
of the systmi function and the FT of the input fimction. Thus, the system function filters 
die input fimction. The output fimction is the flux over all space. However, virtually all 
of the spectrum (information needed to reconstruct the NMR image) of the phantom 
10 exists in the space outside ofUxG phantom because the system fimction is essentially a 
band-pass filter. This can be appreciated by considering the FT, i/[*p,^J , of Eq. (1): 

V 

where is the spatial firequency in the xy-plane or -plane and A, is the spatial 
frequency along the z-axis. H[k^,k^] is a constant for k^ and k^ essentially equal as 
1 5 demonstrated graphically in FIGURE Ic. 

Band-Pass Filter 

When a static magnetic field Ho with lines in the direction of the z-axis is applied 
20 to an object comprising a matoial containing nuclei such as protons that possess 
magnetic moments, the field magn^zes the material. As a result a secondary field 
superposes the applied field as shown in FIGURE 9. In the applied magnetic field, tiie 
magnetic moments of each nuclei processes about tiie ^lied magnetic field. However, 
the magnetization of any one nucleus is not observed from the macroscopic sample. 
25 Rather the vector sum of the dipole moments torn all magnetic nuclei in the sample is 
observed. Thisbulkmagnetizationisdenotedby the vector M. &i thermal equilibrium 
with the primary field , the bulk magnetization M is parallel to . The 
magnetization vector then comprises magnetic dipole m . The secondary field outside of 
the object (phantom) and detected at a detector 301 is that of a series of magnetic dipoles 
30 centered on volmne elements 302 of the magnetized material. In Cartesian coordinates. 



wo 02/1(956 PCT/US01/2S954 

9 

Ihe secondaty magnetic flux, B*. at die point {x,y,z) due to a magnetic dipole having a 
magnetic dipole moment at the position (x,,yo,Zo) is 



where 1, is the unit vector along the z-axis. The field is the convolution of the system 
function, h{x,y,z) orh(p,^,z) (the left-handed part of Eq. (4)), with the delta function 
(the right-hand part of Eq. (4)), at the position (Ai„>'o,Zo)- A very important theorem of 
Fourier analysis states that the Fourier transform of a convolution is the product of the 
1 0 individual Fourier transforms [2]. The Fourier transform of the system function, 
h(x,y,z) or h(p, ^ ,z), is given in APPENDIX V. 

The z-component of a magnetic dipole oriented in the z-direction has the system 
function, h(x,y,z)^ which has the Fourier transform, //[^j^Ar^,*,], which is shown in 
FIGURE Ic. 

47c[Jk/ + AL^] 

15 ^W' n^'.V^X'l " 



1 



(6) 



The output function, the secondary magnetic field, is the convolution of the 
system function, A(x,y,z)— the geometric transfer function for the z-component of a z- 
oriented magnetic dipole with the input fimction-a periodic array of delta functions each 
20 at the position of a magnetic dipole corresponding to a magnetized volume element 

(2z^-x^-/) 
[x -^y +z J 

The Fourier transform of a periodic array of delta functions (the right-hand side of Eq. 
(7)) is also a periodic array of delta functions in k-space: 



(8) 
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By tiie Fourier Theorem, the Fourier transform of the spatial output function, Eq. (7), is 
the product of the Fourier transform of the system fimction given by Eq. (6), and the 
Fourier transform of the input function given by Eq. (8). 




S In flie special case that 

(10) 

the Fourier transform of the system function (the left-hand side of Eq. (9)) is given by 

H^4n (11) 
Thus, the Fourier transform of the system function band-passes the Fourier transform of 

10 the input function. Both the input function (the right-hand part of Eq. (7)) and its Fourier 
transform (the right-hand part of Eq. (9)) are a periodic array of delta functions. No 
frequencies of the Fourier transform of the input function are attenuated; thus, no 
information is lost in the case where Eq. (10) holds. 

In an embodiment of the present invention, the magnetization vector is rotated 

15 into the transverse plane by an additional RF field Hj . The magnetization vector then 
comprises a rotating magnetic dipole m in the transverse plane. The NMR image may 
be reconstructed by sampling the external field from a series of RF dipoles rather than 
that from a series of static dipoles. hi this case, the Fourier transform of the system 
function also band-passes the Fourier transform of the input function. Thus, the 

20 resolution of the reconstructed NMR image is limited by the spatial sampling rate of the 
secondary RF magnetic field according to the Nyquist Sampling Theorem. 

Reconstruction 

25 The NMR image may be reconstituted using a Fourier transform algorithm. The 

algorithm is based on a closed-form solution of the mverse problem — solving the spatial 
distribution of an array of magnetic dipoles from the measured extrinsic secondary (RF) 
field that is transverse to the magnetic flux that magnetizes flie voxels. The transverse 
RF magnetic field of each voxel is that of a dipole, the maxunum ampUtude is given by 

30 Eq. (1) wherem the Larmor frequency of each voxel is essentially the same, and , the 
magnetic moment along the z-axis, of Eq. (1) corresponds to the bulk magnetization M 
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of each voxel. In terms of the coonlinates of Eq. (l), an array of miniature RF antemias 
point samples the maximum dipole component of the RF signal over the sample space 
such as the half ^ace above (below) the object to be imaged wherein each RF signal is 
has a imique phase shift due to the relative spatial relationship of each voxel and each 

5 detector. The phase shift of each conq)onent and the relative spatial relationship of the 
detectors is used to assign a component from each detector to a set Each set of 
conqionents associated by phase comprises the spatial variations of the intensity of the 
transverse RF field of any given voxel. The intensity variation over the sample space is 
used to determine the coordinate location of each voxeL hi the limit, each volume 

0 element is reconstructed independently in parallel with all odier volume elements such 
that the scan time is no greater than the nuclear free induction decay (FID) time. 

The NMR scan performed on the object to be imaged including a human 
comprises the following steps: 

5 • The magnetic moments of nuclei including protons of the object to be imaged that are 
aligned by the primary field are fiuther aligned by a radio frequency (RF) pulse or 
series of pulses. 

• The strength and duration of the rotating Hj (RF) field that is resonant with the protons 
0 of the magnetized volume and is oriented perpendicularly to the direction of the 

magnetizing field is applied such that the final precession angle of the magnetization 
is 90 ^ (0H, = 90"* ) such that the RF dipole is transverse to the primary magnetizing 
field and perpendicular to the RF magnetic field detector. 

5 • NMR pulse sequences which provide the signals for a or image may be £q)plied. 
For sample a 90 pulse maybe followed by a series of 180 pulses. One sequence 
is the Carr-Purcell-Meiboom-Gill (CPMG) sequence [3], 

• The firee induction decay-signals are recorded. 



• The time dependent signals are Fourier transformed to give the intensity and phase of 
each component. The NMR signal of each voxel at any given detector gives rise to a 
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Fourier component with a imique phase angle relative to the Fonrier component of 
any other voxel of the phantom at that detector. 

• The matrix of Fourier components that correspond to the NMR signal of a given voxel 
5 over the detectors is determined. This may be achieved by using a first component 

having a phase angle and calculating the phase angle as a function of spatial position 
of the first detector relative to any other detector and identifying the component at 
each detector having the calculated phase angle. The matrices are determined for all 
of the voxels. The measurements of the spatial variations of the transverse RF field 
10 of a given matrix is used to determine the coordinate location of each voxel. Thus, 

each matrix of components associated by phase comprises the intensity variation over 
the sample space of the RF field of the bulk magnetization M of each voxel. 

• The Foxuier transfonn algorithm is performed on each set of components over the 
15 detector array to map each bulk magnetization M corresponding to a voxel to a 

spatial location over the image space. The bulk magnetiization map (NMR image; 
.also the input fimction) is given as follows. With respect to the coordinate system of 
Eq. (1), (x, y, and z are the Cartesian coordinates, /n, , the magnetic moment along 
the z-axis, of Eq. (1) corresponds to the bulk magnetization M of each voxel, and B 

20 is the magnetic flux due to the magnetic moment shown in FIGURE 9; the 

relationship to the NMR coordinate system is given in the Reconstruction Algorithm 
Section) tiie origin of the coordinate system, (0,0,0), is the center of the upper edge of 
the phantom. The phantom occupies the space below the plane x, y, z=0(z^0in 
the phantom space), and the sampling points lie above the plane ( z > 0 in the 

25 sampling space). The magnetic flux in the sampling space is given by multiplying 

the convolution of the input function with the system function by the unitary function 
(one for z ^ 0 and zero for z < 0). The input function can be solved in closed-form 
via the following operations: 

30 1. Record the RF NMR signal at discrete points in the samplmg space. Each point is 
designated (x, y, z, RF) and each RF value is an element in matrix A. The time 
dependent signals are Fourier transfonned to give the intensity and phase of each 
component. The NMR signal of each voxel at any ^ven detector gives rise to a Fourier 



wo 02/16956 PCT/USOl/25954 

13 

component with a unique phase angle relative to the Fourier component of any other 
voxel of flie phantom at that detector. The matrix of Fourier components that correspond 
to the NMR signal of a given voxel over the detectors is deteraained. This may be 
achieved by using a first component having a phase angle and calculating the phase angle 
5 as a function of spatial position of the first detector relative to any other detector and 
identifying the component at each detector having the calculated phase angle. The 
matrices are determined for all of the voxels. The measurements of the spatial 
variations of the transverse RF field of a given matrix is used to determine the coordinate 
location of each voxel. Thus, each matrix of components associated by phase comprises 
10 the intensity variation over the sample space of the RF field of the bulk magnetization M 
of each voxel. 

2. Discrete Fourier transform each matrix to obtain each matrix . 

15 3. Multiply each element of each matrix 5, by the corresponding inverse (reciprocal) 
value of the Fourier transform of the system Amotion, Eq. (1), evaluated at the same 
frequency as die element of the matrix A„, This is matrix C.. 

4. Generate matrix by taking the discrete inverse Fourier transform of matrix C„ . 

20 

5. Multiply «ach element of each matrix D„ by the distance squared along the z-axis to 
which the element corresponds to generate the position of the bulk magnetization M of 
voxel n. (This corrects the limitation of the sample space to z ^ 0). 

25 6, In one embodiment, the point spread of the reconstructed voxel is corrected by 

assigning one voxel above a certain threshold with the bulk magnetization M . The other 
voxels are assigned a zero value. 

7. This procedure is repeated for all matrices A^. In the limit with sufficient phase 
30 resolution, each volume element is reconstructed independently in parallel wititi all other 
volume elements such that the scan time is no greater than the nuclear firee induction 
decay (FID) time. 
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- BRIEF DESCRIPTION OF THE DRAWINGS 



FIGURE la is a schematic diagram of an q)paratiis according to the teachings of 
the invOTtion in v/bich a magnet magnetizes a body to be imaged, an RF generator 
S excites the body, and an array of detectors detect radiation emitted by the body in 
response to the RF excitation radiation; 

FIGURE lb is a schematic diagram of an apparBtas according to the teachings of 
the invention that enq>Ioys a three-dimensional array of detectors to detect radiation 
emitted by nuclei of a body to be imaged in response to excitation of the nuclei by a 
1 0 radiation source; 

FIGURE Ic is a plot of the Fourier transform H[k^,ky,k^] of the system function 

K^*yf2) (Eq, (1)) corresponding to the z-component of a magnetic dipole oriented in the 
z-direction in accordance with the invention; 

FIGURE 2 is the plot of the field of a ring of dipoles of radius R and magnetic 
1 5 moment m = 10^ Gcm^ given by Eq. (L 14) as a function of radius R where the position 
of the center of the ring relative to flie detector is the point (0,0,10) in accordance with 
the invention; 

FIGURE 3 is the plot of the field of a ring of dipoles of radius R = 0.2 cm and 
magnetic moment /w = 10* Gcrn^ given by Eq. (1.14) as a function of the distance 
20 between the detector at the origin and the center of the ring at the points (0,0, 
cmtoz = l5cm)m accordance with the invention; 
FIGURE 4 is the plot of the field of a shell of dipoles of radius R and magnetic 
moment m = 10* Gcm^ given by Eq, (11.17) as a function of radius R where the position 
of the center of the shell relative to the detector is the point (0,0,10) in accordance with 
25 the invention; 

FIGURE 5 is the plot of the field of a shell of dipoles of radius R = 0.2 cm and 
magnetic moment m = 10* Gcm^ given by Eq, (11.17) as a function of the distance 
between the detector at the origin and the center of tiie shell at the points (0,0, 
z=4 cm to z = 15 cm) in accordance with the invention; 
30 FIGURE 6 is the plot of tiie field of a sphere of dipoles of radius R and magn^c 

moment m = 10* Gcm^ given by Eq. (IV.16) as a function of radius R where the 
position of the center of the sphere relative to the detector is the point (0,0, 10) in 
accordance with the invention; 
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FIGURE 7 is the plot of the field of a sphere of dipoles of radius R = 0.2 cm and 
magnetic moment m = 10^ Gcm^ given by Eq. (WA6) as a function of the distance 
between die detector at the origin and the center of the sphere at the points (0,0, 
z = 4 €iR z = 15 cm ) in accordance with the invention; 
5 FIGURE 8 shows a typical the nuclear magnetic resonance (4D-MRI) apparatus 

in accordance with the invention; 

FIGURE 9 shows the coordinate system (x, y, z) of Eq. (1) with a primary field 
Hq and the corresponding magnetic dipole both oriented paraUel to the z-axis wherein 
the z-component of the flux due to the z-oriented dipole is measured at a detector 
10 according to Eq. (1) and shows the distances bom the voxel to the detector in accordance 
with the invention; 

FIGURE 10 shows the stationary coordinate system of the nuclear magnetic 
resonance (4D-MRI) ^paratus of FIGURE 8 corresponding to the coordinate system of 
FIGURE 9 in accordance with the invention; 

15 FIGURE 1 1 shows the rotating NMR coordinate system (x^ , yR , Zr ) and the 

stationary coordinate system (x, y, z) of the NMR detector corresponding to the 
coordinate system of FIGURE 9 and FIGURE 10 of a primary field Ho oriented parallel 
to the -axis and the z-axis and the corresponding transverse RF magnetic dipole 
oriented in the Xi^y,^ -plane and periodically parallel to the y-axis wherein the spatial 

20 variation of the RF y-component of the flux due to the RF dipole is measured at a 
detector according to Eq. (1) and shows the distances fijom the voxel to the detector in 
accordance with the invention; 

FIGURE 12 is a schematic of the three dim^ional detector (antemiae) array 
with respect to the stationary NMR coordinate system of FIGURE 10 which corresponds 

25 to tiie coordinate systems shown in FIGURES 9 and 1 1 in accordance with the invention; 
FIGURE 13 shows the general process of reconstruction by reiteration, and 
FIGURE 14 shows the stationary coordinate system (x, y, z) of the NMR detector 
corresponding to the coordinate system of FIGURE 9 and FIGURE 10 of a primary field 
Hq oriented parallel to the -axis and the z-axis and the corresponding transverse RF 

30 magnetic dipole oriented in the x^y^ -plane and periodically parallel to the y-axis 

wherein the spatial variation of the RF y-component of the flux due to the RF dipole is 
measured at a detector according to Eq. (1) and shows the distances and angles between a 
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voxel linear to a first detector, a second nonaligaed voxel, and a second detector in 
accordance with die invention; 

Further details regarding specific derivations, and calculations are provided in the 
attached appendices, wherein: 
5 APPENDIX I is the field produced by a ring of dipoles according to the present 

invention; 

APPENDIX n is the derivation of the field produced by a shell of dipoles 
according to the present invention; 

APPENDIX in is the mathematical proof that the field produced by a shell of 
10 magnetic dipoles is different fiom that of a single dipole according to the present 
invention; 

APPENDIX IV is the derivation of the field produced by a sphere of dipoles 
according to the present invention; 

APPENDIX V is the derivation of the Fourier transform of the system function of 
15 the z-component of the magnetic field fi"om a dipole oriented in the direction of the z- 
axis used in the reconstruction process according to the present invention; 

APPENDIX VI is the derivation ofS = HF® U(kJ convolution torn Eq. (55) 
used in a reconstruction process according to the present invention, and 

APPENDIX Vn is the derivation of the Inverse Transform of Eq. (69) to Give 
20 Inverse Transform 1 , Eq. (69), used in a reconstruction process according to the present 
invention. 



DETAILED DESCRIPTION OF THE INVENTION 



25 An exemplary embodunent of a nuclear magnetic resonance (4D-MRI) apparatus 

1 10 accordmg to the teachings of the present invention is shown in FIGURE 1 a The 
apparatus 1 10 includes a magnet 112, such as a superconducting magnet, that provides a 
primary or magnetizing field, to magnetize a body 1 14 to be imaged. A magnetometer 
1 16 determines the primary or magnetizing field in the volume to be occupied by the 

30 body, i.e., in the image space, in flie absence of tibe body. One practice of the invention 
utilizes a magnetometer that employs NMR of protons in water for determining the 
primary radiation field at multiple points in the image space. In an embodiment, the 
primary field is uniform as recorded by magnetometer 1 16. 
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A radiation source 118, such as a radio fiequency generator, applies an RF pulse 
to the body in combination with RF antennae 1 1 8a when the body is placed in die 
magnetic field, to excite and thereby polarize selected nuclei of the body. The excited 
nuclei emit an RF radiation that a plurality of detectors 120, disposed in a plane above or 
S below the object, can detect The excitation pulse can be selected to rotate the 

magnetization of tiie nuclei, preferably by 90 degrees, with flie respect to the primary 
field, hi such a case, the RF radiation that the excited nuclei emit is primarily along a 
direction perpendicular to the plane of the detectors, hi an embodiment using techniques 
known by those skilled in the art, RF pulse sequraces are qiplied to genmte the data for 

10 r, and/or NMR images. The detectors can be selected to respond only to components 
of a magQetic field perpendicular to the plane in which they reside. Thus, such detectors 
can detect the emitted RF field without interference fi-om other components of the 
magnetic field to permit a unique reconstruction of the NMR image. In an embodiment, 
a drive mechanism 122 moves the detectors 120 in a direction perpendicular to the plane 

15 of the detectors, to sample the external RF field over a three-dimensional volume. The 
separation of the detectors and the step size of the movement of the detectors along a 
direction perpendicular to the plane are selected such that the detectors sample the 
external RF field over a three-dimensional volume, i.e., the sample volume, at least at the 
Nyquist rate. Preferably, a three dimension array of detectors is used to sample the RF 

20 field at the Nyquist rate. Such a sampling advantageously allows obtaining the NMR 
image of the body. An embodiment, employs an impedance-matched array of RF 
antennas that are time-multiplexed to reduce cross talk among them. The RF field may 
be sampled synchronously or the field may be sampled at known times so that the phase 
at any given detector may be related to that at any other detector. 

25 The process of magnetization of selected nuclei of the body can be better 

understood by referring to FIGURE 1 1 that shows a magnetizing field H^, applied to a 
voxel 14a of tiie body 1 14. The RF excitation field is selected to be in a direction 
perpendicular to the field Hg , and its magnitude is designated as . NMR active nuclei 
of the voxel, such as protons, possess both angular momentum and a magnetic moment. 

30 Thus, the vector sum of the magnetic moments of all such NMR active nuclei present in 
the voxel 14a give rise to the bulk magnetization of the voxel 14a. The bulk 
magnetization vector M (not shown) designates the collective contribution of a selected 
type of NMR active nuclei to the magnetization of the voxel 14a which corresponds to a 
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RF parallel magnetic moment m . The RF ^citation field is selected to be in resonance 
with the selected type of nuclei, e.g., protons, of the voxel 14a, to rotate the bulk 
magnetization vector M . In a rotating fitune representation, designated in FIGURE 1 1 
as(X|i, y^y ) and weU known to those sldUed in NMR, (he magnetization vector M 
5 rotates about H, so long as the RF radiation is present The rate of rotation of the 
magnetization vector about the ^plied field depends on the gyromagnetic ratio of 
the a£fected nuclei and the magnitude of the excitation field, Hj . The duration of the RF 
radiation can thus be selected to cause a rotation of the magnetization vector, initially 
aligned along H^, onto the Xi^yi^^lane, for example by a 90** rotation. After the RF 

10 excitation field is turned off, the rotated magnetizafion M precesses about at the 
position of the voxel 14a. 

The precession of the magnetization M of each voxel about produces a 
radiating dipolar field corresponding to a magnetic moment m that has its maximum 
intensity along the y^ direction. The precession fi'equency of the magnetic moment 

1 3 called the Larmor fi'equency is dependent on the magnetic flux of the primary field 

which is uniform in a preferred embodiment A detector 20a can be, for example, an RF 
antenna pointing along the y direction to respond selectively to the radiating RF field 
from all of the voxels. This external RF field is recorded as a function of time at each 
over the entire sample space. The NMR signal at each detector as a fimction of time is 

20 processed by a method such as a Fourier transform operation to give a plurality of 
Fourier components each having an intensity and a phase angle. The NMR signal of 
each voxel at any given detector gives rise to a Foxuier component with a unique phase 
angle relative to the Fourier component of any other voxel of the phantom at that 
detector. Each set of components of the NMR signal over the sample space due to a 

25 given voxel is determined 6om the phase data and the detector positions. The spatial 
variation of the NMR signal over the sample space is used to determine the location of 
the voxel in the phantom. This is repeated for all sets of componoits, each 
corresponding to a voxel to reconstruct the NMR image. Preferably, the position of each 
voxel is reconstructed independentiy in parallel with all other voxels such that the scan 

30 time is no greater the time for the excited nuclei to return to their unexcited state called 
the nuclear fi-ee induction decay (FID) time. This may be achieved by using sufBciently 
high spatial and temporal sampling rate such that each set of components associated by 
phase corresponds to a single voxel. 
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With further reference to FIGURE la, an analog to digital converter ("A/D 
converter") 124 converts the analog signal outputs of flie detectors 120 into digital 
signals. A preferred embodiment of the invention utilizes at least a 12 bit A/D converter 
to digitize the output signals. A reconstruction processor 126 receives the digital signals, 
5 and determines the NMR image of the body. The reconstruction processor 126 includes 
a Fourier transform processor 126a that obtains the phase components of the external RF 
field. In addition, a spatial locator 126b, which is a part of the reconstruction processor 
126, employs the variation of the maximum intensity of the external RF field over the 
sample space (the three dimensional space sampled by the detectors) to locate the 

10 positions of the voxels in the image space producing a particular set of components of 
the external RF field associated by phase, in a manner described in detail below (See thb 
Reconstruction Algorithm Section). An display processor 128 displays a two- 
dimensional or a three-dimensional image corresponding to the NMR image of the body. 
The spatial locator I26b employs an algorithm, described in detail below in the 

15 Reconstruction Algorithm Section, to determine the positions of the voxels of the 

magnetized body that produce each set of components of the external RF field associated 
by phase. It is demonstrated in the Uniqueness of Reconstruction Section that a RF 
magnetic field produced by a geometric distribution of dipoles is imique. Therefore, a 
imique spatial distribution of magnetic dipoles, such as those corresponding to the bulk 

20 magnetization of each voxel due to processing nuclei, gives rise to a unique magnetic 
field. Thus, the measured external RF field can provide a unique solution for the spatial 
distribution of magnetic dipoles, i.e., magnetized voxels comprising excited nuclei, in the 
body. 

Referring again to FIGURE la, the spatial locator preferably en^)loys a Fourier 
25 Transform Reconstruction Algorithm, described in detail in the Reconstruction 

Algorithm Section, to map a bulk magnetic moment M corresponding to a particular set 
of components associated by phase of the external RF field onto a spatial location or 
locations in the body. The latter case appHes if more than one location in the body gives 
rise to a particular phase component of the external RF radiation. 
30 As shown in FIGURE 1 a, a preferred embodiment of the invention employs an 

open design magnet, such as a Helmholtz coil design, to allow positioning the array of 
detectors 120 close to the body 1 14 to be imaged. The orientation of the magnetic field 
with respect to the body can be selected to optimize the signal to noise ratio of the 
signals detected by the array of detectors 120. For example, in case of imaging a patient 
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body, the piimary magnetic field can be selected to be coaxial witti die body, or it can 
alternatively be perpendicular to the body axis. 

The 4D-MRI s^paratus of FIGURE la provides a numb^ of advantages. In 
particular, because all of the data is acquired at once rather than over hundreds of ' 
5 repeated sequences of pulses as is the case with NMR systems of the prior art, the 
apparatus of the invoition is particularly suited for imaging cardiopulmonary and 
vascular systems. Further, since the qiparatus may ina:ease the number of detectors to 
increase die resolution, the apparatus of the invention may achieve a hi^er resolution 
such as 10"^ COT* ; th^by, pennitting physicians to view human anatomy and pathology 
10 in a manner not available with conventional imaging techniques. Further, the present 
technique may provide a three dimensional image that can be displayed from any 
perspective. 

An alternative embodiment of the nuclear magnetic resonance (4D-MRI) 
apparatus of the invention, shown in FIGURE lb, employs a three-dimensional array of 

15 detectors 230, spaced apart to detect spatial variations of the emitted RF radiation at least 
at the Nyquist frequency. A magnet 212 provides a magnetizing in a volume to be 
occupied by a body, i.e., image space. A magnetometer 216 measures the magnetizing 
field at a plurality of positions in the image space in the absence of the body to determine 
the uniformity of the primary field. As in embodiment of FIGURE la, the body 214 to 

20 be imaged is placed in a magnetizing field provided by the magnet 212. An RF 

generator 218 in combination with an RF antenna 218a ^ly an RF pulse or a sequence 
of RF pulses to the body to polarize selected nuclei of the body The three-dimensional 
array of detectors 230 provide output signals in response to the RF radiation emitted by 
the body. The RF signal may be recorded synchronously to permit the relative phases of 

25 the Fourier components comprising the RF signal, hi an embodiment, the detection may 
be synchronized relative to the excitation. A digitizer 224 digitizes the output signals 
and sends the digital signals to a construction processor 226 that determines variations of 
the NMR image of the body in a manuCT similar to that described in connection with the 
embodiment of FIGURE la. A display processor 228 receives the information regarding 

30 the spatial variations of the intensity of the NMR signal &om the construction processor 
226, and provides a two-dimensional or a three-dimensional NMR image of the body. 
Employing a three-dimensional array of detectors advantageously decreases the 
acquisition time because the emitted RF signal over the entire three dimensional sample 
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space is detected at once. The shortening oftheacqidation time in turn reduces motion 
artifacts in the NMR image. 

Uniqueness of Reconstruction 

5 

The nature of the RF field can be detennined fiom Maxwell's equations s^plied 
to a sinusoidal current With a sinusoidal current J(x*) confined to small region 
compared with a wavelength, the solution of the vector potential A(x) is [4] 

10 where it = ^ is the wavenumber, and a sinusoidal time dq>radence is understood. The 

c 

magnetic induction is given by 

B = VxA (13) 
while, outside the source, the electric field is 

E = -J-VxB (14) 
k 

15 For a source of dimension rf, the fields in the far zone defined by «X «rare 

transvme to the radius vector and fall oflFas r'\ typical of radiation fields. For tfie near 
zone where r«k (or Ar « 1), the exponential in Eq. (12) can be replaced by unity. 
Then the vector potential is given by 

Tr=i:4-)=-S;7IT^^Jj(»yX(0'.fVV (15) 

At — > u Ci^^zl-tl r * 

20 This shows that the near fields are quasi-stationary, oscillating hannonically as , but 
otherwise static in character. 

Nuclear magnetic resonance (NMR), which is commonly called magnetic 
resonance imaging (MRI), is a means to measure the primary and secondary magnetic 
fields to provide the input to the NMR reconstruction algorithm. The proton 
25 gyromagnetic ratio y ^ I lit is 

I ln^M.^l(m MHz (16) 
The NMR ficquency / is the product of the proton gyromagnetic ratio given by Eq. (16) 
and the magnetic flux B , 

/ = / 2;cfl = 42.57602 MHz T'B (17) 
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A typical flux for a siq)eroonducting NMR imaging magnet is 0. 25 T. According to Eq. 
(17) tins coiresponds to a radio frequency (JRF) of 10.6 Affife vMch corresponds to a 
wavelength of 28.3 iw . In the present invention, each RF antennae of an array is located 
at a distance of about 10 cm from the voxels within tiie image space. Thus, the RF field 
5 is detected in the near zone where r«A (or Ar «1), and the near fields according to 
Eq. (IS) are quasi-stationary, oscillating harmonically as , but otherwise static in 
character. The transverse RF magnetic field of each voxel is that of a RF dipole, (he 
maximum amplitude is given by Eq. (1) wherein the Larmor fipcquency of each voxel is 
determine by the uniform primary field , and , the magnetic moment along the z- 

10 axis, of Eq. (1) corresponds to the bulk magnetization M of each voxel. 

An object containing nuclei with a magnetic moment, herein called a phantom, 
can be considered as a collection of small volxune elements or voxels. When a static 
magnetic field with lines in the direction of the z-axis is applied to an object 
comprising a material containing nuclei such as protons that possess magnetic moments, 

1 5 the field magnetizes the material. As a result a secondary field superposes the applied 
field as shown in FIGURE 9. In the applied magnetic field, the magnetic moments of 
each nuclei precesses about the applied magnetic field. However, the magnetization of 
any one nucleus is not observed from the macroscopic sample. Rather the vector sum of 
the dipole moments from all magnetic nuclei in the sample is observed. This bulk 

20 magnetization is denoted by the vector M . In thermal equilibrium with the primary field 
H^, the bulk magnetization M is parallel to Hq. In an mibodiment of the present 
invention, the magnetization vector is rotated into the transverse plane by an additional 
RF field H, . The magnetization vector tiien comprises a rotating magnetic dipole m in 
the transverse plane. The NMR image may be reconstructed by samplmg the external 

25 field from a series of RF dipoles. 

T he field strength of a magnetic dipole mommt is a fiinction of the external 
position in space relative to the dipole. For convenience of analysis, the field of a series 
of static dipoles m having the coordinates shown in FIGURE 9 is analyzed for 
uniqueness. (The uniqueness of the field of a set of static dipoles £^lies for the 

30 equivalent case of RF dipoles oriented in the transverse plane.) Considering FIGURE 9, 
the net magnetic field at a point extrinsic to the phantom is a sum of the applied field and 
the contributions of each of the voxels within the object, the secondary field. The field is 
point sampled over a three dimensional space. 
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The secondaty magnetic field due to magnetized tissue has to be modeled as 
nonint^racting dipoles aligned with the imposed field. It is demonstrated below that the 
field of any geometric distribution of dipoles is unique, and the siq)erposition principle 
holds for magnetic fields; therefore, a unique spatial distribution of dipoles gives rise to a 
S unique secondary magnetic field, and it is fiirtfaer demonstrated below that this secondary 
field can be used to solve for the NMR image in closed form. It follows that this map is 
a unique solution. To prove that any geometric distribution of dipoles has a unique field, 
it must be demonstrated that the field produced by a dipole can serve as a mathematical 
basis for any distribution of dipoles. This is equivalent to proving that no geometric 

10 distribution of dipoles can produce a field which is identical to the field of a dipole. 

By symmetry considerations, only three distributions of uniform dipoles need to 
be considered. A magnetic dipole has a field that is cyiindrically synmietrical. A ring, a 
shell, a cjdinder, and a sphere of dipoles are the only cases which have this symmetry. A 
cylinder is a linear combination of rings. Thus, the uniqueness of the dipole field is 

15 demonstrated by showing that it is different fiom that of a ring, a shell, and a sphere. 
The uniqueness of the dipole firom the cases of a ring, a shell, and a sphere of dipoles is 
demonstrated in APPENDIX I, APPENDIX H, and Appendix IV, respectively. The plot 
of the three cases of the field of a ring, shell, and a sphere of dipoles each of radius R 
and magnetic moment /w = 10^ Gcm^ givenbyEq. (I.14),Eq. (II.17),andEq. (IV.16)of 

20 APPENDIX I, APPENDIX II, and APPENDIX IV as a fimction of radius R where the 
position of the center of each distribution relative to the detector is the point (0,0,10) is 
given in FIGURES 2, 4, and 6, respectively. Since the fields vary as a fimction of radius 
R , the dipole field is not equivalent to these distributions of dipoles. It is further 
mathematically proven in APPENDIX HI that the field produced by a shell of magnetic 

25 dipoles is different from that of a single dipole. All other fields are a linear combination 
of dipoles. Thus, the dipole is a basis element for the reconstruction of a NMR image. 
Since each dipole to be mapped gives rise to a unique field and since the total field at a 
detector is the superposition of the individual unique dipole fields, linear independence is 
assured; therefore, the NMR image is unique. In other words, there is only one solution 

30 of the NMR image for a given set of detector values which spatially measure the 

superposition of the unique fields of the dipoles. This map can be reconstructed using 
the algorithms described in the Reconstruction Algorithm Section. 

Eq. (1) is a system fimction which gives the magnetic flux output in response to a 
magnetic dipole input at the origin. The phantom is an array of spatially advanced and 
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delayed dipoles weighted according to the bulk magnetization of each voxel; this is the 
input function. The secondary flux is the superposition of spatially advanced and 
delayed flux, according to Eq. (1); this is the output functioiL Thus, the response of 
space to a magnetized phantom is given by the convolution of Eq. (1) with the serieis of 
S weigihted, spatially advanced and delayed dipoles r^resenting the bulk magnetization 
map or NMR image of the phantom. The discrete signals are recorded by a detector 
array over the sample space comprising the xy-plane and the positive z-axis of FIGURE 
9. In an ^bodiment of the present invention, the magnetization vector is rotated into 
the transverse plane by an additional RF field . The magnetization vector then 
10 comprises a rotating magnetic dipole m in the transverse plane. The NMR image may 
be reconstructed by sampling the external field firom a series of RF dipoles. The discrete 
signals are recorded by a detector array over the sample space comprising the xz-plane 
and the positive y-axis of FIGURE 11, 

15 and Jj NMR Images 

The NMR active nuclei including protons posses both angular momentum and a 
magnetic moment. When nuclei are placed in a static magnetic field Hp , they process 
about the field at a fi:equency proportional to the magnitude of Hq . The bulk 

20 magnetization M of each voxel conq>rises the vector sum of the magnetic moments &om 
all of the nuclei in each voxel. If the precessing nuclei are then subjected to an additional 
rotating (RF) field , which is synchronous with the precession, their magnetic 
moments and thus M will precess about H, and rotate away bom the primary field Hg 
by an angle in a coordinate fi[ame which rotates at the Laimor firequ^cy [3]. The 

25 precession about H, continues as long as Hj exists. The final value of then depends 
on the strength of H, , which determines the precession rate, and the time for which it is 
turned on. The nucld absorb energy as they change their orientatioiL This is known as 
nuclear magnetic resonance (NMR). The temperature of the nuclei or nuclear spin 
system rises dining absorption of energy. When the Hj field is removed, the spin 

30 system cools down until it is thermal equiUbrium with its environment The exponential 
relaxation of the spin system temperature to tiiat of the surrounding lattice is called spin- 
lattice relaxation and has a time constant where a time constant is defined as the time 



wo 02/16956 PCT/USOl/25954 

25 

it takes for 63% of the relaxation to occur. The NMR signal may also decay because the 
nuclei initially in phase following the H| pulse get out of aligoment with each oth^ or 
dq)hase by local interactions with the magnetic fields of neighbor nuclei. The dephasing 
of the NMR signal is due to differing precession rates efiTected by the local interactions 
5 and is described by an exponential time constant Tj also known as spin/spin relaxation. 
The main source of NMR image ( also called magnetic resonance images (MRI)) 
contrast is T, and v/tach depend on tissue types. 

In an embodiment of the present invention, a image is produced by a q)plying 
at least one pulse sequence that inverts the magnetization and records the relaxation, a 

10 technique called inversion recovery. For example, the RF receivers are switched on to 
follow the decay following the nuclear excitation comprising a H, pulse. MX0» 
time dq)radent bulk magnetization m the direction of the primary field (coordinates 
of FIGURE 1 1) is examined at a time after an inverting pulse by applying another H| 
pulse equivalent to a rotation by 90*^ after waiting the time following the initial 

IS inversioiL The 90^ pulse puts the z magnetization Af^C^o) into the transverse plane for 
observation. Changing the waiting time allows for observation of Af,(0 at different 
times during relaxation. 

In an embodiment of the present invention, a image is produced by a applying 
at least one pulse sequence that flips the magnetization vector into the transverse plane 

20 and records the transverse relaxation by producing at least one spin-echo. The image 
depends on the NMR signal decaying because the nuclei initially in phase following the 
pulse get out of alignment with each other or dephase by local interactions with the 
magnetic fields of neighbor nuclei. The dephasing of the NMR signal is due to differing 
precession rates effected by the local interactions and is described by an exponential time 

25 constant 7, . Another unwanted source of dq)hasmg is due to inhomogendties in the 
primary field across the image space that causes an addition contribution to the 
precession rates of the magnetic nuclei. The signal ones observes after flipping the 
magnetization into the transverse plane that includes dephasing &om an inhomogeneous 
Hg field is still known as the firee induction decay (FID), but the is called T2 relaxatioiL 

30 The relaxation may be recovered firom a T^* FID since the inhomogeneity is 
constant and may be reversed. The dephasing due to the static inhomogeneity of 
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may be canceled out by applying a 180** Hj pulse at time along the yR-axis as shown 
in FIGURE 1 1 where the relaxation occurs in time . After an additional time t^^ 
the total time elapsed after the 90^ Hi pulse is It^ = which is the spin-echo time. At 
this time the dq)hasing due to the static inhomogeneity of Hq is exactly canceled out; so, 
5 the relaxation is strictly due to diose processes that create 7, relaxation which is 
recorded. After each 180'' pulse, another spin echo is formed. The envelop of the 
maximum ampUtude of the spin echoes is the relaxation. la an embodiment, a pulse 
sequence to give the data for a image known as the Carr-Purcell-Meiboom-Gill 
(CPMG) sequence [3] conq>rises applying a 90® pulse along the Xj^-axis followed by a 
10 series of 1 SO"* pulses along the yj^ -axis at times + Ini^ where n is an integer including 
zero. 

In other words, when nuclei are placed in a static magnetic field and then 
subjected to an additional rotating (RF) field H, , which is synchronous with their 
precession, M will precess about H, and rotate away from the primary field by an 
IS angle . The magnitude of M is a maximum initially and decays with time. This 

occurs by emission of the same multipolarity radiation that it absorbed and by transfer of 
energy to the surrounding lattice. The intensity of the radiation is a function of M and 
tiie coordinate position relative to the RF emitting voxel. In the present invention, the 
measurement of the intensity of the RF signal is performed over time and space 
following Ti and/or encoding pulses. The signal as a function of time at a given 
detector position is Fourier transformed to give the components each having an 
amplitude and a unique phase. Each set of components of the NMR signal over the 
sample space due to a given voxel is determined from the phase data and the detector 
positions. Hie location of each voxel is determined through the spatial variations of the 
intensity of the transverse NMR field of the set of components associated by phase. 

4D-MRI System 

An NMR apparatus used to generate, and measure the secondary field and 
reconstruct the image is shown in FIGURE 8, and the corresponding coordinate system 
is shown in FIGURE 10. The apparatus comprises 1.) a magnet including a 
superconducting magnet to magnetize a volume of an object or tissue to be unaged, 2.) a 



wo 02/16956 PCTAJSOl/25954 

27 

means (magnetometer) to detennine the primaiy or magnetizing flux over the image 
space in the abs^ce of the object to be imaged. In one embodiment, NMR of a proton 
containing homogeneous phantom such as water is determined on a pomt by point basis 
to map the primary field, or the magnetic flux at multiple points is obtained 
5 simultaneously using the reconstruction algorithm described herein, 3.) a radio frequency 
• (BF) generator and transmitt^ inchidii^ an antennae such as saddle coils to excite the 
protons of the magnetized volume^ 4.) a means including an antennae coil to sample the 
dipole component (z-con^nent in terms of Eq. (1)) of the RF secondary magnetic field 
at the Nyquist rate in time oveac the proton fi^e induction decay, 5.) a detector array of 
10 elements of the means to sample &e time signals including an antennae array which is 
selectively responsive to the dipole component (z-component of Eq. (1)) of the RF 
magnetic field of the magnetic moments of the protons which are aligned along the 
transverse axis to ideally point sample the secondary magnetic field at the Nyquist rate 
over the spzAdl dimensions which uniquely detennine the NMR image which is 
1 S reconstructed from ttie measurements, 6.) an analog to digital converts: to digitize the RF 
signals, 7.) a time Fourier transform processor to convert the signal at each detector over 
time into its Fourio* components, 8.) a processor to associate the Fourier components due 
to each voxel by phase, 9.) a reconstruction algorithm processor including a Fourier 
transform processor to convert each set of components into a voxel location of the bulk 
magnetization in the image space in parallel over all the voxels to form an NMR image, 
and 10.) image processors and a display such that the NMR image can be rotated in 
space to be displayed from any perspective as a three dimensional or two dimensional 
(tomographic) image. 

Magnetizing Field 

In an embodiment of the present invention, the applied magnetizing field which 
permeates the object to be imaged including tissue is confined only to that region which 
is to be imaged The confined field limits the source of signal only to the volimie of 
interest; thus, the volimie to be reconstructed is limited to the magnetized volmne which 
sets a limit to the computation required, and eliminates end effects of signal originating 
outside of the edges of the detector array. In the NMR case, the field having a steep 
gradient at the edges Umits the imaged region by providing a range of Lannor 
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frequencies i^erein the data is comprises a narrow frequency band at a desired Larmor 
frequency* 

Detector Array 

5 

An embodiment of the NMR imager of the present inv^tion comprises a detector 
array of multiple detector elements which are arranged in a plane. The array may be a 
two dimensional detector array which is translated over flie third dimension during the 
scan, or it may be a three dimensional detector array. The individual detectors of the 
10 array may respond to a single component of the secondary magnetic field which is 
produced by the magnetized object including tissue where the conyjonent of the field to 
which the detector is responsive determines the geometric system function which is used 
in the reconstruction algorithm discussed in the Reconstruction Algorithm Section. The 
detectors ideally point sample the secondary magnetic field at the Nyquist rate over the 
1 5 spatial dimensions which uniquely determine the NMR image which is reconstructed 
from the measurements. 

Small antennas may measure the RF signals as point samples without significant 
decrease in the signal to noise ratio relative to large antennas by using impedance 
matching while minimizing resistive losses by using superconductmg reactance 
20 elements, for example. In an embodiment, cross talk between antennas is ameliorated or 
eliminated by time multiplexing the signal detection over the array of antennas. The RF 
field may be sampled synchronously or the field may be sampled at known times so that 
the phase at any given detector may be related to that at any other detector. 

Micromagnetic field sensors that are used to detect the primary field in the 
25 absence of the object to be scanned include NMR detectors and superconducting 

quantum interference devices (SQUIDS), Additional devices have been developed that 
are based on galvanometric effects due to the Lorentz force on charge carrieis. In 
specific device configurations and operating conditions, the various galvanomagnetic 
effects (Hall voltage, Lorentz deflection, magnetoresistive, and magnetoconcentration) 
30 emerge. Semiconductor magnetic field sensors include MAGFETs, magnetotransitors 
(MT), Van der Pauw devices, integrated bulk Hall devices including flie vertical MT 
(VMT), and the lateral MT (LMT), silicon on sapphire (SOS) and CMOS 
magnetodiodes, tiie magnetounijmiction transistor (MUJT), and the carrier domain 
magnetometer (CDNQ, magnetic avalanche transistors (MAT), optoelectronic magnetic 
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field seiisois, and iiiagaetoiesistive magnetic field sei^ InfhecaseofNMR 
measuremeiit of (he secondary field (and/or fhe primary field), the detector array 
comprises RF antennas described in the NMR Primary Magnet, Gradient Magnets, RF 
Generator, RF Transmitter, and RF Receive Section. 

5 

Scanning Methods 

The NMR scan perfbnned on the object to be imaged including a human 
comprises the following steps: 

10 

« The magnetic moments of nuclei including protons of the object to be imaged that are 
aligned by the primary field are furth^ aligned by a radio frequency (RF) pulse or 
series of pulses. 

1 5 • The strength and duration of the rotating (RF) field that is resonant with the protons 
of the magnetized volume and is oriented perpendicularly to the direction of the 
magnetizing field is applied such that the final precession angle of the magnetization 
is 90 (^^ = 90° ) such that the RF dipole is transverse to the primary magnetizing 

field and peq)endicular to the RF magnetic field detector. 

20 

• NMR pulse sequences which provide the signals for a 2J or 7^ image may be applied 

For example a 90 ® pulse may be followed by a series of 1 80 ® pulses. One sequence 
is the Cair-Purcell-Meiboom-Gill (CPMG) sequace [3]. 

25 • The free induction decay signals are recorded. 

" The time dependent signals are Fourier transformed to give the intensity and phase of 
each component. The NMR signal of each voxel at any given detector gives rise to a 
Fourier component with a unique phase angle relative to fiie Fourier component of 
30 any other voxel of the phantom at that detector. 

• The matrix of Fourier components that correspond to the NMR signal of a given voxel 

over the detectors is detennined. This may be achieved by using a first component 



wo 02/16956 PCTAJSOl/25954 

30 

having a phase angle and calculating the phase angle as a function of spatial position 
of the first detector relative to any other detector and identi^^g the oonqranent at 
each detector having the calculated phase angle. The matrices are deteimined for all 
of tiie voxels. The measurements of the spatial variations of ttie transverse RF field 
S of a given matrix is used to determine the coordinate location of each voxel Thus, 

each matrix of components associated by phase comprises the intensity variation over 
the sample space of the RF field of the bulk magnetization M of each voxel. 

• The Fourier transform algorithm given in the Fourier Transfonn Reconstraction 

0 Algorithm Section is performed on each set of components over the detector array to 
map each bulk magnetization M corresponding to a voxel to a spatial location over 
the image space. 

• In one embodiment, the point spread of the reconstructed voxel is corrected by 

5 assigning one voxel above a certain threshold with the bulk magnetization M . The 
other voxels are assigned a zero value. In this case, the reconstruction is digital 
versus analog. In terms of the signal to noise ratio, the superiority of digital over 
analog is generally known to those skilled in the art of signal processing. 

0 • In the limit with sufficient phase resolution, each volume element is reconstructed 

independently in parallel with all other volume elements such that the scan time is no 
greater than the nuclear fi-ee induction decay (FID) time. 

NMR Primary Magnet, Gradient Magnets, RF Generator, RF Transmitter, and RF 

5 

Receiver 

The primary magnet is that typical of a magnetic resonance hnaging apparatus. 
The primary magnetizing field is in the z-direction as shown in FIGURE 10, usually 
0 coaxial with the patient's body. However, in an another embodiment the primary 
magnetizing field is in the direction which is perpendicular to the patient's body. An 
open primary magnet design is preferable such as a Helmholtz coil design to 
accommodate the three dimensional RF antennae array. A fiirther field is applied to have 
a gradient in the z-direction (FIGURE 10). This provides a unique field in a chosen 
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volumeofthepatiffltincIuduigacross-^tioiialsUceoffhepa^ In another 
embodiment, ttie NMR (Laimor) fiequ^y fiom each voxel is detennined by the 
magnitude of flie magnetizing flux at the voxel position. The primary magnetic field has 
a sharp gradient at the edges of the image space whorein the Larmor fi^equencies outside 
S of a selected range are rejected as arising fiom voxels outside of the image volume. 

A rotating RF field Hj , of frequency chosen to cause resonance in the slice is 
thenapphed. Thus, only the molecules in the slice resonate. The resonance signal from 
the slice can then be detected. NMR pulse sequences which provide the signals for a ?J 
or image may be ^plied. For example a 90 pulse may be followed by a series of 
10 180 pulses. One sequence is the Carr-Puicell-Meiboom-Gill (CFMG) sequence [3]. 

An embodiment of the NMR imager of flie present invention comprises a RF 
generator 14 and 18 and RF transmitter coils 9 shown in FIGURE 8 which provide the 
rotating H, (RF) field that is resonant with the protons of the magnetized volume and is 
oriented perpendicularly to the direction of the magnetizing field as shown in FIGURE 
15 11. When the precessing nuclei are subjected to the additional rotating (RF) field Hj, 
which is synchronous with the precession, their magnetic moments and thus M 
precesses about H , and rotate away from the primary field Ho by an angle in a 
coordinate fi^e which rotates at the Larmor frequency [3]. The precession about Hj 
continues as long as H, exists. The final value of then depends on the strength of 
20 H, , which determines the precession rate, and the time for which it is tumed on. In an 
embodiment, the strength and duration of H, is such that = 90** such that flie dipole 
is oriented in the x^y^ -plane of FIGURE 11. A rotating coordinate firame is 
traditionally used to explain the physics of NMR [3]. Thus, is terms of the traditional 
NMR coordinate designation as described by Patz [3] and Hounsfield [7] shown in 
25 FIGURE 11, the Xr -axis and y^-axis rotate about the primary field Ho oriented parallel 
to the -axis at the Larmor frequency relative to the stationary NMR coordinate system 
(x,y,z) shown in FIGURE 10. Thus, the additional rotating (RF) field H, and the 
transverse RF magnetic field are stationary m the rotating NMR coordinate system, but 
rotate at the Larmor frequency in the stationary NMR coorxlinate system. Both the 
30 rotating and stationary coordinate systems are shown in FIGURE 11. In an embodiment, 
the rotating H, (RF) field is along the x^-axis and rotates the magnetization vector by 
the angle 90** to into the x^y^ -plane. In terms of the reconstruction by the present 
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invention, the system function of the field conesponding to the rotated magnetization 
vector is equivalent to that of the detection of the z-component of a z-oriented dipole. 
The y-axis is the unique axis offhe NMR system shown in FIGURE 8 and flie stationary 
NMR coordinates shown in FIGURES 10 and 11. This axis corresponds to the unique 

S axis of the Fourier Transform Reconstruction Algoriflmi Section, the z-axis shown in 
FIGURE 9, and the magnetization axis is the z-axis. The magnetization axis of the NMR 
rotating cooidinates [3,7] is the Zi^-axis shown in FIGURE 11. The measured transv^e 
RF magnetic dipole oriented in the xy-plane is periodically parallel to the y-axis and 
rotates about the primary field H^, oriented parallel to the z-axis at the Lannor firequency, 

0 The RF transmitter includes saddle coils. The RF receiver is a three dimensional 

array, or two arrays 180® fix)m each other, or two orthogonal pairs of arrays, wherein 
each member of a pair is 180® from each other. The detector ideally point samples the 
RF field at the Nyquist rate described in the Reconstruction Algorithm Section. In one 
embodiment shown in FIGURE 1 1, each detector 20a is a coil antennae perpendicular to 

5 the y-axis and selective to the y-component of the RF secondary magnetic field arising 
from each voxel 14a magnetized in the z-axis direction. This two dimensional array is 
translated along the y-axis during a scan where readings of the secondary magnetic field 
are obtained as a fimction of the translation. In another embodiment shown in FIGURE 
12, the array 401 is three dimensional comprising multiple parallel two dimensional 

[) arrays 402 wherein each two dimensional array has a plurality of antennae coils 403 
which detect the RF field fix)m each magnetized voxel 404. 

In another embodiment, the rotating H| (RF) field shown in FIGURE 1 1 is along 
the -axis and rotates the magnetization vector by the angle = 90® to project into 
the x^yR -plane. In t^ms of the reconstiuction by the present invoition, the system 

5 fimction of the field corresponding to the rotated magnetization vector is equivalent to 
that of the detection of the z-component of a z-oriented dipole. The x-axis is the unique 
axis of the NMR system shown in FIGURE 8 and the stationary NMR coordinates 
shown in FIGURES 10 and 1 1. This axis corresponds to the unique axis of the Fourier 
Transform Reconstruction Algorithm Section, the z-axis shown in FIGURE 9, and the 

) magnetization axis is the z-axis. The magnetization axis of the NMR rotating 
coordinates [3,7] is the -axis shown m FIGURE 1 1 . The detector array receives 
similar signals as tiiose produced by the transmitter, and both are perpendicular to the x- 
axis. Each detector is a coil antennae peip^dicular to the x-axis and selective to the x- 
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component of fbe RF secondary magnetic field arising fix>m voxels magnetized in the z- 
axis direction. This two dimensional array is translated along the x-axis during a scan 
where readings of the secondary magnetic field are obtained as a fimction of the 
translation. In another onbodiment, the array is three dimensional comprising multiple 
5 parallel two dimensional arrays 402 as shown in FIGURE 12 with the exception that the 
array 401 is perpmdicular to the z-axis. The measured transvo^e RF magnetic dipole 
oriented in the xy-plane is periodically parallel to the x-axis and rotates about the 
primary field oriented paraOel to the z-axis at flie Lannor firequency. 

A linear combination of the cases of rotation of the RF magnetization along the 
10 x^-axis and the yR-axis is within the scope of the present invention and may be adopted 
in a manner straightforward to those skilled in the NMR art to apply the fields and 
detectors described for this invention. 

The method described herein may be performed on a suitable NMR examining 
apparatus such as that shown in simplified foim in FIGURE 8, and the corresponding 
15 coordinate system is shown in FIGURE 10. Illustrated schematically only are coils 6, 
which provide , the steady primary field 7, which provide , the field gradient in the 
X-axis direction as shown, 8 which provide the , the field gradient in the y-axis 
direction as shown, 9 vMch provide the RF field, and 10, which provide G„ the field 
gradient in the z-axis direction as shown. The coils are driven by i(, , G,, (j^, RF, and 
20 G, drive amplifiers 11,12, 13, 14, and 15 respectively, controlled by ^ , G^, RF, and 
control circuits 16, 17, 18, and 19, respectively. These circuits can take suitable 
forms which will be well known to those with experience of NMR equipment and other 
apparatus using coil induced magnetic fields. The circuits are controlled by a central 
processing and control unit 20 to achieve the desired primary field, field gradients, and 
25 RF field to rotate the magnetization vector such that it is perpendicular to coil 9. 

The RF coils may be two saddle shaped coils 9 which are driven in parallel to 
provide the rotating RF field. The FID signals sensed are received in this example by the 
three dim^ional array of RF coils 30 and are amplified by an RF amplifier 21 before 
being applied to signal handling circuits 22. The three dimensional detector array 30 is 
30 shown in more detail in FIGURE 12 as the three dimensional detector array 401 . The 
circuits 22 are arranged to make any appropriate calibrations and corrections but 
essentially transmit the signals to the processing circuits to provide the required 
representation of the examined volume. These circuits can conveniently be combined 
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wiih the circuits wbidi control flie primary field, field gradioits, and RF field and thus 
are included in the circuits indicated at 20. The picture thus obtained is viewed on a 
display 23, such as a television monitor, and this may include uoputs and other . 
pmpherals 24 for the provision of commands and instructions to the machine, or other 
S forms of ou^ut The display is not limited and includes any medium of conveing the 
image. Exanqples of displays, but not limited to, include printers, cathode ray tube 
di^lays, liquid crystal di^lays, plasma screens, three dimensional modelers, 
holographic displays, laser monitors, and projection monitors. 

The ^aratus also includes field measurenient and error signal circuits 25 which 
10 receive signals via anq>Iifiers 26 fiiom field probes X^, X^, Y^,md shown. 

The patient 27 is inserted in the tubular former of and Gy coils 7, 8 and is 

supported there by a suitable couch or other supporting means. Such supports may be 
readily provided in any suitable foinL 

The coils 7, 8 are two sets of coils axially displaced, each set comprising two 

15 pairs of saddle coils the pair 7 being at 90** to the pair 8. These coils are themselves 
inserted into the central aperture in coils 6 which in an embodiment are wound in 
four parts connected in series to provide an approximately circular configuration which 
is well known to be desirable for production of a uniform field. Further details of the • 
coil winding will not be given since suitable coils can readily be devised, by those with 

20 the appropriate skills, to provide the fields required. 

The appropriate stores provide the amplitude and duration signals which are 
converted to analog form in digital to analog converters (D ACs) and applied to 
respective coil drive circuits x, y, z, RF. The respective drive circuits, which can take 
any form well known for driving field coils, provide the specified currents to the 

25 appropriate coil for the specified duration. The apparatus and circuits described so far 
may be adopted to provide different gradients and RF fields, by appropriately adjusting 
the stored sequences and profile data. Similarly other known NMR apparatus which are 
capable of applying a steady magnetic field, a pulsed RF field, and , , and G, field 
gradients to a body, may be adopted in a manner strai^tforward to those skilled in the 

30 NMR art to apply the fields described for this invention. 

Gradients may be s^pUed in any direction to fiirther enhance the image 
reconstruction or image quality by methods known to those skilled in the art In an 
^bodiment shown in FIGURE 8, the system components which provide imaging 
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enhancing gradients are 7, ^ch provide the field gradient in the x-axis direction as 
shown, 8 which provide the , the field gradimt in the y-axis direction as shown, and 
10, which provide G„ the field gradient in the z-axis direction as shown. For example, a 
magnetic field gradient along an axis in tiie xy-plane of the z-directed primary field 
5 is applied to produce a phase dispersion as a fimction of flie distance along the axis, hi 
an embodiment, the phase due to the gradient may be linear and may be larger than the 
phase of any Fourier conqK)nent in the absence of a gradient During reconstruction, the 
phase mformation due to the gradient may be applied to refine the assignment of the 
position of each voxel with respect to the distance along the gradient axis by methods 

10 known to those skilled in the art 

Wifli the basic signal handling system of the present mvention, the FID signals 
fiom the signal sensing coils of the detector array 30 (shown in more detail in FIGURE 
12 as coils 403 of &e three dimensional detector array 401) are amplified in an RF 
amplifi^ and applied via an analog to digital converter (ADQ to a store such as a 

1 S random access memory (RAM). The data is tiien processed according to the procedure 
given in the Reconstruction Algorithm Section. 



Phase Angle and Associated Set of Fourier Components 

20 

In an enibodim^t, the NMR signal at each detector as a function of time is 
processed by a method such as a Fourier transform operation to give a plurality of 
Fourier components each having an intensity, a phase angle, and the same firequency. 
The NMR signal of each voxel at any given detector gives rise to a Fourier component 

25 with a unique phase angle relative to the Fourier component of any other voxel of the 
phantom at that detector. The set of Fourier components that correspond to the NMR 
signal of a given voxel over the detectors is determined. This may be achieved by using 
a component witii a defined standard phase relative to a first detector. Over the array of 
sampled points, the phase angle of components at each detector may be converted to the 

30 corresponding phase at the position of the component having the standard phase. The 
components that have the standard phase are associated into a set of components 
comprising the components &om a given voxel over the sample space. 
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In an embodiment, the phases of the components from the voxel are determined 
at the other detectors using the relative position of other detectors of the sample space 
relative to ttie first detector and the standard phase of the component at the fiist detector. 
In other words, the set of associated components may also be detennined by 1 .) 
S identifying a first component having a characterizable phase angle, 2.) calculating the 
phase angle as a iimction of spatial position of the other detectors relative to the first 
using the phase angle of the conq)onent at the first detector, and 3.) identifying the 
component at each detector having the calculated phase angle. These components 
identified fi:om the data are associated to form the set 

10 In the present invention, the phase variation of the NMR signal is provided by the 

combination of 1 .) the angle 9 suspended between the direction of the detector and the 
tadial vector, the vector fix)m the dipole to the detector, and 2.) the angle ^ due to a 
sq>aration distance r between a voxel and a detector given by the wavenumber of the RF 
field * times r . The distance r of the detector fiiom the voxel gives rise to the phase 

15 tenn e^^ of the component of the detected RF signal where k is the wavenumber of the 
NMR signal. The harmonic oscillation of each RF dipole is equivalent to the dipole 
rotating in the transverse plane. The detector is responsive to a component in this plane. 
At a point in time, each RF dipole is directed at an angle 0 relative to the direction of 
detection of the detector. The angle 6 of the RF dipole relative to the direction of 

20 detection axis of the detector gives rise to a phase angle term e"'^ . In a preferred 

embodiment, the simi of the phase angles, kr and 6 , are unique for each voxel at each 
detector. The position of each detector relative to a different detector may be used to 
calculated the phase angle of the second relative to the first. This may be repeated over 
all of the detectors to give the set of intensities of the NMR signal over the sample space 

25 due to each voxel. 

In an embodiment, the voxels that are on the same axis with a given detector or a 
plurality of detectors that align vertically with one or more voxels is determined. The 
phase angle is then only a fimction of the distance r of the detector firom the voxel. Such 
components may be identified by the presence of at least one component of a phase 

30 given by ^ + kAr where is the phase angle of the first component and Ar is the 
distance of the second detector relative to the first of the detector. This detector(s) is 
defined as having a zero angle 0 . In this case, the phase angle ^ of a second component 
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at a second detector aligned on the detection axis coiresponding to a first conq}on^t at a 
first detector having a phase angle ^, is given by 



5 torn the voxel are detennined at the other detectors using the relative position of other 
detectors of the sample space relative to the first detector and of the component at the 
first detector. 

FIGURE 14 shows the stationary coordinate system (x, y, z) of the NMR detector 
corresponding to the coordinate system of FIGURE 9 and FIGURE 10 of a primary field 

10 H„ oriented parallel to the z^-axis and the z-axis and the corresponding transverse RF 
magnetic dipole oriaited in the x^y^ -plane and periodically parallel to the y-axis 
wh^in the spatial variation of the RF y-component of the flux due to the RF dipole is 
measured at a detector according to Eq. (1). FIGURE 14 fiirther shows the distances and 
angles between a voxel 801 linear to a first detector 802, a second nonlinear voxel 800, 

1 5 and a second detector 803 in accordance with the inventioit The NMR signal of eadi 
voxel at any given detector gives rise to a Fourier component with a unique phase angle 
relative to the Fourier component of any other voxel of the phantom at that detector. In 
the present invention, the imique phase variation of the NMR signal is provided by the 
combination of 1 .) the angle 0 susp^ded between the direction of the detector and the 

20 radial vector, the vector &om the dipole to the detector, and 2.) the angle 0 due to a 
separation distance r between a voxel and a detector. The phase ^ due to a separation 
distance between a voxel and a detector of r is given by the wavenumb^ of the RF field 
k times r. 



25 where the wavenumber k is given by 




In terms of the NMR coordinates of the detector shown in FIGURE 1 1 and FIGURE 14, 
the phase angel ^ is given by 



^=^,+Mr (18) 
From the components identified as coaxial with a detector, the phases of the components 



(19) 




(21) 
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where Ax , Ay , and Az is the separation distance betwem a voxel and a detector in the 
X, j;,and z directions, respectively. In the case that a two dimensional slice is imaged 
in the transv^e plane, the phase angel ^ is given by 

*=*^ = YV(^/ + (A>')' (22) 

5 All voxels which have the same radial distance firom a given detector will have 

the same phase given by Eq. (19). Howcvct, at each case that the radial distances r of 
two voxels to a given detector are equivalent, a unique angle 0 suspended betwem the 
direction of the detector and the vector along r exists. In an embodiment, the RF dipoles 
of all voxels are time synchronous. The NMR pulsing and detection may be 

10 synchronized. At least one standard dipole may be used as reference phase to set the 
phases of the RF dipoles of the phantom. All voxels that are not collinear with the axis 
of a detector will posses an angle Q suspended betwem the direction of the detector and 
the radial vector, the vector from the dipole to the d^ector. Since the RF dipoles are 
time synchronous, at each point in time, this angle corresponds to a phase angle Q of the 

15 RF signal from each voxel at each detector. In the case that a two dimensional slice is 
imaged in the transverse plane, the phase angle Q is given by 

The total phase angle ^j. of the RF signal from each voxel at each detector is given as the 
sum of ^ and 6 . 

20 0j-=^ + 0 (24) 

The uniqueness of the phase angle a Fourier component of the NMR signal of 
each voxel at any given detector relative to a Fourier component of any other voxel of 
the phantom at that detector is demonstrated by the following cases according to tiie 
coordinates shown in FIGURE 14: 

25 

Case 1 Ajc = 0.1 cm Ay=0.1ciw 

0=sin'^-. f'^ 4>^~J(0lf + (0if 

V(0.iy+(0.iy ^ X ^^" '^ • ^"'^^ 

e^O.nrad il>^ 5.9X10-^ rod 

30 = 0 + ^ = 0.78 + 5.9 10"* rod = 0.7806 rad 
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Case 2 Ax = 0.1 cm 6y=2Qcm 

V(0.1)' +(20/ ^ X""" ^ ^ 

5 e =5 10"' rod * =8.37;iriO"* rod 

=0 + 0 =5 JriO"' + 8.377 JTlO-* rod = S.SSJiriO-^ rod 

Case 3 Ax =20 cm Ay = 0.1 cm 

10 0=sin-'-j == A = — V(20y+(0.lf 

V(2oy+(o.i/ ^ ^ 

9 =1.56 rorf 0 = 8.37XlO~^ra</ 

^j. =0 + 0 =1,56 + 8.377^-10"* ro/ = 1.643 rarf 

Case4 Ax =20 cm A^/s 20 cm 

15 

V(20)' + (20y A ^ ^ 

0 = 0.785 rarf ^ =0.1184 rarf 

0r =^ + 0 = 0-785 + 0.1 184 rarf = 0.9034 rarf 

The condition for uniqueness of the phase ^j. of each voxel at a given detector 
20 dqjoids on the unique sum given by Eq. (24). Since all voxels which have the same 
radial distance from a given detector will have the same phase ^ given by Eq. (19), but a 
unique angle 0 , each voxel vdU have a unique phase ^j. when the following condition is 
satisfied: 

0^, the maximum angle 0, is greater ttian , the maximum angle 0, for any 

25 r. 

Q»«>*o»x (25) 
^mn corresponds to and , the minimum separation of the voxel and the detector 
in the X and y directions, respectively, and is given by Eq. (23). 
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corresponds to , the minimmn sq[)aratioii of the voxel and the detector in the x 
direction and , the maximum separation of the voxel and the detector in the 
direction, and is given by Eq. (22). 

5 ^Y'^i^y -^(ynj' (27) 

Substitution of Eq. (26) and Eq. (27) into £q. (25) gives the condition for uniqueness of 
the phase angle 



sm 



"^T=7>T^i^y*(y^y (28) 



10 For » , Eq. (29) gives 

For = J'mi. . Eq.(30) gives 

A>27rV2>u (31) 
A>8.9y^ (32) 
15 Jn tenns of frequency /, Eq. (32) gives 

^'■^ 

where c is the speed of li^t. For ^20 cm, the condition for uniqueness of the 
phase angle is / =< 168 MHz, 

Consider a component of a voxel at a first detector with a phase ^^ defined as a 

20 standard phase. The corresponding phases of the components from the same voxel are 
determined at other detectors using the relative position of the other detectors of the 
sample space relative to first detector. In the case that a voxel is collinear with the first 
detector, the phase angle 0=0. The phase angle ^ of the corresponding component at 
a second collinear detector can be calculated from the phase angle ^, of a first using the 

25 detector separation along the collinear axis, hi the case that a two dimensional slice is 
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imaged and the second detector is displaced Ay in the >^ direction of tiie sanq)Ie space 
relative to the first detector with a phase angle the angle ^ is given by 

^=V(Wy + (*Ayy =V(W)'+^' (34) 
Consider the phase angle ^„ of a corresponding component at a second detector 
S that is nonlinear with the first detector. The phase angle 0„ of the second component 
due to the same voxel may be calculated fix>m In this case a contribution exists &om 
1 .) fhe angle 0 suspended between the direction of the detector and the radial vector, the 
vector firom the dipole to fhe detector, and 2.) the angle ^ due to a separation distance r 
between a voxel and a detector given by fhe wavenumber of the RF field k times r . In 
10 fhe case that a two dimensional slice is imaged and the second detector is displaced Ax* 
and Ay in fhe jc and y directions of fhe sample space, respectively, relative to the first 
detector with a phase angle 0p the phase 9 is given by 



0 =sm ' — =sm i = sm i ■ 

^ V(Ax'y+(Ayy 4{kt^)'^{wy 



= sm 



(35) 



VOrAx-y + CAAyy + ^f 
The angle ^ is given by 

15 ^ = + + {k^yf = V(*^)' + + (36) 

The phase angle 0„ of the second component due to the same voxel is given by the sum 
ofEq.(35)andEq.(36), 

^,,=fl-hfe=sin-'- . ^ ^4Q^)"^M^^', (37) 

The component at the second detector which has phase angle is associated with the 
20 component from which ^„ was calculated. All phase angles ^^i detectors 
are determined firom each standard phase angle ^, . This is repeated for all using those 
detectors which are determined to are collinear with a voxel to give the sets of 
componrats each comprising fhe spatial variation of the RF field over fhe sample space 
due to each voxel. 

25 
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When a static magnetic iSeld with lines in the direction of the z-axis is applied 
to an object comprising a material containing nuclei such as protons that possess 
5 magnetic moments, the field magnetizes the material. As a result a secondary field 
superposes the applied field as shown in FIGURE 9. Jn the applied magnetic field, the 
magnetic moments of each nuclei processes about the applied magnetic field However, 
the magnetization of any one nucleus is not observed fix>m the macroscopic sample. 
Rather the vector sum of the dipole moments fix)m all magnetic nuclei in the sample is 

10 observed. This bulk magnetization is denoted by the vector M. In themial equilibrium 
with the primary field Hj,, the bulk magnetization M is parallel to Ho . The volume to 
be imaged is divided into volume elements called voxels, and the magnetized voxel 302 
shown in FIGURE 9 with bulk magnetization M is modeled as a magnetic dipole m . 
Consider the case wherein data comprising the z-component of the magnetic field of a 

15 dipole oriented in the z-direction is acquired by detectors 301 in the three dimensional 
sample space comprising the xy-plane and the positive z-axis as shown in FIGURE 9. 
The magnetic moment, , of each voxel within the phantom is a magnetic dipole. And, 
the phantom can be considered to be a three-dimensional array of magnetic dipoles. At 
any point extrinsic to the phantom, the z-component of the secondary flux, B , fi^om any 

20 single voxel is 

2z^-x^-y 

8 = '"x , 2 ^ 2 ^ 2vS/2 (38) 

(x +y +z ) 

where x, y, and z are the distances &om the center of the voxel to the sampling point It 
is shown in APPENDICES I-FV that no geometric distribution of magnetic dipoles can 
give rise to Eq. (I). Therefore, the flux of each magnetic dipole (voxel contribution) 

25 serves as a basis element to form a unique reconstruction of the aiiay of dipoles which 
comprise the bulk magnetization mzp or NMR image of the phantom. 

Eq. (1) is a system fimction which gives the magnetic flux output in response to a 
magnetic dipole input at the origin. The phantom is an array of spatially advanced and 
delayed dipoles weighted according to the buDc magnetization of each voxel; this is the 

30 input function. The secondary flux is the superposition of spatially advanced and 
delayed flux, according to Eq. (1); this is the output fimction. Thus, the response of 
space to a magnetized phantom is given by the convolution of Eq. (1) with the series of 
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weighted, spatially advanced and delayed dipoles rq>iesenting the bulk magnetization 
map of the phantom. The discrete signals are recorded by a detector array over the 
sample space comprising the xy-plane and the positive z*axis. 

In an embodiment of the preset invention, the magnetization vector is rotated 
5 into the transverse plane by an additional RF field H, . The magnetization vector then 
comprises a rotating magnetic dipole m in the transverse plane. The NMR image may 
be reconstructed by sampling the external field fi-om a series of RF dipoles rather than 
that fix)m a series of static dipoles. In this case, the geometric system fimction is given 
also given by Eq. (1), the function of the z-component of the flux B fi-om a z-oriented 

10 dipole at any point extrinsic to the image space, fi-om any single voxel. The geometric 
system function corresponding to a dipole which determines the spatial intensity 
• variations of the RF field is a band-pass for A:^ = A, as shown in APPENDIX V. In the 
liniiit, each volume element is reconstructed independently in parallel with all other 
volume elements such that the scan time is no greater than the nuclear firee induction 

15 decay (FDD) time. 

In an embodiment, the strength and duration of the rotating H , (RF) field that is 
resonant with the protons of the magnetized volume and is oriented perpendicularly to 
the direction of the magnetizing field is applied such that the final precession angle of ttie 
magnetization is 90^ (^^ = 90** ) such that the RF dipole is transverse to the static 

20 magnetizing field and perpendicular to the RF magnetic field detector 20a as shown in 
FIGURE 1 1 . The rotating NMR coordmates and flie stationary NMR coordinates are 
both shown in FIGURE 1 1 . According to Eq. (1), the signal as a function of time which 
is Fourier transformed arises fi^m each transverse RF magnetic dipole oriented in die 
^rYr -plane (rotating NMR coordinates) which is periodicaUy parallel to the y-axis and 

25 rotates about the primary field oriented parallel to the z-axis at the Larmor frequency 
(stationary NMR coordinates). For the stationary NMR case, the y-axis corresponds to 
the z-axis of Eq. (1), and /w, , the magnetic moment along the z-axis, of Eq, (1) 
corresponds to the bulk magnetization M of each voxel. 

In addition, each Foiuier component comprises an additional part that 

30 corresponds to an RF magnetic dipole ^^iaich is periodically parallel to the x-axis. This 
part also rotates about the primary field oriented parallel to the z-axis at the Larmor 
fi-equency wherein the x-axis of the stationary NMR coordinates corresponds to the x- 
axis of Eq. (1). The parts rotate at the same firequency but are orthogonal. The 
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orthogonality corresponds to a phase angle of the Fourier transform of 90^; thus, each 
Fourier con^onent contains a real and an imaginary part. Tlie magnitude of the 
ampUtude of ttie signal due to the dipole oriented parallel to the y-axis (z-axis of Eq. (1)) 
exceeds that of flie dipole oriented parallel to the x-axis (x-axis of Eq. (1)) which 

5 identifies this part of a Fourier component, hi an embodim^t, the spatial variation over 
the detector array of the part of the Fourier component with the maximum anq>litude (the 
signal due to the dipole oriented parallel to tiie y-axis) is used to determine the 
coordinate location of each voxel using the geometric system fimction of the detection of 
the z-component of a z*oriented dipole (Eq. (1)) as given in the Reconstruction 

0 Algorithm Sectioa 

In another embodiment, the spatial variation over the detector array of the lesser 
magnitude orthogonal part of each Fourier component is used to determine the 
coordinate location of each voxel using the corresponding geometric system function as 
described in the Reconstruction Algorithm Section. The system function for the case of 

5 the minor orthogonal part in terms of the coordinates given in FIGURE 9, corresponds to 
the detection of the x-component of a z-oriented dipole. The geometric system function 
is given by the function of the x-^nq>onent of the flux B* fi-om a z-oriented dipole at any 
point extrinsic to the image space, from any single voxel: 

^^'^(PTTW^ ^^^^ 

} where x, y, and z are the distances fiom the center of the voxel to the sampling point and 
is the magnetic moment along the z-axis corresponding to the bulk magnetization M 
of each voxel. 

In the reconstruction process described herein, the secondary field may be in the 
same or in a transverse orientation relative to the primary field. In each case, the 

> orientation of the secondary dipole field and the measured secondary dipole field 
component are according to Eq. (1). The relationship of the coordinate systems of the 
present invention are shown in FIGURES 9, 10, and 11. The primary and secondary 
fields are parallel and stationary in FIGURE 9 versus transverse in the rotating and 
stationary NMR coordinates shown in FIGURES 10, 11, and 12. The relationship of the 

) coordinate system of the reconstruction method of the NMR system used herein to that 
of the coordinate system of Eq, (1) shown in FIGURE 9 is according to FIGURE 11. 
The magnetization axis of Eq (1) is the z-axis; whereas, the magnetization axis of the 
stationary NMR coordinates is the z-axis shown in FIGURE 10, and the magnetization 
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axis of the NMR rotating coordinates [3,7] is tbe -axis shown in HGURB 1 1 , 
Regarding Eq. (1) and the reconstruction algorithm m the NMR case, the z-component of 
Eq, (1) is substituted with the stationary NMR y-component, the y-<x)mponCTt is 
substituted with the stationary NMR z-component, and the x-componoit is substituted 
5 vnftk the stationary NMR z-component 

The reconstruction algorithm can be a reiterative, a matrix inversion, or a Fourier 
Transform algorithm. For all reconstruction algorithms, the volume to be imaged is 
divided into volume elements called voxels, and the magnetized voxel with magnetic 
moment m, is modeled as a magnetic dipole. In an embodiment, the matrix of Fourier 

10 components fliat correspond to the NMR signal of a given voxel over the detectors is 
determined. The matrices are determined for all of the voxels. The measurements of the 
spatial variations of the transverse RF field of a given matrix is used to determine the 
coordinate location of each voxel Thus, each matrix of components associated by phase 
comprises the intensity variation ovct the sample space of the RF field of the bulk 

15 magnetization M of each voxel The matrices are the input for flie reconstruction 

algorithuL This procedure is performed either in parallel or series for each matrix. The 
bulk magnetization m^ (NMR image) is the superposition of the independent images, 
each of A^ch corresponds to a given voxel The siq>erposition of images is plotted and 
displayed. 

20 An embodiment of a matrix inversion reconstruction algorithm comprises the 

steps of 1.) using the geometric system fimction to determine the spatial intensity 
variation of the transverse RF field over the detector array, 2.) inverting the 
corresponding matrix, and 3.) multiplying the signal over the detector array by the 
inverted matrix to give the voxel sources. For example, a matrix inversion 

25 reconstruction algorithm is to determine a coefficient for each voxel mathematically (Eq. 
(41 )) or by caHbration which when multiplied by the bulk magnetization M of each 
voxel is that voxel's contribution to the signal at a given detector witii tiie corresponding 
unique phase at each detector. This is repeated for every detector, and those coeflBcients 
are used to d^ermine a matrix which, when multiplied by a column vector of the bulk 

30 magnetization M values of the voxels, gives the voltage signals at tiie detectors. This 
matrix is inverted and stored in memory. Voltages as a function of time are recorded 
over the detector array. The signal as a function of time is Fourier transforaaed to give 
the Fourier components each having an amplitude and a unique phase. In an 
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embodiineiit, of a matrix inversion reconstruction algoritbni, tbe components are 
multiplied by the inv^e matrix, to generate tiie bulk magnetization M map. 

In another embodiment of a matrix inversion reconstruction algoritiun, the matrix 
of Fourier conq)onents that correspond to the NMR signal of a given voxel over the 
S detectors is detennined. The matrices are determined for all of the voxels. The 
measurements of the spatial variations of the transverse RF field of a given matrix is 
used to determine the coordinate location of each voxel Hius, each matrix of 
conqionents associated by phase comprises the intensity variation over the sample space 
of the RF field of the bulk magnetization M of each voxel. Each matrix of components 
10 are multiplied by the inverse matrix, to generate the bulk magnetization M map. In one 
embodiment, tiie point spread of the reconstructed voxel is corrected by assigning one 
voxel above a certain threshold with the bulk magnetization M . The other voxels are 
assigned a zero value. This procedure may be repeated for all voxels. la the limit with 
sufiBcient phase resolution, each volume element is reconstructed iadepoDidently in 
IS parallel with all other volume elements such that the scan time is no greater than the 
nuclear free induction decay (FID) time. The total bulk magnetization map (NMR 
image) is the superposition of the sqyarate m^s for each magnetic moment which is 
plotted and displayed. 

In an embodiment of a reiterative reconstruction algorithm, the geometric system 
20 fimctionJs used to determine the system of linear equations which gives the intensity, 
spatial variation, and phase of the RF field over the ssmplQ space. The signal as a 
fimction of time is Fourier transformed to give the Fourier components each having an 
amplitude and a unique phase. The system of linear equations gives the voltage and 
phase from each voxel at each sensor based on the bulk magnetization M value of each 
25 voxel and the position of the voxel relative to the sensor. In an embodiment, weighting 
coeflScients are detennined based on these equations. The coefficients may be 
detennined mathematically (Eq. (41)). alternatively, they may be determined by 
calibration. The bulk magnetization M for each voxel is estimated, and the signals at 
each detector are calculated at each phase. In an embodiment, the bulk magnetization M 
30 value of each voxel times its weighting coefficient and its calculated phase at a given 
detector is compared to the measured voltage and phase. A correction is made to M of 
each voxel. This gives rise to a second, or recomputed, estimate for M of each voxel. 
The signal value from this second estimate is computed and conections are made as 
previously described. This is repeated until the correction for each reiteration 
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approaches a predefined limit which serves to indicate that the reconstruction is within 
reasonable limits of eiTor. This procedure is r^eated for all sensors. The final bulk 
magnetization map is plotted and displayed 

The general process of reconstruction by reiteration is shown according to the 

5 steps of FIGURE 13 (and is implemented in processor 20 in FIGURE 8). The image 
displayed according to the process 200 is merely a mq)ping of the bulk magnetization 
M of voxel sections of the object examined. The signal as a function of time is Fourier 
transformed to give the Fourier components each having an amplitude and a unique 
phase. The bulk magnetization M at a given location determines the intensity and 

0 spatial variation of the RF field over the sample space. Accordingly, signals produced 
by the RF magnetic sensors 1 10, in terms of volts as a function of phase, are a direct 
result of the bulk magnetization M of the voxel elements. Therefore, a reference voltage 
is generated at 210 £ix>m which the actual or measured sensor voltage at a given phase is 
subtracted at 220. The reference voltages are modeled by assxmiing a voltage with a 

5 unique phase fi-om each voxel at each sensor. Therefore, the voltage and phase fi-om 
each voxel is determined at a sensor according to a weighting based on the position of 
the voxel relative to the sensor. The weighting may be given by Eq. (1) in terms of the 
coordinates x, y, z of FIGURE 9. The phase may be determined as given in the Phase 
Angle and Associated Set of Fourier Components Section. The resulting modeled or 

D calculated voltage signals are compared at step 220, providing a difference or A{x,y, z) 
signal, weighted at step 230 to produce a weighted difference signal, which is then added 
to the previously estimated bulk magnetization M value for each voxel element at step 
240. The resulting level, available in three dimensions corresponding to the axes x, y, 
and z, is selectively displayed on the display at step 250. Having adjusted the estimated 

5 bulk magnetization M for each voxel, the calculated bulk magnetization M is 
recalculated at step 260, the resulting estimated sensor voltage at each phase is then 
compared to the actual sensor voltage at each phase at step 220, the process 200 being 
repeated until the difference is reduced to a tolerable value. This procedure is repeated 
for all sensors. The final bulk magnetization map is plotted and displayed. 

) In terms of the coordinates of Eq. (1) and FIGURE 9, the reconstraction 

algorithm using Fourier Transforms (FFT) involves exploiting the FFT to solve Eq. (41) 
given below. For the case that follows, data is acquired in the x, y, and z-directions, but 
in general, data is acquired over the dimensions which uniquely determine the bulk 
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magaetizadon map (NMR image). Also, the present analysis is for measuring the z- 
conqx>neDt of the magnetic field of a dipole onented in the z-direction; however, the 
analysis ^plies to fhe other two orthogonal components where the geometric system 
function for the z-component of the z-oriented dipole is replaced by the geometric 
S system fimction for the x or y-components of the magnetic field produced by the dipole 
where the geometric system fimction is defined below as the impulse response of the 
detector to the given component of the field of a dipole of given orientation. The sample 
space, or space over which the secondary RF field is measured, is defined in the present 
example as the three-dimensional space comprising the entire xy-plane and the positive 

10 z-axis, as shown in FIGURE 9. Other sample spaces are valid and each requires special 
consideration during the reconstruction as described below. The discrete voltages 
recorded fi-om an infinite detector array in the xy-plane which is translated to infinity 
along the z-axis starting fi:om the origin where the detector array is responsive to the z- 
component of the secondary magnetic field is given by Eq. (1), where the voltage at any 

15 point in space produced by dipoles advanced in the z-direction and advanced or delayed 
in the X and y-directions relative to the origin is given by the following Eq. (40), where 
the voltage response is Q times the secondary magnetic flux strength in the case shown 
in FIGURE 9. 

In a preferred embodiment, the data comprises the RF field in the transverse 
20 plane over the sample space shown in FIGURE 1 1 . The relationship of the coordinate 
designations of Eq. (1) and FIGURE 9 versus the NMR system and the relationship 
between the magnetic moment of a static secondary field and tfie bulk magnetization M 
of the RF field is given in the Reconstruction Algorithm Section and in the Fourier 
Transform Reconstruction Algorithm which follows. 

25 

Fourier Transform Reconstruction Algorithm 

In terms of the coordmates of Eq. (1) and FIGURE 9, the volume to be imaged is 
30 divided into volume elements called voxels and the magnetized voxel is modeled as a 
magnetic dipole m, . Data comprising the z-component of the magnetic field of a dipole 
oriented in the z-direction is acquired in the three dimensional sample space comprising 
the xy-plane and the positive z-axis. 
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The phantom can be considered to be a tfaree-dimendonal airay of magnetic 
dipoles. At any point extrinsic to the phantom, the z-conaponent of the secondary flux, 
B' , fiom any single voxel is 

S where x,y, and z are flie distances from the center ofthe voxel to the sampling point It 
is shown in APPENDICES I-FV that no geometric distribution of magnetic dipoles can 
give rise to Eq. (1). Therefore, the flux of each magnetic dipole (voxel contribution) 
forms a basis set for the flux of the array of dipoles which comprise the bulk 



10 Eq. (1) is a system fimction which gives the magnetic flux output in response to a 

magnetic dipole input at the origin. The phantom is an array of spatially advanced and 
delayed dipoles wei^ted according to the bulk magnetization of each voxel; this is the 
input function. The secondary flux is the superposition of spatially advanced and 
delayed flux, according to Eq. (1); this is the output function. Thus, the response of 

IS space to a magnetized phantom is given by the convolution of Eq. (1) with the series of 
weighted, spatially advanced and delayed dipoles fq)resenting the bulk magnetization 
mapoftfaephantoHL Thediscrete voltages recorded by a detector array ov^ the sample 
space comprising the xy-plane and the positive z-axis are 

20 (41) 
where the flux magnetizing each voxel is unity, the volume element is made unity, and 
the voltage response of each detector is Q times the secondary magnetic flux strength. 
The variables of Eq. (40) are defined as follows: 

^2.18 magnetization of the voxel located at the 

position defined by the Dirac delta function. 



30 



dipole spacing in the x, y, and z-directions, 
respectively 
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0ie dimensions of ttie image space, the total distance in flie x, y, and z-directions, 
respectively, for whidi the bulk magnetization of the voxels is nonzero 

the detector spacing in the x, y, and z-directions, respectively 

5 

Hie voltage signals recorded at the detector array ov^ the sample space is 

OO 00 00 



ai3«0 



(42) 

Eq. (41) can be represented symbolically as follows 
10 ^ = Qfex[*®/|xtt(z)] (43) 

where Q is the proportionality constant between the signal voltage and the output flux 
strength, s is (he discrete function of voltage signals recorded of the secondary flux over 
the san^le space, g is the secondary magnetic flux sampling function given as follows: 

IS In Eq. (42), h is the system fimction which is also defined as the geometric system 
function given as follows: 

^ (45) 

The system function is the inq)ulse response of the detector array where the external 
magnetizing flux is set equal to unity. In Eq. (42), / is the bulk magnetization function 
20 or NMR image function given as 

In Eq. (42), u(z) is the unitary z function which is one for positive z and zero otherwise. 
The function g discretizes the continuous voltage function V given by Eq. (40) which is 
the function h convolved with the fiinction / then multiplied by the function u(z). The 
25 discrete voltages recorded over the sample space are used in a Fourier transform 
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algoritlim d^ved herein to reconstruct the NMR image. Consider the fimction s of Eq. 
(41) which is given as follows: 



S by the fimction u(z). The function S which is the Fourier transform of is given as 



The function S is equivalent to the resultant fimction of the function ff , the Fourier 
transform of the system fimction-A , multiplied by die function F, the Fourier transform 
10 of the bulk magnetization function- /» convolved with the function UQc^X die Fourier 
transform of the unitary zfimction-u(z). The Fourier transform of the bulk 
magnetization function— / (Eq. (45)) 



j = (A®/)xw(z) (47) 
Eq. (46) is equivaloit to the function h convolved with the function / then multiplied 



follows: 



5 = (//xF)®(7(jtJ 



(48) 



/ = 




(49) 



IS 



15 




(50) 



where x^=njc^; y„ = n^k^ ; z„ = -n^k^ . The Fourier transform of 
u(z) = lforz>0 and u(z) = 1 /ar z < 0 is [8] 



(51) 



The Fourier transfonn of flie system fimction (See APPENDIX V) 



20 




(52) 



is 




(53) 




(54) 



25 where 
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X 
y 

z 



(55) 



The resultant fimctioii of the product of the functioiK H (Eq. (52)) and F (Bq. (49)) th^ 
convolved with the function (Eq. (SO)) is 

5 The result is given as follows, and the derivation appears in APPENDIX VI. 

1 +/a/2 +/,/2 , . r _9 ttP ' 

..£;/2 :.ii2* l2[ii^+y*py j 

(57) 

Substitution of the Fourier transform of the system function, H (Eq. (S2)X the Fourier 



transform of the bulk magnetization function, F (Eq. (49)), and factoring out in the 



1 0 second term gives 

, e 



5 = 



r 

+f»/2 +/|/2 



1 



2il+r 



(58) 

Multiplication by the complex conjugates gives 
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r 




L K. 





1 









) 

(59) 



Factoring out the common tenns gives 
+*i r.— 



(60) 



5 Combiniiig exponential tenns gives 
S- 



<0 



(61) 



The function S divided by tiie function H is 



(62) 



10 The inverse Fourier transfonn of the function S divided by the function H is given as 
follows, where the symbol is defined as the inverse Fourier transform of the 
function Q, 



.r^i_„ '^..-o .... ^^^^ ^^^^ 

l^e-*'''-' + ^^''''■l - I 



The inverse Fourier transfonn of the function F is 

15 ^\F)^f=Y. E E ;t.,,,,5[x-w-j.„,z-z.] 



(64) 
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Specifically, (he inverse Fourier transfonn of the function F in tenns of the dipole 
spacing is 



t-'j H/2 +4/2 , , V 

Z Z E -'•^'"^ ' 



(65) 



Consider the general Fourier transform pair in cylindrical coordinates: 
5 «(r,e,z)oG(p,^,a))(66) 

Under circular symmetry, that is when g is independent of 9 (and hence G independent 
of ^ ) the inverse Fourier transform is [9] 



0 — 



(67) 



The inverse Fourier transform of Eq. (62): 



10 F'M 
I 

I 



h z;.<0 



(68) 



J 



given by Eq. (64) and Eq. (66) is 



O«,=0 y2* iH— /,/2)t 



-so 0 



+27r J Je*'^''-'y.[V)]*//*''*, 

-cXJ 0 



15 



(69) 

The solution of Eq. (68) ^peais as follows, and the derivation appears in APPENDIX 

vn. 
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[»,,«0 y^^htl x,^U2 I ^p^ -!|J 



Mverse Transfonn 1 



S[z] 



(70) 



S Thus, the invoke Fouiio- transfonn of — (Eq. (62)) is givoi by Eq. (69). 

H 



= F- 



/*».=0 y,— /a/2 jt.o-(y2 L ^P^ 



10 The convolution replaces each coordinate with a spatially shifted coordinate. 

[SI K Z 



(71) 



(72) 



The result of the evaluation of Eq, (71) at each coordinate x^^y^^z^ is 
Solving for the bulk magnetization of each dipole gives 



(73) 



15 



2^ 



(74) 



The solution of the bulk magnetization of each dipole follows from Eq. (73). 
Discrete values of the voltages produced at the detector array due to the secondary 
magnetic field are recorded during a scan which represent discrete values of function s 
(Eqs. (41-42)); thus, in the reconstruction algorithm that follows, discrete FoutIct and 
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Inverse Fourier transforms replace tiie corresponding continuous fiinctions of Eq. (73) of 
the previous analysis. 

Discrete values of H of £qs. (52-54), die Fourier transform of the system 
function, rq)lace die values of the continuous function. Furtiiennore^ eadi sanq)le 
5 voltage of an actual scan is not truly a point sanqile, but is equivalent to a sample and 
hold which is obtained by inverting the grid antennae matrices or which is read directiy 
fiom a micro-antennae as described in the Finite Detector Length Sectioit The spectrum 
of a function discretely recorded as values, each of which is equivalent to a sample and 
hold, can be converted to the spectrum of the function discretely recorded as point 
10 samples by dividing the former spectrum by an appropriate sine function. This operation 
is performed and is described in detail in the Finite Detector Length Section From these 
calculated point samples, the bulk magnetization fimction (NMR image) is obtained 
following the operations of Eq. (73) as given below. 

15 Reconstruction Algorithm 

1) Record the RF NMR signal at discrete points ovot the sample space. Each point is 
designated (x, y, z, RF) and each RF value is an element in matrix A. 

2) 20 Invert the detector grid matrices defined by the noncommon areas of the overlapping 

elements of the detector array described in the Finite Detector Length Section to obtain 
the sample and hold voltages which form Matrix A' (if micro-antennas are used, form 
Matrix A' of the recorded voltages directly). 

3) 25 Fourier transform the time dependent signals to give the intensity and phase of each 

component. The NMR signal of each voxel at any given detector gives rise to a Fourier 
component with a unique phase angle relative to the Fourier component of any other 
voxel of the phantom at that detector. The matrix of Fourier components that correspond 
to the NMR signal of a given voxel over the detectors is determined. This may be 
30 achieved by using a first component having a phase angle and calculating the phase angle 
as a fimction of spatial position of the first detector relative to any other detector and 
identifying the component at each detector having the calculated phase angle. The 
matrices are determined for all of the voxels. The measurements of the spatial 
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variations of the transverse RF field of a given matrix is used to determine tiie coordinate 
location of each voxel. Thus, eadi matrix of components associated by phase comprises 
the intensity variation over the sample space of the RF field of the bulk magnetization M 
of each voxel. 

5 

4) Three-dim^onally Fourier transform each matrix , using a discrete Fourier 
transform formula such as that which appears in McC Siebert [10] to obtain each matrix 

of elem^ts at frequencies corresponding to the spatial recordmg interval in each 
direction. 

10 ^„(W)=>J».(M,A) 

5) Multiply each element of each matrix by a value which is the inverse of the Fourier 
transform of a square wave evaluated at the same fipequency as ttie element where the 
square wave corresponds to a sample and hold operation performed on the continuous 

1 S voltage fimction produced at the detector array by die secondary field. This 

multiplication forms each matrix B^ * which is the discrete spectrum of the continuous 
voltage fimction discretely point sanq>led (See the Finite Detector Length Section for 
details of this operation). 

' smc(fc„A^,*:,) ^ 

20 

6) Multiply each element of each matrix B„ * by the value which is the inverse (reciprocal) 
of the Fourier transform of the system fimction evaluated at the same fi^equency as the 
element to form each matrix B^ 

25 

7) Inverse three-dimensionally Fourier transform each matrix B^** using the discrete 
inverse Fourier transform formula such as that which appears in McC. Siebert [10] to 
form each matrix whose elements correspond to the bulk magnetization of the dipoles 
at the points of the image space of spatial interval s^propriate for the frequency spacing 

30 of points of matrix 5„ * * . 

^«**(*,.^.*,)=>C,(x,;/,z) 
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TV 

8) Divide each element of each matrix C, by — r to coirect for the restriction that the 

sample space is defined as z greater than zero (z > 0). This operation creates each matrix 
which is tiie bulk magnetization M mip ^]MR image). 

2^1 



In the NMR case shown in FIGURES 8, 10, 11, and 12, Step 7 of the 
Reconstruction Algorithm of the Fourier Transform Reconstruction Algorithm Section is 
10 equivalent to the graeral case except for the coordinate designations and die relationship 
between die magnetic moment of a static secondary field and the bulk magnetization M 
of die RF field wherein , the magnetic moment along the z*axis, of Eq. (1) 
corresponds to the bulk magnetization M of each voxel: 



Z'^S Divide each element of each matrix C„ by to correct for the restriction that the 

sample space is defined as y greater than zero, y > 0. This operation creates each matrix 
which is tiie bulk magnetization M raap. 

In one embodiment, the point spread of the reconstructed voxel is corrected by assigning 
20 one voxel above a certain threshold with the bulk magnetization M . The other voxels 
are assigned a zero value. This procedure may be repeated for all voxels. In the limit 
with sufficient phase resolution, each volume element is reconstructed indq)endently in 
parallel with all other volume elements such that the scan time is no greater than the 
nuclear fi^ee induction decay (FID) time. 

25 

9) Superimpose the separate bulk magnetization M maps corresponding to each matrix A„ 
Plot the bulk magnetization M with the same spatial interval in each direction as the 
sampling interval in each direction over all matrices (i.e..tfae total bulk magaetization 
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M map which is plotted and displayed is die superposition of tiie separate msps of each 
magnetic moment corresponding a matrix A^). 

(The above stq)s for each matrix relate generally to the program implementation shown 
5 in the listing of the Exemplary Reconstruction Program as follows. The above Steps 1 
and 2 relate to the Data Statements beginning at lines SO; and Step 4 relates to the Z 
and Y FFT operations of lines 254, 455 and 972, respectively. Steps 5 and 6 are 
implemented by the processes of lines 2050, 2155 and 2340; and Step 7 relates to the X, 
Y and Z inverse transform of lines 3170, 3422 and 3580, respectively. Step 8 and 8* . 
10 relates to the correction and normalization process of line 4677.) 

Exemplary Reconstruction Program 

10! 4D-MRI ALGORITHM 

25 OPEN "4D-MRIPROTOTyPBI.LIS" FOR OUTPUT AS #1 
15 26 C».05 

27 DIN X(9,9,9) 

28 FOR W=l TO 9 STEP 1 

29 FOR T=l TO 9 STEP 1 

30 FOR S=l TO 9 STEP 1 
20 31 X{W,T,S)»0 

32 NEXT S 

33 NEXT T 

34 NEXT W 

35 X(5,5,6)=l 
25 36 DIM DI(9,9) 

37 PRINT #1, "DIPOLB PHANTOM" 

40 FOR Q»l TO 9 STEP 1 

41 FOR R«l TO 9 STEP 1 

42 FOR U=l TO 9 STEP 1 
30 43 LET H=X(U,R,Q) 

44 LET DI(U,R)=H 

45 NEXT U 

46 NEXT R 

47 MAT PRINT #1,DX, 
35 48 NEXT Q 

50 DATA 0,0,0,0,0,0,0,0,0,0,0,0,0 

52 DATA 0,0, 0,0,0, -.6, -.9, -1.1, -1.2. -1.8, 5. 7. 3. 5,. 4, -.4, -.5, -.7, -.6,. 7 
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54 DATA 4.2, 19.7,7.8,-1, -.9, -.3, -.4, -.7, .7,13.9,25.65,10.1, .9, -.8, -.3, 
.6 

56 DATA 4, 2. 7, 10. 7, 12. 4, 5. 6, 2. 6, 5, -. 3, -.5, -.6, -.6, .3,3.7, .9,-. 5,- 
1.0 
5 58 DATA 

0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0 
60 DATA 0,0,-.5,-.8,-.8,-1.0,-1.2,-.8,6.9, .2,-.2,-.4,-.6,- 
.6,1.2,4.1,12.7,0, .6 

62 DATA 6, -.3, -.5, -.7, -.3, 3. 7, 9. 05, 7. 6, 1.3, -.5,-. 3,-. 4,- 
10 .2,2.7,4.3,10.6, .8 

64 DATA 1.7,-.4,-.3,-,5,-.4,-.l,l.l,l.l,1.7, .2,- 
.6,0,0,0,0,0,0,0,0,0,0,0,0,0 

66 DATA 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0,0, -.5, -.6, -.6, -.5, 
.5 

15 68 DATA .3,2.5, .l,-.l, -.3, -.4, -.3, 1.3, 3. 9, 8. 1,1. 7, .9, -.3,-. 3,-. 4,-. 6, - 
.8,1.5 

70 DATA 1.65, 4. 7, 1.3, -.3, -.2, -.2, .1,2.1,3.5,7.1, .7, .7, -.3, -.2, -.4,- 
.2, .2,1.2 

72 DATA 1.1,1.4, .4,- 
20 .3,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0 

74 DATA 0,0,0, 0,0, 0,0, 0,0,0, -.4, -.5,0,-. 1, -.1,1.7,1.5, .2,-. 1,-. 2, - 
.2,0,1.3,3.1 

76 DATA 4. 1,3. 4, .9,0, -.2, -.3, -.4, -.4, 1.1, -65, 2. 9, 1.1,-. 2,-. 1,- 
.1, .2,1.2,1.6 

25 78 DATA 4.9, .1, .5, -.3, -.2, -.3, -.1,0, .3,1.3, .6, .1, - 
.1,0,0,0,0,0,0,0,0,0,0,0,0 

80 DATA 0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0,0, -.3,- 
.2, .1, .9,1.6,1.6 

82 DATA . 4,0,0, -.l,-.l,. 2,1.0,2. 2, 2. 7, 3. 7, .5,0, -.2,-. 2,-. 2,- 
30 .1, .6, .25,1.6, .7 

84 DATA -.2,-.l,-.l,.l,.7,.9,2.8,0,.3,-.2,-.l,-.2,-.l,0,.2,.3,.3,0,- 

.1,0,0,0 

85 DATA 0,0,0,0,0,0,0,0,0,0,0,0,0,0,0 

86 DATA 0, 0^0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0,0, -.2,- 
35 .1, ,1,1.0,1.4,1.3, .2,0,0 

88 DATA ^.1,0, .2, .7,1.5,2.0,2.9, .3,0,-. l,-.l, .1,0, .7, .25, . 9,-4, - 
.1,0,0, .1, .5 

90 DATA .5,1.3,0, .2, -.1,0, -.1,0,-1, ,2,. 1, .2,0, - 
.1,0,0,0,0,0,0,0,0,0,0,0,0,0,0 
40 92 DATA 0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,- 
.1,0,.1,.8,.9,1.0,.1,0,0 



wo 02/16956 PCT/USOl/25954 

94 DATA 0,0, .2, .5,1.1,1.3,2.0, .3, 0,0, -.1,0,0, .3, .15, .2, - 
.1,0,0, .1, .4, .3, .9 

96 DATA 0, .2,-.l,0,- 

.1, .1, .1, .1, .1, .2,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0 
5 98 DATA 0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,- 
.1,0, 0, .7, .7, .6, .1,0, .1 
100 DATA 0, .1, .2, .3, .6, .8,1.0, .2,0,0, - 
.1,0,0, .1, .15, .3,.1,0,0,0,0,.2,.2,.5 
102 DATA 0, .1, - 

10 .1,0,0,0, .1, .1, .1, .1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0 
104 DATA 

0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,.5,.4,.4,.l,0,0,0r0,.l,.l 
106 DATA .4, .5, .6, .1,0,0,0,0,0,0, .05, .1,0,0,0,0,0, .1, .1, .3,0, , 1,- 
.1,0,0,0,0,0 

15 108 DATA .1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0 

120 DIM V(9,9,9) 

130 FOR Z=l TO 9 STEP 1 

140 FOR Q=l TO 9 STEP 1 

145 FOR 0=1 TO 9 STEP 1 
20 147 READ V(Z,Q,U) 

222 NEXT U 

223 NEXT Q 

224 NEXT Z 

225 PRINT #1, "VOLTAGE DATA" 
25 226 DIM V0(9,9) 

227 FOR Q=l TO 9 STEP 1 

228 FOR R=l TO 9 STEP 1 

229 FOR U=l TO 9 STEP 1 

230 LET H»V(U,R,Q) 
30 231 LET VO(U,R)=H 

232 NEXT U 

233 NEXT R 

234 MAT PRINT #1, VP , 

235 NEXT Q 

35 254 IFFT THE ROWS OF SAMPLED VOLTAGES IN THE X DIRECTION 

255 DIM MR (9) 

256 DIM MI (9) 

257 DIM R{9) 

258 DIM RV(9,9,9) 
40 259 DIM IV(9,9,9) 

260 FOR V=0 TO 9 STEP 1 
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270 FOR M=l TO 9 STBP 1 
280 FOR N=l TO 9 STBP 1 
290 LET Y=V(M,N,V) 
300 IiBT R(»)=y 
5 305 NEXT N 

310 EXTERNAL SUB PFT(DIM{) ,DIM{) ,OTM() ) 
320 CAIiL FFT(R(),MR(),MI{)) 
330 FOR W=l TO 9 STBP 1 
340 LET Y=MR(N) 
10 350 LET RV(M,M,V)«Y 
360 KEXT N 

370 FOR N»l TO 9 STBP 1 

380 LET G»HI(N) 
•390 LET IV(M,N,V)=G 
15 400 NEXT N 

410 NEXT M 

415 NEXT V 

417 DIM RV0(9,9) 

41*9 DIM IV0{9,9) 
20 420 PRINT #1,»»RV" 

421 FOR Q=l TO 9 STEP 1 

422 FOR R=l TO 9 STEP 1 

423 FOR U=l TO 9 STEP 1 

424 LET H=RV(U,R,Q) 
25. 425 LET RVO(U,R)=H 

426 NEXT U 

427 NEXT R 

428 MAT PRINT #1, RVO, 

429 NEXT Q 

30 430 PRINT #1, "IV" 

431 FOR Q=sl TO 9 STEP 1 

432 FOR R=l TO 9 STEP 1 

433 FOR U=l TO 9 STBP 1 

434 LET H=IV(U,R,Q) 
35 435 LET IVO{U,R)=H 

436 NEXT U 

437 NEXT R 

438 MAT PRINT #1, IVO, 

439 NEXT Q 

40 455 J PPT THE COLUMNS OF THE SAMPLED VOLTAGE IN THE Z DIRECTION 
460 DIM RRV(9,9,9) 
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470 DIM IRV(9,9,9) 
475 FOR V=l TO 9 STEP 1 
480 FOR N=l TO 9 STEP 1 
490 FOR M=:l TO 9 STEP 1 
5 500 LET Y=RV(M,N,V) 
510 LET R(M)=y 
520 NEXT M 

530 CALL FFT(R(),MR(),NI()) 
540 FOR M=l TO 9 STEP 1 
10 550 LET H=MR(M) 

560 LET RRV(M,N,V)=H 
570 NEXT M 

580 FOR M-1 TO 9 STEP 1 
590 LET 6»MI(M) 
15 600 LET IRV{M,N,V)=G 
610 NEXT M 
620 NEXT N 

625 NEXT V 

626 DIM RRVO(9,9} 
20 627 DIM IRV0(9,9) 

630 PRINT #1,"RRV« 

631 FOR Q=l TO 9 STEP 1 

632 FOR R=l TO 9 STEP 1 

633 FOR 0 =1 TO 9 STEP 1 
25 634 LET H=RRV(U,R,Q) 

637 LET RRVO(U,R)«H 

638 NEXT U 

639 NEXT R 

640 MAT PRINT #1, RRVO, 
30 641 NEXT Q 

650 PRINT #1, "IRV" 

651 FOR Q»l TO 9 STEP 1 

653 FOR R=l TO 9 STEP 1 

654 FOR U=l TO 9 STEP 1 
35 655 LET H= IRV{U,R,Q) 

656 LET IRVO(U,R)=H 

657 NEXT U 

658 NEXT R 

659 MAT PRINT #1, IRVO, 
40 660 NEXT Q 

661 DIM RIV(9,9,9) 
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662 DIM IIV(9,9,9) 
685 FOR V=l TO 9 STEP 1 
690 FOR N=l TO 9 STEP 1 
700 FOR M«l TO 9 STEP 1 
5 710 LET y=IV(M,N,V) 
760 IiET R(M)»Y 
770 NEXT M 

830 CALI, FFT(R(),MR(),MI()) 
840 FOR M=l TO 9 STEP 1 
10 850 LET H=MR(M) 

860 LET RIV(M,N,V)=H 
870 NEXT M 

872 FOR M= 1 TO 9 STEP 1 

873 LET G»MI(H) 

15 874 LET IIV(M,N,V)=G 

875 NEXT M 

876 NEXT N 

877 NEXT V 

878 DIM RIV0{9,9) 
20 879 DIM IIVO(9,9) 

880 PRINT #1, "RIV" 

881 FOR Q=l TO 9 STEP 1 

882 FOR Rsl TO 9 STEP 1 

883 FOR U«l TO 9 STEP 1 
25 884 LET H=RIV(U,R,Q) 

885 LET RIVO{U,R)«H 

886 NEXT U 

887 NEXT R 

888 MAT PRINT #1, RIVO, 
30 889 NEXT Q 

890 PRINT #1, "IIV" 

891 FOR Q=l TO 9 STEP 1 

892 FOR R=l TO 9 STEP 1 

893 FOR U=l TO 9 STEP 1 
35 894 LET H=IIV{U,R,Q) 

895 LET IIVO(U,R)«H 

896 NEXT U 

897 NEXT R 

898 MAT PRINT #1, IIVO, 
40 899 NEXT Q 

900 DIM RVA(9,9,9) 
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901 DIM IVA(9,9,9) 
904 DIM RVAO(9,9) 
906 DIM IVAO(9,9) 
908 FOR Q=l TO 9 STEP 1 
5 910 FOR R«l TO 9 STEP 1 

911 FOR Usl TO 9 STEP 1 

912 IiBT H«IIV(a,R,Q) 

913 H»(-1)*H 

914 LET G=RRV(U,R,Q) 
10 915 LET L==G+H 

916 LET RVA(U,R,Q)oL 

917 NEXT V 

918 NEXT R 

919 NEXT Q 

15 920 FOR Q«l TO 9 STEP 1 

921 FOR R=l TO 9 STEP 1 

922 FOR U=l TO 9 STEP 1 

923 LET H=IRV(U,R,Q) 

924 LET L=RIV(U,R,Q) 
20 925 LET IVA(U,R,Q)»H4L 

927 NEXT U 

928 NEXT R 
930 NEXT Q 

932 PRINT #1, "RVA" 
25 934 FOR Q=l TO 9 STEP 1 
936 FOR R«l TO 9 STEP 1 
938 FOR U=l TO 9 STEP 1 
940 LET H=RVA(U,R,Q) 
942 LET RVAO(U,R)=H 
30 944 NEXT U 
946 NEXT R 

948 MAT PRINT #1, RVAO, 

950 NEXT Q 

952 PRINT #1, "IVA" 

35 954 FOR Q=l TO 9 STEP 1 
956 FOR R»l TO 9 STEP 1 
958 FOR U=l TO 9 STEP 1 
960 LET H=IVA(U,R,Q) 
962 LET IVAO(U,R)=H 

40 964 NEXT U 
966 NEXT R 
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968 MAT PRUTT #1, IVAO, 
970 KEXT Q 

972IFFT THE ROWS OP THE SAMPLED VOLTAGES IN THE Y DIRECTION 
980 DIM RVAY(9,9,9) 
5 990 DIM IRVAY(9,9,9) 
992 DIM RVAYO(9,9) 
994 DIM IRVAY0(9,9) 

1012 FOR M=l TO 9 STEP 1 

1013 FOR N«l TO 9 STEP 1 
10 1014 FOR Vol TO 9 STEP 1 

1015 LET Y=RVA{M,N,V) 

1016 LET R(V)=Y 

1017 NEXT V 

1018 CALL FFT(R() ,MR() ,MI()) 
15 1019 FOR V=l TO 9 STEP 1 

1020 LET H=MR(V) 

1021 LET RVAY {M,N,V)=H 

1022 NEXT V 

1023 FOR V=l TO 9 STEP 1 
20 1024 LET G»MI(V) 

1025 LET IRVAY(M,N,V)»G 

1026 NEXT V 

1027 NEXT N 

1028 NEXT M 

25 1030 PRINT #1,"RVAY" 

1032 FOR Q=l TO 9 STEP 1 

1034 FOR R=l TO 9 STEP 1 

1036 FOR U=l TO 9 STEP 1 

1038 LET H=RVAY(U,R,Q) 
30 1040 LET RVAYO(U,R)=H 

1042 NEXT U 

1044 NEXT R 

1046 MAT PRINT #1,RVAY0, 

1048 NEXT Q 
35 1050 PRINT #1, "IRVAY" 

1052 FOR Q=l TO 9 STEP 1 

1054 FOR R=l TO 9 STEP 1 

1056 FOR U=l TO 9 STEP 1 

1058 LET H=IRVAY(U,R,Q) 
40 1060 LET IRVAyO(U,R)=H 

1062 NEXT U 
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1064 NEXT R 

1066 MAT PRINT #1,IRVAY0, 
1068 NEXT Q 
1070 DIM RIVY(9,9,9) 
5 1080 DIM IIVY(9,9,9) 

1085 DIM RIVY0(9,9) 

1086 DIM IIVYO(9,9) 
1090 FOR M»l TO 9 STEP 1 
1100 FOR N«l TO 9 STEP 1 

10 1138 FOR V=:l TO 9 STEP 1 

1139 LET y=IVA(M,N,V) 

1140 LET R(V)«y 

1141 NEXT V 

1142 CALL FFT(R() ,MR() ,MI()) 
15 1143 FOR V=l TO 9 STEP 1 

1144 LET H=MR(V) 

1145 LET RIvy(M,N,V)aH 

1146 NEXT V 

1147 FOR V=l TO 9 STEP 1 
20 1148 LET G=MI(V) 

1149 LET IIVY(M,N,V)=G 

1150 NEXT V 

1151 NEXT N 

1152 NEXT M 

25 1153 PRINT #1. "RIVY" 

1160 FOR Q=l TO 9 STEP 1 

1162 FOR R=l TO 9 STEP 1 

1164 FOR U=l TO 9 STEP 1 

1166 LET H=RIVy(U,R,Q) 
30 1170 LET RIVyO(U,R)=H 

1172 NEXT U 

1174 NEXT R 

1178 MAT PRINT #1, RIVYO, 

1180 NEXT Q 
35 1185 PRINT #1, "IIVY" 

1190 FOR Q=l TO 9 STEP 1 

1200 FOR R=l TO 9 STEP 1 

1210 FOR Usl TO 9 STEP 1 

1212 LET H= IIVY(U,R,Q) 
40 1214 LET IIVYO(U,R)=H 

1216 NEXT U 
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1218 NEXT R 

1220 MAT PRINT #1,IIVY0, 

1222 NEXT Q 

1230 DIM YRVA(9,9,9) 
5 1240 DIM yiVA{9,9,9) 

1241 FOR Q=l TO 9 STEP 1 

1250 FOR R«l TO 9 STEP 1 

1260 FOR U=l TO 9 STEP 1 

1270 LET LoIIVY(0,R,Q)^ 
10 1280 LET B=RVAY(U,R,Q) 

1290 LET YRVA(U^R,Q)«B-L 

1300 NEXT U 

1310 NEXT R 

1320 NEXT Q 
15 1330 FOR Q=l TO 9 STEP 1 

1340 FOR R=l TO 9 STEP 1 

1345 FOR U=l TO 9 STEP 1 

1350 LET.H=RIVY(U,R,Q) 

1360 LET L=IRVAY{U,R,Q) 
20 1370 LET yiVA(U,R,Q)«L+H 

1380 NEXT U 

1390 NEXT R 

1400 NEXT Q 

1410 PRINT #1, "YRVA" 
25 1412 DIM YRVAO(9,9) 

1414 DIM YIVAO(9,9) 

1420 FOR Q=l TO 9 STEP 1 

143 0 FOR R=l TO 9 STEP 1 

1440 FOR U=l TO 9 STEP 1 
30 1450 LET H=YRVA(U,R,Q) 

1460 LET YRVAO(U,R)=H 

1470 NEXT U 

1480 NEXT R 

1490 MAT PRINT #1,YRVA0, 

35 1500 NEXT Q 

1510 PRINT #1, "YIVA" 
1520 FOR Q«l TO 9 STEP 1 
1530 FOR R=l TO 9 STEP 1 
1540 FOR U=l TO 9 STEP 1 

40 1545 LET H=YIVA(U,R,Q) 
1550 LET YIVAO(U,R)=H 
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1560 NEXT U 
1570 NEXT R 

1580 MAT PRINT #1,YIVA0, 
1590 NEXT Q 

5 2050 IGENERATE THE DISCRETE SPECTRUM OF THE SYSTEM FUNCTION AND THE 
SINC 

2055 IFONCTION OF THE SAMPLE AND HOLD CORRESPONDING TO THE FINITE 
DETECTOR 
2057 (DIMENSIONS 
10 2060 DIM SFH(9,9,9) 

2061 DIM SFHO(9,9) 

2062 DIM SINC (9,9,9) 

2063 DIM SINC0(9,9) 
2065 LET SFH(4,4,4)o4*PI 

15 2066 FOR M=-4 TO 4 STEP 1 
2070 FOR Vb-4 TO 4 STEP 1 

2080 FOR N=-4 TO 4 STEP 1 

2081 H=.04 

2082 J^ABS (M) +ABS (N) 4>ABS (V) 
20 2085 IF JbO then GO TO 2098 

2090 T=4*PI* ( (2*PI*N/9*1/ . 1) *2+ ( (2*PI*V/9*l/ . 1) ^2) ) 
2092 B= (2*PI*N/9*1/ . 1) "^2+ (2*PI*M/9*l/C) ^2+ (2*PI*V/9*l/ . 1) ^2 
2095 LET SFH(M+5,N+5,V+5)»T/B 
2098 G=ABS (N) +ABS (V) 
25 2100 IF 0=0 THEN GO TO 2114 
2101 A«l 

2105 IF N=0 THEN GO TO 2107 

2106 A=SIN(2*PI*10*N/9*.1)/{PI*N/9*10) 

2107 B=l 

30 2108 IF V=0 THEN GO TO 2110 

2109 B=SIN {2*PI*10*V/9* • 1) / (PI*V/9*10) 

2110 H»A*B 

2111 IF N=0 THEN H=.2*H 

2112 IF V=0 THEN H=.2*H 

35 2114 LET SINC(M+5,N+5,V+5)=H 

2130 NEXT N 

2131 NEXT V 

2132 NEXT M 

2135 PRINT #1, "SFH" 
40 2136 FOR 0=1 TO 9 STEP 1 
2137 FOR R=l TO 9 STEP 1 



wo 02/16956 PCT/US01/259S4 

70 

2138 FOR U=l TO 9 STBP 1 

2139 LET S»SFH(U,R,Q) 

2140 LET SFHO(U,R)«S 

2141 NEXT U 
S 2142 HEXT R 

2143 MAT mim #1, SFHO, 

2144 NEXT Q 

2145 PRINT #1, "SINC" 

2146 FOR Q=l TO 9 STBP 1 
10 2147 FOR R«l TO 9 STEP 1 

2148 FOR U=l TO 9 STBP 1 

2149 LET S»SINC(U,R,Q) 

2150 LET SINCO(U,R)sS 

2151 NEXT U 
15 2152 NEXT R 

2153 MAT PRINT #1, SINCX), 

2154 NEXT Q 

2155 ! INVERSE THE DISCRETE SPECTRUM OF THE SYSTEM FUNCTION AND THE SINC 
2157 ! FUNCTION 

20 2160 DIM HR(9,9,9) 

2170 DIM HRO(9,9) 

2171 DIM SINCR(9,9,9) 
2175 FOR V=l TO 9 STBP 1 
2180 FOR M=l TO 9 STEP 1 

25 2190 FOR N=l TO 9 STEP 1 

2200 LET yoSFH(M,N,V) 

2201 LET H=SINC(M,N,V) 

2210 IF Y=0 THEN GO TO 2212 

2211 Y=l/Y 

30 2212 IF HaO THEN GO TO 2230 

2221 H»l/H 

2230 LET HR{M,N,V)n:Y 

2235 LET SINCR(M,N, V) =H 

2240 NEXT N 
35 2250 NEXT M 

2260 NEXT V 

2310 PRINT #1,"HR" 

2311 FOR Q=l TO 9 STEP 1 

2312 FOR R=l TO 9 STBP 1 
40 2313 FOR U«l TO 9 STEP 1 

2314 LET S=HR(U,R,Q) 
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2315 LET HRO(U,R)sS 

2316 NEXT U 

2317 NEXT R 

2320 MAT PRINT #1,HR0, 
5 2321 NEXT Q 

2340 I DIVIDE THE TRANSFORMED DATA BY THE TRANSFORM OF THE SYSTEM 
FDNCTION 

2345 lAND THE SINC FUNCTION 
3030 DIM FR(9,9,9) 
10 3050 DIM FI(9,9,9) 
3052 DIM FRO (9, 9) 
3054 DIM FIO(9,9) 

3056 FOR V»l TO 9 STEP 1 

3057 FOR M«l TO 9 STEP 1 
15 3058 FOR N=l TO 9 STEP 1 

3059 T=YRVA(M,N,V) 

3060 S=HR(M,N,V) 
3062 Ii^SINCR(M,N,V) 
3065 K;=S*T*L 

20 3066 LET FR{M,N,V}=K 

3067 NEXT N 

3068 NEXT M 
3070 NEXT V 

3080 FOR V=l TO 9 STEP 1 
25 3081 FOR M*l TO 9 STEP 1 

3082 FOR N»l TO 9 STEP 1 

3083 H=SINCR(M,N,V) 

3093 K= HR(M,N,V) 

3094 L« YIVA(M,N,V) 
30 3095 E«K*L*H 

3096 LET FI(M,N,V)=E 

3097 NEXT N 

3098 NEXT M 
3100 NEXT V 

35 3130 PRINT #1, "FR" 

3131 FOR Q»l TO 9 STEP 1 

3132 FOR R»l TO 9 STEP 1 

3133 FOR U«l TO 9 STEP 1 

3134 LET H== FR(0,R,Q) 
40 3135 LET FRO(U,R)=H 

3136 NEXT U 
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3137 NEXT R 

3138 HAT PRINT #1, FRO, 

3140 NEXT Q 

3141 PRINT #1, "FI" 

5 3142 FOR Qal TO 9 STEP 1 

3143 FOR Ral TO 9 STEP 1 

3144 FOR U»l TO 9 STEP 1 

3145 LET H»FI(U«R,Q) 

3146 LET FIO(U,R)»R 
10 3156 NEXT XJ 

3157 NEXT R 

3158 MAT PRINT #1, FIO, 
3160 NEXT Q 

3170 {INVERSE TRANSFORM THE ROWS IN THE X DIRECTION 
15 3180 DIM RP{9,9,9) 

3185 DIM RPO(9,9) 

3187 DIM MIFO{9,9) 

3190. DIM MIF(9,9,9) 

3195 DIM I (9) 
20 3196 FOR Vol TO 9 STEP 1 

3200 FOR M»l TO 9 STEP 1 

3210 FOR N=l TO 9 STEP 1 

3220 LET Y=PR(M,N,V) 

3230 LET R(N)=Y 
25 3240 NEXT N 

3250 FOR N=l TO 9 STEP 1 

3260 LET y=FI(M,N,V) 

3270 LET I(N)=Y 

3280 NEXT N 

30 3285 EXTERNAL SOT IPT (DIM () ,DIM () ,DIM () , DIM {) ) 

3290 CALL IFT(R() , I () ,MR() ,MI () ) 

3300 FOR N»l TO 9 STEP 1 

3310 LET Y>=MR(N) 

3320 LET RF(M,N,V)=Y 
35 3330 NEXT N 

3340 FOR N=l TO 9 STEP 1 

3350 LET Y»MI(N) 

3360 LET MIF<M,N,V)=Y 

3370 NEXT N 
40 3380 NEXT M 

3385 NEXT V 
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3390 PRINT #1, "RF« 

3391 FOR Q«l TO 9 STEP 1 

3392 FOR R=l TO 9 STEP 1 

3393 FOR U=l TO 9 STEP 1 
5 3394 JjET H«RF(U,R,Q) 

3395 I«BT RFO{U,R)eH 

3396 NEXT U 

3397 NEXT R 

3400 MAT PRINT #1, RFC, ^ 
10 3405 NEXT Q 

3410 PRINT #1," MIF " 

3411 FOR Q=l TO 9 STEP 1 

3412 FOR R»l TO 9 STEP 1 

3413 FOR U=l TO 9 STEP 1 
15 3414 LET H=s MIF(U,R,Q) 

3415 LET MIFO(U,R)=H 

3416 NEXT U 

3417 NEXT R 

3420 MAT PRINT #1, MIFO, 
20 3421 NEXT Q 

3422IINVERSE TRANSFORM THE ROWS IN THE Y DIRECTION 
3430 DIM RFY(9,9,9) 

3432 DIM MIFY(9,9,9) 

3433 DIM RFYO(9,9) 
25 3434 DIM MIPYO(9,9) 

3435 FOR M=l TO 9 STEP 1 
3440 FOR N=l TO 9 STEP 1 
3450 FOR V=l TO 9 STEP 1 
3460 LET y«RF(M,N,V) 
30 3470 LET R{V)=Y 
3480 NEXT V 

3490 FOR V=l TO 9 STEP 1 
3500 LET Y=MIF(M,N,V) 
3510 LET I(V)=Y 
35 3520 NEXT V 

3525 EXTERNAL SUB IFT(DIM() ,DIM() ,DIM () ,DIM () ) 

3526 CALL IFT{R{) , I <) ,MR () ,MI {) ) 

3527 FOR V=l TO 9 STEP 1 

3528 LET y=MR(V) 

40 3529 LET RFY(M,N,V)=Y 
3530 NEXT V 
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3531 FOR V=l TO 9 STEP 1 

3532 LET Y=MI (V) 

3533 LET MIFY{M,M,y) =Y 

3534 NEXT V 
5 3535 NEXT N 

3536 NEXT M 

3537 PRINT #1, "RPY" 

3538 FOR Q=l TO 9 STEP 1 

3539 FOR R=l TO 9 STEP 1 
10 3540 FOR U=l TO 9 STEP 1 

3541 LET H=RFY(U,R,Q) 

3542 LET RFYO(U,R)«H 

3543 NEXT U 

3544 NEXT R 

15 3545 MAT PRINT #1, RFYO, 

3546 NEXT Q 

3547 PRINT #1, "MIFY" 
3550 FOR Q=l TO 9 STEP 1 
3555 FOR R«l TO 9 STEP 1 

20 3560 FOR Ual TO 9 STEP 1 

3565 LET H=MIFY (U, R, Q) 

3566 LET MIFYO(0,R)=H 

3567 NEXT U 

3568 NEXT R 

25 3570 MAT PRINT »1,MIFY0, 
3575 NEXT Q 

3580! INVERSE TRANSFORM THE COLUMNS IN THE Z DIRECTION 
3581 DIM F(9,9,9) 
3590 DIM FO(9,9) 
30 3592 FOR V=l TO 9 STEP 1 

3593 FOR N=l TO 9 STEP 1 

3594 FOR MbI TO 9 STEP 1 
3600 LET Y=RFY(M,N,V) 
4546 LET R (M) «Y 

35 4547 NEXT M 

4548 FOR M=l TO 9 STEP 1 

4549 LET Y*MIFY{M,N, V) 

4550 LET I(M)=Y 
4552 NEXT M 

40 4553 EXTERNAL SUB IFTZ (DIM () ,DIM() ,DIM() ,DIM() ) 
4554 CALL IFT2 (R() , I {) ,MR() ,MI () ) 
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4556 FOR M«l TO 9 STEP 1 

4557 LET Y=MR(M) 
4559 LET P(M,N,V)=Y 
4600 NEXT M 

5 4602 NEXT N 

4604 NEXT V 

4605 PRINT #1, "F» 

4610 FOR Qol TO 9 STEP 1 

4620 FOR R=l TO 9 STEP 1 
10 4630 FOR U=l TO 9 STEP 1 

4635 LET H=P{U,R,Q) 

4640 LET FO(U,R)»H 

4650 NEXT U 

4660 NEXT R 
15 4666 MAT PRINT #1,P0, 

4670 NEXT Q 

4677 ! CORRECT FOR THE U(Z) CONVOLOTION AND NORMALIZE THE RECONSTRUCTION 

4678 DIM CF(9,9,9) 
4780 DIM CF0(9,9) 

20 4809 FOR V=l TO 9 STEP 1 

4810 FOR N=l TO 9 STEP 1 

4871 LET K=P(1,N,V) 

4873 K=K/10E8 

4974 LET CF(1,M,V)=K 
25 5915 NEXT N 

6616 NEXT V 

6617 FOR V=l TO 9 STEP 1 
6620 FOR M=2 TO 9 STEP 1 
6630 FOR N=l TO 9 STEP 1 

30 6640 LET K=P(M,N,V) 

6642 P=0.5*PI/(M-1)*2/C''2 
6650 K«K/P 

6660 LET CF(M,N,V)=K 
6670 NEXT N 
35 6680 NEXT M 
6681 NEXT V 

6690 PRINT #1, "RECONSTRUCTION" 

6691 FOR Q=l TO 9 STEP 1 

6692 FOR R=l TO 9 STEP 1 
40 6693 FOR U=l TO 9 STEP 1 

6694 LET H=CP(U,R,Q) 
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6695 H=H/CF(5,5,6) 

6696 LET CPO(U,R)=H 

6697 NEXT U 

6698 NEXT R 

5 6700 MAT PRINT #1,CP0, 
6705 NEXT Q 
6710 END 

6722 SUB FFT(R() ,MR() ,MI()) 
6726 FOR M=l TO 9 STEP 1 
10 6730 A=0 

6740 FOR N=l TO 9 STEP 1 
6750 LET H=R{N) 

6760 B»H*C0S(2*PI*(M-5)*(N-5)/9) 
6770 A»A+B 
15 6780 NEXT N 
6790 A=A/9 
680O LET MR(M}=A 
6810 NEXT M 

6820 FOR M»l TO 9 STEP 1 
20 6830 A»0 

6840 FOR N=l TO 9 STEP 1 
6880 LET H=R(N) 
6890 H=-H 

6900 B=H*SIN(2*PI*{M-5)*(N-5)/9) 
25 6910 A»A+B 

6920 NEXT N 

6930 A=A/9 

6935 LET MI (M) =A 

6940 NEXT M 
30 6950 END SUB 

6960 SUB IFT(R(),I(),MR{),MI(}} 

6970 DIM MRR(9) 

6980 DIM MRI(9) 

6990 DIM MIR (9) 
35 7000 DIM Mil (9) 

7010 FOR N=l TO 9 STEP 1 

7020 A«0 

7030 FOR M=l TO 9 STEP 1 
7040 LET G=R(M) 
40 7050 B=G*C0S(2*PI* (M-5)*(N-5)/9) 
7060 A=A+B 
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7070 NEXT M 

7080 LET MRR(N)=A 

7090 NEXT N 

7100 FOR Nal TO 9 STEP 1 
5 7110 A=0 

7120 FOR Mai TO 9 STEP 1 
7130 IJST OsR(M) 

7140 B»G*SIN(2*PI*(M-5)*{N-5)/9) 
7150 A«A+B 
10 7160 NEXT M 

7170 LET MRI(N)=A 
7180 NEXT M 

7190 FOR N»l TO 9 STEP 1 
7200 AaO 
15 7210 FOR M=l TO 9 STEP 1 

7220 LET G=I(M) 

7230 B=G*COS(2*PI*(M-5)*{N-5)/9) 
7240 A=A+B 
7250 NEXT M 
20 7260 LET MIR(N)<iA 
7270 NEXT N 

7280 FOR N«l TO 9 STEP 1 
7290 A=0 

7300 FOR M=l TO 9 STEP 1 
25 7310 LET G=I(M) 

7320 B=G*SIM{2*PI*(M-5)*(N-5)/9) 

7330 A=A+B 

7340 NEXT M 

7350 LET MII(N)=A 
30 7360 NEXT N 

7365 MAT MIIe(-l)*MII 

7375 MAT MR«MRR+MII 

7385 MAT MI=MIR+MRI 

7400 END SUB 
35 7410 SUB IFTZ(R() ,10 ,MR() ,MI{}) 

7420 DIM MRR(9) 

7430 DIM MRI(9) 

7440 DIM MIR (9) 

7450 DIM Mil (9) 
40 7460 FOR N=l TO 9 STEP 1 

7470 A=0 
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7480 FOR M»l TO 9 STEP 1 
7490 LET G=R(M) 

7500 B=G*COS(2*PI*{M-5)*(-N+l)/9) 
7510 A=A+B 
5 7520 NEXT M 

7530 LET MRR(N)=A 
7540 NEXT N 

7550 FOR N«l TO 9 STEP 1 
7560 A=0 
10 7570 FOR M=l TO 9 STEP 1 
7580 LET G»R(M) 

7590 B=G*SIN(2*PI* (M-5) * {"N+l)/9) 
7600 A»A+B 
7610 NEXT M 
15 7620 LET MRI(N)sA 
7630 NEXT N 

7640 FOR N=l TO 9 STEP 1 

7641 A=0 

7642 FOR M»l TO 9 STEP 1 
20 7650 LET G=I(M) 

7660 B=G*COS (2*PI* (M-5) * (-N+1) /9) 
7670 A=A+B 
7680 NEXT M 
7690 LET MIR(N)=A 
25 7700 NEXT N 

7710 FOR N»l TO 9 STEP 1 
7720 A=0 

7730 FOR M=l TO 9 STEP 1 

7740 LET G=I{M) 
30 7750 B=G*SIN(2*PI*(M-5}*(-N+l)/9) 

7760 AsA+B 

7770 NEXT M 

7780 LET MII(N)»A 

7790 NEXT N 
35 7795 MAT MII=(-1)*MII 

7800 MAT MR«MRR+MII 

7810 MAT MI«MRI+MIR 

7880 END SUB 

$ 

40 
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In an ^bodiment, the matrix inversion and/or the reiterative algorithms are used 
in combination with the Fouri^ Transfonn Algorithm. For example, an NMR image 
created by the Fourier Transfonn Algorithm may used as input for tiie first reiteration of 
the reiterative algorithm. 

5 

The Nyquist Theorem Wifh The Determination Of The Spatial Resolution 

The dmvation of Eq. (41) demonstrates that the system function behaves as a 
filtCT of the spectrum of the bulk magnetization function (NMR image). It is well known 

10 in the art of signal processing that if such a filter passes all fi:equencies for which an 
input function has significant energy, then the input function can be recovered 
completely fiom samples of the filtered function taken at the Nyquist rate. This premise 
embodies the Nyquist Sampling Theorem. The spectrum of the system function (Eq. (1)) 
is derived in APPENDIX V and shown in FIGURE Ic. This function is a band-pass for 

15 all fi'equencies of the bulk magnetization function where and are comparable. 
Thus, the bulk magnetization function can be recovered by sampling the continuous RF 
field function given by Eq. (40) at the Nyquist rate, twice the highest firequency of the 
bulk magnetization function, in each spatial dimension over the sample space for which 
the function has appreciable energy. Sampling operations other than the present 

20 operation and the negligible error encountered by not sampling over the entire sample 
space are discussed in McC. Siebert [11] and the references therein disclosed which are 
all incorporated herein by reference. In the absence of noise, the spectrum of the bulk 
magnetization function can be completely recovered if the detector spacing frequency is 
equal to the Nyquist rate which is twice the highest frequency of the bulk magnetization 

25 function, and this represents the limit of resolution. However, the density of the detector 
. spacing is limited by noise. The three-dimensional bulk magnetization map is a 
reconstruction from independent recordings at independent detector spatial locations 
relative to the voxels of the image space where two detector signals are independent if 
they are sufficiently spatially displaced from each other such that the difference in signal 

30 between the two detectors due to a test voxel is above the noise level. The resolution 
based on signal-to-noise arguments is discussed in the Contrast and Limits of Resolution 
Section. 
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Contrast and Limits of Relation 

S The ability to visualize a structure in a noise-iCree environment dq>endSy among 

other fectors, on the local contrast C, which is defined as 
•A/ 

C~ (75) 

where / is the average background intensity and AI is the intensity variation in the 
region of interest The main source of NMR image (also called magnetic resonance 

10 images (MBS)) contrast is 7, and 7, ^ch depend on tissue types. The contrast may be 
increased by using RF pulse sequences to polarize the protons which encode 7| 
(spin/lattice) and 7, (spin/^in) relaxation time information in the data of the secondary 
RF field in a manner straightforward to those skilled in the NMR art. 

Contrast, is not a fundamental limit on visualization since it can be artificially 

1 S ' enhanced by, for example, subtracting part of the background or raismg the intensity 
pattern to some power. Noise, however, represents a fimdamental limitation on the 
ability to visualize structures. The signal-to-noise ratio, a basic measure of visualization, 
is detennined by the ratio of the desired intensity variations to the random intensity 
variations, which are governed by the statistical properties of the system. The signal-to- 

20 noise ratio (SNR) is defined as 

SNR = ^=C^ (76) 

where 0, is the standard deviation of the background intensity representing the nns 
value of the intensity fluctuations. The noise properties of the 4D-MRI imager involve 

additive noise only principally &om thermal (y) noise in the RF measurement circuits of 

the antennas of the detector array and firom the fluctuations of the primary magnetic 
field. 

A feature of superconducting magnets is their extreme stability. The object to be 
imaged is magnetized with a highly stable magnet such as a superconducting magnet In 
this case, a magnetic field stability of 10^ % over a month's time is feasible. Small 
antennas measure the RF signals as point samples without significant decrease in the 
signal to noise ratio relative to large antennas by using impedance matching while 
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mmimizing resistive losses by using superconducting reactance elemGOts, for example. 
In an embodiment, cross talk between antemias is ameliorated or eliminated by time 
multiplexing the signal detection over the array of antennas. External sources of RF 
noise can be ameliorated by placing the 4D'MRI scanner in a shielded room (Faraday 
5 cage). 

The quality of the image (i.e. the signal to noise ratio of the image) can be 
increased by repeating the reconstruction over multiple time points wherein each data set 
of a given time point is the set of matrices of the intensity variation over the sample 
space of the RF field of the bulk magnetization M of each voxel having die 

0 corresponding magnetic moment at that synchronized (common) time point 

Sources of detector error are random, and the noise averages out as the rnunber of 
detectors increases. Typically, the noise is suppressed by a factor of the inverse square 
root of the number of detectors. Thus, the effective limit of RF field detection is 
increased by a factor of the square root of the number of detectors for a constant SNR. 

5 The resolution depends on the extent that the field of the ring, shell, or sphere of 

dipoles differs fiTom that of a single dipole at the center. The plot of the three cases of the 
field of a ring, shell, and a sphere of dipoles each of radius R and magnetic moment 
m = 10' Gem' given by Eq. ^14), Eq. (n.l7), and Eq. (IV.16) of APPENDIX I, 
APPENDIX n, and APPENDIX IV, respectively, as a function of radius R where the 

0 position of the center of the ring, shell, and sphere relative to the detector is the point 
(0,0,10) is given in FIGURES 2, 4, and 6, respectively. The plot of the three cases of the 
field of a ring, shell, and sphere of dipoles of radius R = 0.2 cm and magnetic moment 
/M = 10' Gem' given by Eq. (1.14), Eq. (n.l7), and Eq. ([V.16) of APPENDIX I, 
APPENDIX n, and APPENDIX IV, respectively, as a fimction of the distance between 

5 the detector at the origm and the center of the ring, shell, and sphere at the points (0,0, 
z = 4 cm /a z = 15 COT ) is shown in FIGURES 3, 5, and 7, respectively. From FIGURES 
2-7, it can be appreciated that the detector must be able to resolve three to four 
significant figures in order to reconstruct a ms^ of 0.2 cm resolution with a field depth of 
15 cm. 

3 

Finite Detector Length 



The system fimction, A , of the Reconstmction Algorithm Section is the impulse 
response for a point detector. The following analysis will concern the impulse response 
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for an NMR detector which has finite dimensions. Consider a detector airay 401 as 
shown in FIGURE 12 comprising multiple parallel planes 402 whmin each plane has a 
plurality of antennae coils 403 wherein each has an area A. The signal at any coil 403 is 
given by the integral of Eq. (40) over the area of the coil. This is the impulse response 
5 which is die syst^ function which replaces h for a finite area detector. The Fourier 
Transform of this system Amotion contains an argument of a product of the detector area 
A and the spatial firequency variables. Reconstruction could be perfonned as previously 
described in the Reconstruction Algorithm Section wliere this system fimction is 
substituted for the system fimction for a pomt detector. In the limit of zero area, the 

10 measurement is that of a point sample. Thus, another approach is to use the linearity of 
the superposition of RF signal to gain a higher signal to noise ratio advantage by taking 
the difference of relatively large signals as opposed to performing the measurements 
with miniature anteimas. In this case, the set of n coils each of area A of each plane 
comprise a grid of n blocks each of area a«A formed by the absence of common area 

1 5 overlap of one or more of the coils. The signal of each element of the grid of area a can 
be solved by the set of linear difference equations of the signals of coils each of area A 
that corresponds to the noncommon area for the overl^ping coils. The solution of the 
signals due to the grid elements can be obtain using a computer by matrix inversion. 
The resulting values represent the average signal for each grid center location. 

20 The effect of this data processing operation on the spectrum can be modeled as a sample 
and hold, vdiere the voltages at the centers of the grid elements are sampled by 
multiplying by a picket fence of delta functions spaced s units apart which are convolved 
with d , a square wave function in the x-direction and the y-direction of width s units 
where the coordinates of Eq. (1) and FIGURE 9 are used. In the firequency domain, this 

25 data processing operation causes the spectrum of the signal fimction s to be multiplied 
by D, the Fourier Transform of the square wave function of widdi s units, to form 
function S* , If this multipUcation does not multiply S, the Fourier Transform of the 
signal function, s , by zero for any firequency less than its bandwidth, then S can be 
recovered firom 5* by multiplying 5* with the inverse of the Fourier transform of the 

30 sample and hold square wave function, a two dimensional sine fimction for the x and y- 
directions. This analysis applies to all axes in which direction the detectors have finite 
length. Furthermore, z-sampling is achieved by translating the array in the z-direction by 
interval distances at which points discrete signals are recorded or by using multiple 
parallel plane detector arrays spaced at the sampling interval along the z-axis. However, 
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if llie signals are not sampled at discrete z-points, but each sample point is the integral 
resultant of the signal acquired continuously over a z-displacement of q units i^ch is 
much greater than the dimension of the detector in the z*direction» then the 
corre^nding sample and hold square wave has a widfii of q units. 

5 

It will thus be seen that the invention efficiently attains the objects set forth 
above, among those made apparent from the preceding description. Since certain 
changes may made in the above constmctions without departing from the scope of the 
invention, it is intended that all matter contained in the above description or shown in the 
10 accompanying drawings be interpreted as illustrative and not in a limiting sense. In 

particular, other methods such as electronic and optical methods of detection of magnetic 
resonance are within the scope of the preset invention. 

While the claimed invention has been described in detail and with reference to 
15 specific embodiments thereof, it will be apparent to one of ordinary skill in the art that 
varioiis changes and modifications can be made to the claimed invention without 
departing from the spirit and scope thereof. 
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APPENDIX I 

Integration of a RING of Dipoles 

Derivation of the Field Produced by a Ring of Magnetic Dipoles 

The z -component of the magnetic field due to a dipole or a loop of current of 
radius R with dipole moment m = in^ at the origin is given as follows: 



10 B =m 



(LI) 



The z -component of the magnetic field at the position (x,y,z) due to a ring of dipoles of 
radius R with dipole density cmtered at the ori^ is 

27C/? 



_ ;„ ^ '(2(z-zf-(x-.-y-Cy-yyV 
' 2«*i ((x-;Cy+(y-yy + (z-fyf 

wherein (x*,/^*) is a variable which conesponds to die position of each individual dipole 
IS of the ring. The relationship between Cartesian cooixlinates and cylindricd(X)ordinates 
wifli f =0 is 



;i' = ^cos^ a-3) 
y = i{sin^ 

Substitation of Eq. (L3) into Eq. (1.2) ^ves 

'^(2Z^ -(X-/?COS^)'-0>-Sin»)')/^ 

Bi =:r- 7 ; ; — w 

2 « 0 ((x - i?cos 0)' + (y - iJsin + z^j 

20 Multiplying out terms gjves 

^ ^ '(2z^ - (c^ + CPS' 0 - 2x/?cos^)- -t-i;^ sin^ ^ - 2;<Rsin^))/^ ^^^^ 
'~2>rj + cos'^-2xRcos^)+^' + /iW*-2>'Jlsin0)+z'y" 

Siibstitution of Eq. (13) in the denominator and using the assodative relationship and the 
trigonom^c identity, cos* a + sin* a = 1 in the denominator gives 

n ^^(2 z* -(t* + /?'cos^^-2xgcos0)-^ +^^sin^0 -2y^sin^))f0 

2 w J (r* + - 2;{(xcos ^ + ysm 0)) 
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The doiominator can be qipioximated using tbe following relationship 
(a + 6)"sa'+«6a"-* 

whoe 

a = r' + je 

b - -2i?(x cos ^ + ^'sin ^) 

and 

n = -5/2 

This case gives the far field where a«b and/or r»R. Usmg Eqs. (I.7-I.10) gives 

p + ^ -2ii(xcos^ +yw^^^"^ 

1 5i? 



a7) 

(1.8) 

a?) 
a-io) 



10 Substitution of Eq. (LI 1) into Eq. (L6) gives 



B. =- 



... J?'cos'^ 2xficos^ 

2a J. Si^Vcos^^sin^ lO;c^j?^cos^0 IQjcyi^^cos^sin^ 
w' f / . R^sux^<t> 2yRsm(l> 

2n I . SJfysin^^ lOJt^.xycos0sin^ 10^/sin^^ i 



(LIl) 



(L12) 



The integration of Eq. (1.12) gives 

f 4flz' 2 roc' wg' lOrtx'j?^ 1 



B,=— ^ 

2n i 2gy' 7tR' lOw^V 



I 

J 



(1.13) 
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Collecting tenns gives 



5 

Ring of Dipoles: 




a.i4) 



10 : 

Eq. (1.14) d^onstrates that the magnetic field of a ring of dipoles is not equal to 
that of a single dipole at the origin. The ring radius, R , spears in the denominator of 
eachtenn. The first term is the magnetic field of a dipole at the origin only when the 
15 variablecorrespondmgtotfaeradiusof tfaering, A,iszm. The second and Aird terms 
are additional perturbations of the field of a dipole at the origin whose magnitude is a 
fimction of the radius of the ring. The second and third terms vanish only when the 
radius of the ring is zero. A ring of dipoles has a field that is cylindrically symmetrical. 
A shell, a cylinder, and a sphere of dipoles are the only other cases which have this 
synmietry. Acylinder is a linear combination of rings. Thus, the uniqueness of the 
dipole field is demonstrated by showing that it is different fix>m that of a ring, a shell, and 
a sphere. The present result that the field of a dipole is different for that of a ring of 
dipoles as well as the same result in the cases of a shell and a sphere of dipoles shown in 
APPENDIX n and APPENDIX IV, respectively, demonstrate that dipole field is unique. 
All other fields are a linear combination of dipoles. Thus, the dipole is a basis element 
for the reconstruction of a NMR image. The resolution depends on the extent that the 
field of the ring of dipoles differs fit>m that of a single dipole at the origm. The plot of 
the field of a ring of dipoles of radius R and magnetic moment m^lO^ Gem given by 
Eq. (1. 14) as a function of radius R is given in FIGURE 2. The position of the center of 
the ring relative to the detector is the point (0,0,10). The plot of the field of a ring of 
dipoles of radius R = 0.2 cm and magnetic moment m = 10^ Gem given by Eq. (1. 14) 
as a function of the distance between the detector at Uie origin and the center of the ring 
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at ttie points (OA z=4 ciRroz = lScm)isshowninFIGUSE3. From FIGURES 2 and 
3 it can be q>preciated that the detector must be able to resolve three to four significant 
figures in order to reconstruct a map of 0.2 cm resolution with a field depth of IS cm. 
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APPENDIX n 

Dategration of a Spherical Shell of Dipoles 

Dcaivation of Ihe Field Produced by a Shell of Magnetic Dipoles 

The z -component of the magnetic field due to a dipole or a loop of current of 
radius R with dipole moment m = iTc^ at the origin is given as follows: 



10 B,=/B 



OLl) 



The i -component of the magnetic field at the position (x,y^) due to a shell of dq)oles of 

radius R vnSn dipole density entered at the onan is 

4kR^ 

wherein (x'.y^*) is a variable which corresponds to the position of each individual dipole 
IS offliesheU. The relatioiiship between C^esian coonlinates and spherical coordinates is 

X =psm^cosO jc'= iJsin^cos0 (IL3) 

y=psin^sin0 y = 5sin^sin0 

Substitution of Eq. (IL3) into Eq. (112) gives 

^ m y|> ^(z-i?cos<^y -(x-jgsin^cos0y -(y- jgsin^ sine)')??' sin ^rf^rfO 
0 0 ((x~i?sin^cose)' + ();-/?sin^sin0)V(2-iJcos^/J' 

(n.4) 

20 Multiplying out terms gives 

g m Y( ~ 4zRcos ^ + 2R^ cos^ 4>-x^+ 2xR sin ^ cos 0 

' 4;ri? J J(t^--2xi?sin^cos0 + i?sin^^cos^0+/~2yi?sin^sin0 

-JZ^sin^ ^ cos^ 0 + 2>/;tsin 0siQ0 - i?^ sin^ ^ sin^ 0)/? sin 
+iP sin^ ^ sin^ 0 + - 2zi?cos ^ + ii^ cos^ ^ 

(n.5) 

Multiplying out the j?^ sin ^ term gives 
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g m ^? f (2 jg^z' sin ^ - 4zR^ sin ^cos^ 2R^ cos^ ^ sin^ - x^^ sin^ -f 2jcR^ sin' ^ cos Q 
'''47C^ll ^ (c^-2jcRsin^cose+i?^sin^^cos^0+/ -2yRsin^sine 

-/t' sin' ^cos' e-y sin^ + 2yR^ sin' ^sin Q - it' sin' 0sin' e)/^rfe 
+J?'sin'0sin*0+z' -2zi!cos^ + i?' cos'^jf'* 

(IL6) 

Substitution of Eq. (11.3) in ttie denominator gives 

^ in '? j(2ig^z'sin ^ - 4^^^ sin ^^cos^ + 2R^ cos' ^ sin^ - jt^R? sin^ + 2x/g' sin' ^ cosO 
'^4;ii? (p' -2xRsin^cos0 + i?sin'^cos'0-2>'i?sin^sin0 

-R^ sin' ^cos' 9-/^ sin^ -hlyR^ sin' ^sin g - J?' sin' ^sin' e)/^rfO 
+ii'sin' ^sin' 0 - 2zJl cos^ + /{' cos' ^ 
5 (n.7) 

Using the associative relationship and using the trigonometric identity, 

cos'a + sin' a = 1 in the d^ominator gives 

m '? f ^ j?'z'sin ^ - 42i?'sin ^cos^ + 2R^ cos' ^ sin0 - x^I^ sin^ + 2jcig^ sin' <^ cos 6 
' j] (>' +/?' -2xRsin*cos0 

-/g^ sin' ^cos' e-/Ii? sin^ + 2yR' sin' ^sin 0 -R^ sin' ^sin' e)/^Je 
-2>i{sin^sin6-22i{cos^) ^ 

(n.8) 

1 0 The denonunator can be approximated using the following relationship 

(a + &)"sa"+/i6a""' (n.9) 

where 

a-p' + /?' (n.io) 
6 = -2xR sm ^ cos 0 - 2yR sin^sin0-2ziicos0 (n.11) 

15 and 

« = -5/2 (n.l2) 
This case gives the far field where a«b and/or p»R, Using Eqs. (11.8-11. 12) gives 

(p' +i?' -2xRsin^cos0-2>;i?sin^sin0-2zficos0)'^'^ 

1 5/? , . , ^ . ^ . ^ ^x(n.l3) 

(jcsm ^cos 0 + ysm 0 sm 0 + zcos^) 



\7/2 y 



Substitution of Eq. (11.13) into Eq. (IL8) gives 

« 'f ff2/?^z'sin^ lOi^z'sin^. . ^ „ .t • « a\ 

' J J| ^^2 ^ + ^ (Jcsm^ COS0 +ysm^ sin0 +zcos ^) 
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4z^smdicosA 20z/{^sin6cos^y . , . • . . « i\ 

- . , ^v5>i^ / 1 ^ (xsm^cose+j'sm^sme + zcos^) 

2Itcos^6sm6 lOR'co^isinir . . « • ^ • « .\ 
(p+tr) (p 

- . . .J/2 - / , l,J)i (jgsin^cose^-ysui^sine + zcos^) 

+ ^2 4. jj^y" (p^ + ie^'^ — (jcsin^cosO+ysm^smO+zcos^) 

i?*sin'^cos^0 5/?' sin' ^ cos* 0, . . „ • a • « a\ 
— (p* + /pj"^ ^-j-^;y72— (*sm^cose+ysm^sm0+zcos^) 

"71 — ;SW~71 — ;T777iVcsm^cos9+;/sin^sm0 + zcos0) 

2 v/? sin* A sine I0v/?*sin* Asin0/ . . « .1 • /> a\ 
+ ^^^^yh +— ^T-j-^yTT— (^sin^cose+ysin^sine +zcos0) 

i?^sin'isin*0 5i?' sin' ^ sin* 0 , . , . • . • « ^Jj^jo /iti>i\ 
— (p* + ^y" (p* (Jtsin^cose +;>sm^sme + zcos^)|/^f/0 (n.l4) 



10 The integral of a sum is equal to the sum of the integrals. 



Inn] 
0 I 



J J / 2' ,.2\^^^ •*• / 2 „,v^2 (Jgsm0cose+;;sm^sme+zcos») 

7fr4zR'sin^cos^ 20z/{*sin^cos^/ . , . ■ , ■ n xVjxjq 
" j; [ (p^ + J^f (p* + f^y" (xsm^cose +j/sm^sm0 + zcos^) #t/0 

'f"(j 2/t*cos* ^sin^ 10/?'cos*^sinA X . . • a • a, x\ jajo 

J J ~71 — +— 71 — ;^b7r^(rsm#cos0+3;sm^sm0+zcos^) d^O 

iiy (p'+R') (p'+^) J 

'f'J x^R^anA Sx^^amk ^ . , ^ . , . « j\ ,« 

^fV 2x/?'sin* Acos0 1 Ox/?* sin* A cos 0, . . . ■ . • « jl\ 
15 +11—7-^ — r.s/2 + — 71 — ;t^772 — (rsm^cos0+j/sm^sm0+zcos^) 

*ff r/?*sin' ^cos*0 5/?' sin* ^ cos* 0. . . . • x • « A^ jajo 
~Jj /2 r.2\i^^ + /2 (rsm^cos0+ysm0sm0+zcos^)U/^0 
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(x sin^ cos d + j'sin ^sind + z cos ^) </^d 



11 

775— (jcsin^cosfl +>'sin^sm0 + zcos^) Mdd}^ 

i J 



(IL15) 



The integration of Eq. (IL15) gives 



B = 



\7/2 + 



^577+ 



772 



' 4n^[(p' + ief 3(p'+R')" 3(p'+R'f ^^+^ 3^7^ 

(n.i6) 



10 (Collecting terms gives 



Shell of Dipoles: 

15 



20 Eq. (n. 1 7) demonstrates that the magnetic field of a shell of dipoles is not equal 

to that of a single dipole at the origin. The shell radius, R , appears in the denominator of 
the first term. The first term is the magnetic field of a dipole at the origin only when the 
variable corresponding to the radius of the shell, R , is zero. The second term is an 
additional perturbation of the field of a dipole at the origin whose magnitude is a 

25 fimction of the radius of the shell. The second term vanishes only when the radius of the 
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shell is zero. Thus, the dipole is a basis element for the reconstruction of a NMR image. 
The resolution depends on the extent that the field of the shell of dipoles differs fiom ttiat 
of a single dipole at the origh. The plot of the field of a shell of dipoles of radius R and 
magnetic momrait m = 10* Gcm^ ^ven by Eq. (EL. 17) as a function of radius R is given 

S in FIGURE 4. The position of the center of the shell relative to the detector is the point 
(0,0,10). The plot of the field of a shell of dipoles of radius = 0.2 cm and magnetic 
moment m = 10* Gcm^ ffvea by Eq. (11.17) as a fimction of the distance between tiie 
detector at the origin and the center of the shell at the points (0,0, z) is shown in FIGURE 
5. From FIGURES 4 and 5 it can be appreciated that the detector must be able to resolve 

10 three to four significant figures in order to reconstruct a map of 0.2 cm resolution with a 
field depth of IS cm. 
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APPENDIX m 

Proof that the Field Produced by a Shell of Magnetic Dipoles is Different from ttiat of a 
S Single Dipole 

Consider Eq. (11.5). 

g m ^rr (22^ -4z/tcos^ -h2iZ'cos^ ^-x^ -f 2xgsin^cos0 
'"47ri? jJ(x'-2jcRsin^cosO + iesin^^cos'0+/~2>!«sm^sin^ 

-ii^sin* 0 cos^ 9 + 2>;i{sin 0sin0 ~ sin^ 0sin^ e)t^ sin^^rffl 
+1? sin^ ^ sin^ 8 + - 2ziicos^ + cos^ 

(in.l) 

10 The conditions for this integral to equal that of the field of a dipole at the origin are 

-4z/?cos^ + 2i2^ cos^ ^ +2x;?sin^cos0 - ii^sin^ ^cos^ 0+2j;J?sin^ sin0 -ii^sin^ ^sin^O =0 

(in.2) 

and 

-2xi?sin<^ COS0 + sin^ 0 cos^ 0 - 2yi?sin ^sin0 + i?^sin^ ^sin^ 0 - 2zi?cos ^ + R^ cos^ ^ = 0 
15 

(ffl.3) 

Thus, Eq. (in.2) must equal Eq. (111.3). 
-4zifcos^ + 2i? cos^ ^ +2xi?sin^cos0 - 1? sin^ ^cos^ 0+ 2yifsin^ sinO - 1? sin' 0sin'0 
= -2x/isin^cos0 + ii'sin' ^ cos' 0 - 27i?sin 0sin0 + if sin' 0sin' 0 - 22i?cos^ + cos' ^ 

20 (ni.4) 
Therefore, 

-2zRcos^ + cos^ ^ - 2/?^ sin^ ^ cos' 0 + 4xsin^ cos0 + 43/i?sin^sin 0 - 2if^ sin^ 0sin' 0=0 

001.5) 

25 Using the associative relationship and the trigonometric identity, cos' a + sin' a = 1 as 
well as factoring out an in Eq. (III.5) gives 



30 



-22i?cos^ + cos' - 2*' sin' ^ + 4xsin ^cos0 + 4yRsm(^ sin0 = 0 (III.6) 
-2zcos 0 + 3/? cos^ 0 - 2/? + 4x sin^ cos 0 + 4ysm ^ sin 0 = 0 (HI.?) 
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For X9fe0,>'9ft0, z^ftO, Eq. (HLT) is true only if 15 = 0 which proves that the field of a 
shell of dipoles centered on the oiig^ is differrat from that of a smgle dipole at ttie 
origia 

Eq. (mS) can not be integrated in closed form; however, the integral can be 
5 approximated so that the detector tolerances for a given image resolution can be 
determined. The denominator of Eq. (ULS) can be ^^proximated using the following 
relationship 

(a+by^a^^nba''^' (m.8) 

where 

10 a^p^ + R^ (m.9) 

b = -2x/?sin ^cosO - lyRsin^ sin 0 - 2zi2cos 0 (III.10) 

and 

rt = -5/2 (in.ii) 
The result given by Eq. (IL17) is 



15 B,=-J^^4nJJ^ 



(ni.i2) 



Higher ord^ terms of the s^>proximation given by Eq. (IIL8) would contain mixed 
products of the coordinate variables which would increase the deviation of the derivation 
of the field a shell of dipoles from that of a single dipole at the origin. 
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APPENDIX IV 
Megratioii of a Sphere of Dipoles 
Derivation of fhe Field Produced by a Sphere of Magnetic Dipoles 



The z -componrat of the magnetic field due to a sphere of dipoles is derived 
using the equation for a shell of dipoles with die substitution of die dipole density of the 
3/w 

10 sphere — where the radius of the sphere is R . The field of a shell of dipoles is given 

471K 

^proximately by Eq. (HIT) of ttie Derivation of the Field Produced by a Shell of 
Magnetic Dipoles Section 



B = 



Anr 



AOnr* 



2z^ -x^ 



(IV.l) 



where the radius of the shell is r (The parameter r replaces the constant R of Eq. (11.17) 
15 as the radius of the shell.). For small r , the first term of Eq, (IV.l) dominates and is 

used to approximately calculate the z -component of the magnetic field due to a sphere 

of i -oriented dipoles as follows: 

A sphere of dipoles is equivalent to the integral over concentric shells of dipoles 

each of radius r where 0^r<R. Thus, the z -component of the magnetic field due to a 
20 sphere of dipoles is the integral of the field of the shells given by Eq. (IV.l). 



B = 



m 



\r^dr 



Let 



25 



- = tan0;rfr = psec^ecfe 
P 

and using the trigonometric identity 

l + tan^0=sec^0 
Eq. (TV.2) becomes 

« 2 2 2\ 1* 47itan^0 2nj^ 

3 



(IV.2) 



(IV.3) 



(IV.4) 



(rv.5) 



Substitution in of 
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and 



tan0 = 



sece = 



sinO 

COS0 

1 



COS0 

intoEq. (TV.S) gives 



3 



Let 



= sin 9; = cosOJd 
Substitution of Eq. (IV.9) into Eq. (IV.8) followed by integration gives 

10 Substitution of Eq. (IV.9) into Eq. (IV.IO) ^es 

3 



FromEq.(IV.3) 



— = tanO; .sinQ = 
P 



Substitution of Eq. (IV. 12) into Eq. (IV. 11) gives 



15 



B,=(2z^-x^-/) 



m 



1 + 



m 



Evaluation at the linuts of the integral gives 



B,=(2z^-;c^-/) 



m 



11+ 



Multiplication of Eq. (IV. 14) by 1 = - 



37210 



1 + 



1 + 



gives 



(IV.6) 

Civ.7) 

(IV.8) 

av.9) 

(FV.IO) 



(IV.ll) 



(IV.12) 



(IV.13) 



(IV.14) 
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10 



Multiplying out the term in the denominator gives 




Sphere of Dipoles: 



(WAS) 



(IV.16) 



Eq. (IV.16) demonstrates that the magnetic field of a sphere of dipoles is not 
equal to that of a single dipole at the origin. The shell radius, R , appears in Hhe 
15 denominator of the first tenn. The first term is the magnetic field of a dipole at the origin 
only when the variable corresponding to the radius of the shell, R , is zero. The second 
term is an additional perturbation of the field of a dipole at the origin whose magnitude is 
a fimction of the radius of ttie qihere. The second term vanishes only when the radius of 
the sphere is zero. Thus, the dipole is a basis element for the reconstruction of a NMR 
I image. The resolution depends on the eictent that die field of the sphere of dipoles differs 
&om that of a single dipole at the origm. The plot of the field of a sphere of dipoles of 
radius R and magnetic moment m = 10^ Gcrn given by Eq. (IV. 16) as a fimction of 
radius R is given in FIGURE 6. The position of the center of the sphere relative to the 
detector is the point (0,0,10). The plot of the field of a sphere of dipoles of radius 
R=^02 cm and magnetic moment m = 10^ Gcm^ given by Eq. (IV.16) as a fimction of 
the distance between the detector at the origin and the center of the sphere at the points 
(0,0, z) is shown in FIGURE 7. From FIGURES 6 and 7 it can be appreciated that the 
detector must be able to resolve three to four significant figures in order to reconstruct a 
map of 0.2 cm resolution with a field dq>th of 15 cm. 
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Fourior Transform of tbe System Function 



5 The system function, h(p,^,z),ia cylindrical coordinates is 

The spacetime Fourier transfoim in three dimensions in cylindrical cooidinates, 
^^(*p» is given [9] as follows: 

a» 2jv<o 

Hik^,<I>A) = J J /A(p.^,z)exp(-,27r[A:ppcos(0- ^)+ K^'^pd^dz (V.2) 

-flO 0 0 

10 With circular symmetry [9] 

=2;r/ \h(p,z)J,{}c^)e-'"pdpdz (V.3) 
o-« 

The Fourier transform of the system function is given by the substitution of Eq. (V. 1) 
into Eq. (V.3). 

L^^i ^f^^?f2 Jo[k^]pdpe^''dz (V.4) 

1 5 Consider the integral of Eq. (V.4) with respect to dp only. Factorization of A(p, ^, z) 
gives 

Consider the definite integral 

^ [/'+z'r* 2T[« + 1] ' 
20 and the modified Bessel fiinction of the third kind relationship, 

K^,M^K,[x] (V.7) 

3 

The first factor of Eq, (V.5) is the same form as Eq. (V.6) with v = 0; « = - , flius, 

(V.8) 
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where iCa/iCVl = ^viiK^] (Eq. (V.7)). The second factor of Eq. (V.5) can be made 
mto fhe same form as Eq. (V.6) using the Bessel function of the first kind recurrence 
relationship 



S Consider fhe second factor of the integral of Eq. (V.S). 
Eq.(y.9)witho = lis 

10 Let 

Substitution of Eq. (V.13) into Eq. (V.12) is 



Substitution of Eq. (V.IO) into Eq. (V.14) is 



(V.9) 

(V.IO) 

(V.ll) 
(V.12) 

(V.13) 

(V.14) 



15 



(V.15) 

The first factor of the right hand side of Eq. CV.15) is the same fonn as Eq. CV.6) wifli 
3 

v=l; « = -,thus. 



(V.16) 

where K.y^lkz] = Kti^lkpZ] (Eq.(V.7)). The second fector of the right hand side of Eq. 



(V.15) is the same fonn as Eq. (V.6) with v= 2; « = — , thus. 
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(V.17) 

Combining the parts of the integration wiUi r^pect to dp of Eq. (V.4) by adding Eq. 
(V.8), Eq. (V.16), and Eq. (V.17) gives 

£1 LnF^^'f*'^^-S^ wi^rF^^ 



1/2 f 3/2 



The modified Bessel function of the third kind formiilae is 
Substitution of Eq. (V.13) into Eq. (V,19) with u = 1 is 



10 



L2V. 



1/2 



Substitution of Eq. (V.13) into Eq. (V.19) with o = 0 is 
Substitution of Eq. (V.20) and Eq. (V.21) into Eq. (V.18) is 

r 



(V.18) 
(V.19) 

(V.20) 
(V.21) 



£1 



|[ r[5/2] ['^2k^2 ^^J r[5/2] "'(2''^)r[5/2]X2A:,zJ j 
15 (V.22) 



^{r[5/2] " r[5/2]2 



r[5/2] 



r[5/2p " J 



20 Collecting terms gives 

With r[3] =2 and r[5/2]=3 /4«"', Eq. (V.24) is 



(V.23) 
(V.24) 

(V.25) 
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2nk^[f'''e-^''dz (V^7) 

4nk^ ^e'^V^^'dz (V^8) 

Ank^^e'^-^^'^dz (V.29) 
S Integration of Eq. (V.29) with respect to dz gives 

Multiplication of Eq. (V.3 1) by 
10 gives 



(V.32) 



4»t.^-^j (V.33) 

The system function (Eq. (V.l)) is an even function; thus, the spacetime Fourier 
transform in three dimensions in cylindrical coordinates, H{k^^k^)y is given by taking the 

real part of Eq. (V.33) [8]. 
15 H[k^,k^ = J^ (V.34) 

The spacetime Fourier transform in three dimensions in Cartesian coordinates, H{kpyk^)y 
is 

where the relationship betwera the wavenumbers and the spatial Cartesian coordinates is 
20 as follows: 
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2jc 2n 



(V.38) 
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APPENDIX VI 
Derivation of 5 = HF® uQc,) from Eq. (55). 



*.»0 ^— h/2 X,«-4/2 (VII) 

y/c, 'V ^ J'.— '2/2 X,— 

where the Fourier transform of m(z) = 1/or z S: 0 aw^/ «(z)= \forz<0 [8] is 

10 and where H is given by Eq. (52) and F is given by Eq. (49). The convolution integral 
of the second term of Eq. (VI. 1) is 

Collecting exponential tarns gives 

2 ,,»o J'.— la/2 M2 J — ^xpp '^^J J 

15 (VI.4) 
Expansion of the draominator of the second term gives 

(VI.5) 

20 The factorization of the denominator of the second term, kI - 2k^K^ + ^ + using the 
quadratic formula is 
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2 

^^K^!t!t^iz^ (VL6) 
Substitution of Eq. (VL6) into Eq. (yi.5) is 

5 (VI.7) 
The expansion of Eq. (VI.7) by the method of partial fractions is 



2 ,,=0 J'.'-IJ2 ».=-/,/2 



— e. I — + + dK. 



(V18) 



The factors A, B, C of Eq. (VL8) are detomined as follows: 

4K*-2^fc,+A:;+*^]+4K,Jc.-^-y*J+C[K,f:.-^+y*,]=l (VI.9) 
10 Let K, = 0. Then, Eq. (VI.9) gjves 

A^l+k^yi (VLIO) 

Thus, 

Let K^ = k,+ jkp. Then Eq. (VI.9) is 

-2^ft, + +y*p-^ -y*p] 

+ Cfi Jt, +yAp ]= 1 

The term in B vanishes. Substitution of Eq. (VI. 11) into Eq. (VI. 12) gives 
l^^p +2yAA -K -^l - Wp ^kiycJ^jkJc^ -2kl]= 1 (VLB) 

where the term in A also vanishes. Solving for C gives 

C=r T-^ i CVI-14) 

[-2*^ + 2^*,] 

20 Let = jkp . Then, Eq. (VI.9) is 



(VL12) 
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The torn in C vanishes. Substitution of Eq. (VLll) into Eq. (V115) gives 

v^eie the term in A also vanishes. Solving for B gives 

5 B- ^ 

\-2k^-2jkJc^] 

Substitution of , , and C into the convolution integral (Eq. (VL8) gives: 



*,oO y«"-/a/2 x«— 1|/2 



7^1 [K,-*.+y*pI-2*;-2Mf'''i 



The first convolution integral is of the fonn given by Mc Seibert [8]. 



— o sgnz', sgnz„ = -1 where z„ <0 
JK 



10 Thus, 



rfK. = —rn — rr=i where <0 



A change of variable in the second convolution integral given by letting 
K, = -k^ -jk^ ; rfK, = -dk^ 

in 



15 



z!l52 f — £ 



IS 



r P i f ^-F T-MJ where z„ <( 



(VL15) 



(VL16) 



(VI.17) 



(VL18) 



(VI.19) 



(VI.20) 



(VL21) 



(VI.22) 



(VI.23) 
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1 J ^-r-dk, (VI.24) 



= -rr^ where <0 (VL25) 

A change of variable in the thud convolution integral givoi by letting 
5 K, = -K +J% . dK, = -dk, (VL26) 

in 



IS 



r , " i f ^^^-F T-[-<«^] where z, <( 



10 



[-2A^ +2yAit,].i 2yt. 



nk. 



2 



(VL27) 



(VL28) 



^ -V'^''-' where z. <0 (VI.30) 



15 



-JkJ^p] 

Combining Eqs. (VI.20), (VL25), and (VL30) gives the convolution of Eq. (VLl). 

, -/i +/j/2 *ij2 T. . ^ -T— SirA* 



(VL31) 



wo 02/16956 



107 



PCTAJS01/2S954 



APPENDIX Vn 



Derivation of the Inverse Transfonn of Eq. (68) to Give Inverse Transfonn 1, Eq. (69). 



(VILl) 



10 » ». J-. *. 



where <0 

Consider the definite integral: 



-00 0 I 

-» 0 J 

(VIL2) 



(VIL3) 



In the case of the first integral of Eq. (VIL2), the parameters of Eq. (VII.3) are 

v = 0 



15 



b = p 



(VIL4) 



Thus, 



where 



2kirry2 



(VIL5) 



(vn.6) 
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Consider the following relationship of die modified Bessel function of flie third kind: 

7Lr^[z]=e 2 K^^e ^j + c ^ iSC^j^ze (VIL7) 
In the case that v = 0, Eq. (VII.7) is 

5 Thus, 

0 0 

0 ** 0 . 

(vn.io) 

Witib a change of variable as foUows: 
10 First 
Integral: Second Integral: 

and an expansion of the complex exponential factor, Eq. (VILIO) becomes 



= T 1 [cosV.> jsinft,2.]f cosfk,z,]-ysin[v. J:.[*,P>p (Vn.12) 



0 



15 The transform is given by Bateman [12]. 
Thus, 



(vn.ii) 



(vn.i3) 



Je*>-'y,^,p]*,=p-:^ (vn.i4) 

0 L^n +P J 



0 

The Fourier transform relationship between a function /[x] and its Fouria: transform 
20 g[y] given by McC. Siebert [8] is 
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^^]« (vn.i5) 

J ok 

Consider flie following Inverse Fourier transform integral fiom Eq. (VIL2) 

0 

Eq. (Vn.l6) is equivalent to the following Inverse Fourier transform 
From Eq.(VII.15) with 

Eq. (Vn.16) is 



(vn.i8) 



(vn.i9) 



10 Taking the imginaiy part because Eq.(VIL 19) has odd sjvnmetrygive^ 

^'^''■'y:,|^^"kfe„ = TT^N^ ; z„ <o (VIL20) 



0 

In the case of the third and fourth terms of Eq. (Vn.2), consider the following Fourier 
transform given by Gray [13]: 

Je-"y,[6x]c&=p-^ (VIL21) 

15 Thus, 

J ^''^V.fir^Jt, =Y7?f'' ^ ^^^^ 

The Inverse Fourier transform with respect to is given by McC. Siebert [8], In the case 
of the first two terms of Eq. (Vn.2), 

eQ 

Je*^'<fi^=5[z] (VIL23) 

— CD 

20 In the case of the third and fourth terms of Eq. (Vn.2), 

CO 

\ jkfi*^'-'dk,^h[z] (vn.24) 
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Combining transfonns fiom Eqs. (VIL6), (VII.20), and (VIL23) in the case of the first 
and second terms of Eq. (Vn.2) and combining transforms fit)m Eqs. (Vn.l4) (VIL22) 
and (Vn.24) in the case of flie third and forath terms of Eq. (Vn.2) gives: 

l'»».-o ^.-«,p L 'v^ ilJ 



kl kl 



whwe ^ <0 

Substitution of Eq. (48) into Eq. (Vn.25) gives 



(VIL25) 



r' 



<>x.»0 j..o-y2 ;«,=^</2 L 'V*" "^IJ 



8 1,^.0 ifj=-/i/2* ii,«-/,/2t 



10 



where <0 

(Vn.26) 

Taking the derivative given by the doublet fiinction gives 



l<>,,-0 jii«^ft/2 jr. --4/2 



15 



^ijfk ^yflk +Ixt2k 



Inverse Transform 1 



(Vn.27) 
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1 . Apparatus for providing a multidimensional image of inanimate or animate objects 
conqprising: 

S a magnetization source for producing a magnetic field to magnetize a volume of 

an object to be imaged, the volume comprising a plurality of voxels; 

a radiation source for applying a radiation field to the object to be imaged when 

the object is disposed in the magnetic field to produce a secondary radiation field 

emanating bom the object; 
0 an output signal detector for producing output signals in response to said 

secondary radiation at a plurality of spatial locations outside of flie object as a 

function of time; 

a first processor for determining a plurality of Fourier components, each having 
the same fi-equency, an intensity and a phase angle; 
5 a second processor for associating the Fomier components due to each voxel of 

the object to be imaged by phase to form a set corresponding to a spatial variation of 
intensity of the secondary radiation due to each voxel at said plurality of spatial 
locations outside of the object; 

a third processor for converting each said set of components into a voxel 
0 location; and 

an image processor for producing an image based on the secondary radiation 
intensity firom each voxel and the location of each voxel. 

2. The apparatus according to claim 1, further comprising a display for displaying the 
5 image. 

3. The apparatus according to claim 2, characterized in that said display comprises at 
least one selected firom the group consisting of printers, cathode ray tube displays, 
liquid crystal displays, plasma screens, three dimensional modelers, laser monitors, 

0 projections monitors and holographic displays. 

4. The apparatus according to any one of the preceding claims, characterized in that one 
processor handles the processing functions of at least two of said first processor, 
second processor, third processor and image processor. 
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S. The apparatus according to any one of the preceding claims, characterized in that one 
processor handles the processing functions of at least three of said iSrst processor, 
second processor, third processor and image processor. 

5 6. The apparatus according to any one ofthe preceding claiins, characterized in that the 
third processor comprises a reconstruction processor. 

7. The ^aratus according to any one of the preceding claims, characterized in that the 
image is a magnetic resonance image. 

0 

8. The apparatus according to any one of the preceding claims, charactCTzed in that the 
magnetization source comprises a superconducting magnet 

9. The ^paratus according to any one of the preceding claims, characterized in that the 
5 radiation source comprises a radio frequency generator and transmitter. 



0 



10. The ^aratus according to claim 9, characterized in that the radio frequency 
generator and transmitter comprises an antennae to excite protons in a magnetized 
volume. 

1 1. The apparatus according to claim 10, characterized in that the antenna comprises two 
saddle coils driven in parallel to provide a rotating radio frequency field. 



12. The apparatus according to any one of the preceding claims, characterized in that the 
source of radiation produces a component of secondary radiation which is transverse 
to a direction of magnetization of a volume. 

13. The apparatus according to claim 12, characterized in that the secondary radiation 
which is transverse to the direction of the magnetization of the volume comprises 
radio frequency dipole radiation. 



14. The apparatus according to claim 13, characterized in that the radio frequency dipole 
radiation comprises near field radiation. 
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15. The apparatus according to any one of the preceding claims, characterized in that the 
output signal detector comprises a plurality of detectors, 

16. The q>paratus according to claun IS, characterized in &at each of the detectors 
S conq>rises an associated antenna. 

17. The ^aratus according to claim 16, characterized in that flie plurality of detectors 
comprises a fliree dimensional array of antennas. 

0 18. The apparatus according to claim 17, characterized in that each antenna of the three 
dimensional array of antennas responds to a secondary radio frequency field as a 
function of time. 

19. The apparatus according to claim 18, characterized in that the plurality of detectors 
5 are located at spatial locations outside of an object to be imaged. 

20. The apparatus according to any one of the preceding claims, characterized in that the 
output signal detector produces output signals in response to secondary radiation at a 
plurality of spatial locations as a function of time samples at a Nyquist rate in time 

0 over the proton free induction decay, 

21. The apparatus according to any one of the preceding claims, characterized in that the 
first processor further comprises an analog to digital converter to digitize signals of 
secondary radiation. 

5 

22. The apparatus according to any one of the preceding claims, characterized in that the 
first processor comprises a time Fourier transform processor to convert secondary 
signals from the signal output detector over time into Fourier components, each 
having the same frequency, an intensity and a phase angle. 

0 

23. The apparatus according to any one of the preceding claims, characterized in that the 
plurality of Fourier components each having the same frequency, an intensity and a 
phase angle are such that the phase of each component is unique. 
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24. Hie apparatus according to claim 23, characterized in that the unique phase of each 
componeat is provided by a combination of the angle 6 susp^ded between each 
spatial location at which the secondary field is detected and a radial vector, a vector 
from each voxel to each spatial location, and the angle ^ due to a separation distance 

S r between each voxel and each spatial location given by a wavenumber of detected 
field A: times r. 

25. The apparatus according to any one of the preceding claims, characterized in that 
said first processor uses a first component having a phase angle, calculates the phase 

10 angle as a function of spatial position of a first spatial position relative to any other 
spatial location and identifies a component at each spatial location having the 
calculated phase angle. 

26. The apparatus according to claim 25, characterized in that a separate detector 
15 samples each separate spatial locatioa 

27. The apparatus according to claim 26, characterized in that each detector comprises a 
point detector. 

20 28. The apparatus according to any one of the preceding claims, characterized in that the 
source of radiation comprises at least one radio firequency gmerator selected fiom the 
group consisting of klystrons, backward wave oscillators, Gunn diodes, and traveling 
wave tube ampUfiers. 

25 29. The ^paratus according to any one of the preceding claims, characterized in that the 
third processor uses a spatial variation of the secondary radiation over a sample 
volume comprising a plurality of spatial locations at which the radiation is detected 
to detennine a location of the voxels. 



30 30. The apparatus according to claim 29, characterized in that the location of each voxel 
is detennihed independently by the third processor. 



wo 02/16956 PCT/USOl/25954 

116 

31. The qypaiatus according to claim 2, wherdn the image processor and display can 
rotate and display an image from a plurality of desired perspectives. 

32. The apparatus according to claim 2, characterized in ttiat &e image processor and 
5 display can display a reduced dimensional image. 

33. The apparatus according to any one of the preceding claims, characterized in that the 
apparatus provides a two dimensional or three dimensional nuclear magnetic 
resonance image of inanimate or animate objects. 

10 

34. The apparatus according to claim 33, characterized in that the output signal detector 
comprises a nuclear magnetic resonance detector for detecting a time dependent 
•transverse radio frequency field component as a function of time of over a three 
dimensional sample volume outside of an object 

15 

35. The apparatus according to claim 34, characterized in that secondary radiation 
emanating from a magnetized object conq)rises the transverse, time hamionic nuclear 
magnetic resonance magnetic dipole radiation. 

20 36. The ^paratus according to claim 34, characterized in that the time dependent 

transverse radio frequency field component comprises a harmonic oscillation of each 
radio frequency dipole that is equivalent to a dipole rotating in a transverse plane. 

37. The apparatus according to claim 34, characterized in that the apparatus can generate 
25 an image from a three-dimensional map of a transverse resonant radio frequency 
magnetic flux external to an object and said nuclear magnetic resonance detector 
comprises a plurality of detectors at spatial locations, and wherein 

a nuclear magnetic resonance signal of each voxel at a given detector gives rise to 
a Fourier component with a unique phase angle relative to a Fourier 
30 component of other voxels of the object at that detector; 

a set of Fourier components that correspond to the nuclear magnetic 
resonance signal of a given voxel over the detectors is determined; and 

an intensity variation of a transverse radio frequency field over space is 
used to determine a coordinate location of each voxel. 
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38. The apparatus according to claim 37» characterized in that the each detector is 
selectively responsive to a component of the secondary field. 

5 39- The ^paratus according to claim 38, characterized in that a transverse radio 

frequency field is a near field which is a dipole that s^es as a basis element to form 
a unique reconstruction. 

40. The apparatus according to claim 39, characterized in that the eadi detector is 
10 responsive to a component of radio frequency in a transverse plane. 

41. The apparatus according to claim 40, characterized in that each detector is 
responsive to a component of the radio frequency dipole field parallel to the radio 
frequency dipole. 

15 

42. The apparatus according to claim 41, characterized in that each detector is 
responsive to the component of the radio frequency dipole field parallel to the radio 
frequency dipole as a fimction of time and detects a maximum intensity of the signal 
during each time harmonic cycle, 

20 

43. The apparatus according to any one of the preceding claims, characterized in that the 
third processor determines a location of a voxel from signals recorded of a transverse 
dipole radio frequency field. 

25 44. The apparatus according to claim 43, characterized in that a geometric system 
function corresponding to a dipole determines a spatial intensity variations of the 
radio frequency field. 

45. The apparatus according to claim 44, characterized in that the system function is a 
30 band-pass when =k^, the frequency of spatial variation of the recorded secondary 
radiation in the rho direction is equal or comparable to that in the z direction of the 
three dimensional space comprising said plurality of spatial locations. 
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46. The apparatus according to claim 45, characterized in that each voxel is 
reconstnicted indq>6adently in parallel with other voxels using each unique set of 
components associated by phase such that flie scan time is no greater than die nuclear 
fiee induction decay time. 

5 

47. The apparatus according to claim 46, characterized in that all voxels comprising an 
image are reconstructed in parallel. 

48. The apparatus according to claim 46, charactoized in that the apparatus can graerate 
10 high resolution, three-dimensional, real-time anatomical images and ims^es based on 

physiological parameters with Kttle or no deterioration fiom motion artifact. 

49. The apparatus according to claim 41, characterized in that at each point in time, each 
radio frequency dipole rotating in the transverse plane is directed at an angle 6 

1 5 relative to the direction of detection of each detector. 

50. The apparatus according to claim 49, characterized in that the phase angle 0 of the 
radio frequency dipole relative to the direction of detection axis of Qie detector gives 
rise to a phase angle term e"'^ . 

20 

5 1 . The apparatus according to claim 50, characterized in that at each point in time, each 
radio frequency dipole rotating in the transverse plane with a separation distance r 
between each voxel and each spatial location at a detector gives rise to an angle ^ 
given by a wavenumber of the detected secondary radiation * times r directed at the 

25 angle 0 relative to a direction of detection of each detector. 

52. The s^paratus according to claim 51, characterized in that the phase angle ^ of the 
radio frequency dipole relative to a direction of detection axis of a detector gives rise 
to a phase angle term e"^* directed at angle 0 relative to the direction of detection of 

30 each detector and wherein the angle 0 gives rise to a phase angle term e . 

53. The s^aratus according to claim 52, characterized in that a sum of the phase angles, 
kr and 6 , is unique for each voxel at each detector. 
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54. The ^pamtus according to claim 37» characterized in that the unique phase of each 
component is provided by a combination of 

angle 9 suspended between each spatial location at which a secondary nuclear 
5 magnetic resonance field is detected and a radial vector &om each voxel to each 
spatial location, and 

angle ^ due to a separation distance r between each voxel and each spatial 
location given by the waveniunber of detected nuclear magnetic resonance field k 
times r. 

10 

55. The apparatus according to claim 37, further comprising a processor, characterized 
in that a set of Fourier components that correspond to a nuclear magnetic resonance 
signal of a given voxel over detectors is determined by the processor using a position 
of a first detector relative to a different detector to calculated the phase angle of the 

15 different detector relative to the first detector. 

56. The apparatus accoiding to claim 55, characterized in that the processor deteraiines 
the set of Fourier components tiiat correspond to the nuclear magnetic resonance 
signal of a given voxel over the detectors by using the position of each first detector 

20 relative to a different detector to calculated the phase angle of tiie different detector 
relative to the first detectors, and the process is repeated over all of the detectors to 
give the set of intensities of the nuclear magnetic resonance signal over the plurality 
of detectors due to each voxel. 

25 57. The apparatus according to claim 56, characterized in that the set of intensities of the 
nuclear magnetic resonance signal over the plurality of detectors due to each voxel is 
processed to determine a coordinate location of each voxel. 

58. The apparatus according to any one of the preceding claims, fiirther comprises coils 
30 for providing 

, a steady primary field, 
G^, a field gradient in an x-axis direction, 
Gyy a field gradient in a y-axis direction. 
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a radio frequency field, and 

Cr,, a field gradient in a z-axis direction.. 

59. The apparatus according to claim 58« associated amplifiers for driving the coils and 
5 associated control circuits for controlling the coils which provide T, and pulse 

sequences to alter an intensity of the secondary radio fi-equency field as a function of 
time to provide corresponding signals to be processed. 

60. The q)paratus according to claim 9, characterized in that nuclei present in object are 
10 subjected to an additional rotating radio frequency field , which is synchronous 

with precession and their magnetic moments, and magnetization M is caused to 
process about H, and rotate away from the primary field by an angle ^ in a 
coordinate firame i^ch rotates at the Larmor frequency, the precession about 
continues as long as exists, and a final value of 0 which dq)mds on the strength 
15 of H J , which detemiines the precession rate, and the time for which it is turned on is 
such that ^ = 90° that a dqiole is peipendicular to a detector array tiiat is responsive 
to the nuclear magnetic resonance radiation at a plurality of spatial locations extemal 
to the object 

20 61. The apparatus according to any one of the preceding claims, characterized in that the 
magnetizing source is constmcted to provide a magnetic field which permeates the 
object to be imaged and is confined to a volume to be imaged, and the confined 
magnetic field limits the source of signal to the volume of interest, characterized in 
that the volume to be reconstructed is limited to the magnetized volume which sets a 

25 limit to required computations and reduces end efifects of signal originating outside 
of edges of the output signal detector. 

62. The apparatus according to claim 61, characterized in that the apparatus is 
constructed to provide a magnetizing field which limits an imaged volume of an 
30 object by providing a range of Larmor frequencies and data is rejected at a threshold 
Larmor frequency. 
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63. The apparatus according to any one of die preceding claims, characterized in that the 
s^aratus is constructed to provide a magnetizing field which limits an imaged 
volume of an object by providing a range of Lannor fiequencies and data is rejected 
at a threshold Larmor frequency. 

5 

64. The apparatus according to any one of the preceding claims, diaracterized in that the 
output signal detector comprises detector elements disposed on a first axis and 
having corresponding output signals. 

10 65. The apparatus according to claim 64, characterized in that the output signal detector 
includes detector elements disposed on a second axis and having a corresponding 
output signal, and the detectors on the second axis are responsive to the 
corresponding oriented field components of the secondary radiation. 

1 5 66. The apparatus according to any one of the preceding claims, characterized in that 
the third processor comprises at least one reconstruction processor selected from the 
group consisting of a detector processor for processing the detector output signal 
with an inverse matrix, a reiterative processor, a time Fourier transform processor, a 
sine fimction filter processor for detecting the detector output signals, a space and 

20 time Fourier transform processor, and a compensating processor for compensating 
for a distance from the detector along a detector responsive axis. 

67. The apparatus according to any one of the preceding claims, characterized in that the 
signal output detector further comprises a line or plane of detectors and a drive 

25 mechanism that moves the detectors to sample an external secondary radiation over a 
three-dimensional volume. 

68. The apparatus according to claim 67, characterized in that the secondary radiation is 
a radio frequency magnetic field due to nuclear magnetic resonance of nuclei in a 

30 magnetic field produced by the magnetization source and excited by a radiation field 
from the source of radiation: 

69. The apparatus according to claim 68, characterized in that the nuclear magnetic 
resonance active nuclei are protons. 
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70. The apparatus ^xording to claim 19, characterized in that each antenna of the three 
dimensional anay of antennas are impedance matched. 

5 71. The apparatosaccoiding to claim 70, characterized in that each anten^ 

dimensional array of antomas fiu:Qier comprises siq>erconducting reactance elements 
which minimize resistive losses. 

72. The ^paratus according to claim 19» characterized in that cross talk between 

10 antennas is ameliorated or eliminated by time multiplexing the signal detection ov^ 
the array of antennas. 

73. The apparatus according to claim 19, characterized in that the signal output detector 
comprises a synchronous detector wherein the radio frequency field is sampled 

1 5 synchronously so that a phase at any given detector may be related to that at any 
other detector. 

74. The apparatus according to claim 19, characterized in that the signal output detector 
comprises a synchronous detector wherein the radio frequency field is sampled at 

20 known times so that a phase at any given detector may be related to that at any other 
detector. 

75. The apparatus according to claim 19, characterized in that each antenna of the three 
dimensional array of anteimas responds selectively to a component of the secondary 

25 radio frequency field at a spatial location outside of the body as a fimction of time. 

76. The apparatus according to claim 19, characterized in that each antenna of the three 
dimensional array of anteimas responds selectively to a parallel component of the 
secondary radio frequency field at a spatial location outside of an object as a fimction 

30 oftime. 

77. The apparatus according to claim 19, characterized in that each antenna of the three 
dimensional array of anteimas responds selectively to a perpendicular component of 



wo 02/16956 PCT/USOl/25954 

123 

the secondary radio fireqaoicy field at a spatial location outside of an object as a 
function of time. 



78. Apparatus for providing a multidimensional image of inanimate or animate objects 
5 comprising: 

means for producing a magnetic field to magnetize a volume of an object to be 
imaged, the volume comprising a plurality of voxels; 

means for applying a radiation field to the object to be imaged when the object is 
disposed in the magnetic field to produce a secondary radiation field emanating fit>m 
10 the object; 

means for producing ou^ut signals in response to said secondary radiation at a 
plurality of spatial locations outside of the object as a function of time; 

first processor means for determining a pluraUty of Fourier components, each 
having the same iSrequency, an intensity and a phase angle; 
1 5 second pxwessor means for associating the Fourier components due to each voxel 

of the object to be imaged by phase to form a set corresponding to a spatial variation 
of intensity of the secondary radiation due to each voxel at said plurality of spatial 
locations outside of the object; 

third processor means for converting each said set of components into a voxel 
20 location; and 

image processor means for producing an image based on the secondary radiation 
intensity from each voxel and the location of each voxel. 

79. An apparatus for providing a two dimensional or three-dimensional nuclear 
25 magnetic resonance image of an object comprising: 

means for generating an image fi^om a three-dimensional map of a transverse 
resonant radio frequency magnetic flux external to an object having a volume to be 
imaged, the volume comprising a plurality of voxels; and 

processor means for determining the location of each voxel which is a source of 
30 the radio frequency field through spatial variations of an intensity of the nuclear 
magnetic resonance field of a set of components associated by phase. 

80. An ^paratus for providing a two dimensional or three-dimensional nuclear 
magnetic resonance image of an object comprising: 



wo 02/1^956 PCT/US01/2S954 

124 

a magaetizing source for magaetiziiig a volume an obj ect in a first direction, the 
volume comprising a plurality of voxels; 

means for causing nuclear magnetic resonance magnetization to rotate into a 
plane transverse to the first direction; 
5 an output signal detector for detecting die nuclear magnetic resonance field 

rotated into a transverse plane external to the object at a plurality of spatial 
locations; 

a fiirst processor for det^mining a set of components of spatial variation of the 
nuclear magnetic resonance field over a sample volume by association of the 
10 con^nents by phase; 

. a second processor for determining a location of each voxel which is a source of 
the radio fiiequency field through each set of components representative of spatial 
variations of an intensity of the nuclear magnetic resonance field over the sample 
volume, and 

IS an image processor for displaying the voxel locations rq>resentative of 

anatomical images of the object based on selected physiological 



8 1 . An q)paratus for providing a two dim^ional or three-dimensional nuclear 
magnetic resonance image of an object corqprising: 
20 means for magnetizing a volume an object in a first direction, the volume 

' comprising a plurality of voxels; 

means for causing nuclear magnetic resonance magnetization to rotate into a 
plane transverse to the first direction; 

means for detecting the nuclear magnetic resonance field rotated into a transverse 
25 plane external to the object over a sample volume comprising a plurality of spatial 
locations; 

processor means for determining a set of components of spatial variation of the 
nuclear magnetic resonance field over a sample volume by association of the 
components by phase; 

30 processor means for determining a location of each voxel which is a source of the 

radio frequency field through each set of components representative of spatial 
variations of an intensity of the nuclear magnetic resonance field over the sample 
volume, and 
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processor means for displaying die voxel locations representative of anatomical 
images of (he object based on selected physiological parameters. 

82. An qiparatus for providing a two dimensional or three-dimensional nuclear 
S magnetic resonance image of an object comprising: 

means for generating an image from a three-dimensional map of a transverse 
resonant radio fiequency magnetic flux external to an object comprising a volume 
to be imaged, the volume comprising a plurality of voxels; and 
processor means for determining locations of voxels which is a source of a radio 
10 fiequency field througji spatial variations of an intensity of a nuclear magnetic 
resonance field of a set of components associated by phase. 

83. A method for providing a multidimensional image of inanimate and animate objects 
comprising the steps of: 

1 5 placing an object haviag a volume in a magnetic field, the volume comprising a 

plurality of voxels; 

applying a first source of radiation to the object; 

recording a secondary source of radiation emanating from the object at a plurality 
of spatial locations; 

20 forming components of the spatial variation of the secondary source of radiation 

external to the object due to voxels of the object; 
determining the location of voxels ftom the intensity variation of the components; 
generating and displaying the image from the locations and the intraisity of the 
secondary radiation of each voxel. 

25 

84. The method according to claim 83, characterized in that the secondary recorded 
radiation is a three-dimensional transverse resonant radio frequency magnetic flux 
extemal to the object recorded as a ftinction of time. 

30 85. The method according to claim 84, characterized in that the three-dimensional 

transverse resonant radio frequency magnetic flux extemal to the object is nuclear 
magnetic resonance radiation. 
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86. The method accoidmg to any one of claims 84 or 85, further conq>ridng recording 
the secondary source of radiation emanating fiom the object at a plurality of spatial 
locations using a plurality of detectors at each spatial location. 

S 87. The method according to any one of claims 83-86, characterized in that the recorded 
signals fiom each voxel contains components each having a unique phase angle 
relative to other conq>onents fiom difiarrat voxels. 

88. The method according to any one of claims 83-87, characterized in that a time 
10 dependent signal at each detector may be transformed into a series of components 

having intensity and phase data. 

89. The method according to any one of claims 83-88, characterized in that each set of 
components of the nuclear magnetic resonance signal over a sanople volume due to a 

15 given voxel is deteraiined from phase data and detector positions. 

90. The method according to any one of claims 83-89, characterized in that an intensity 
variation of a transverse radio frequency field over the sample volume is used to 
determine the coordinate location of each voxel. 

20 

91. The method according to any one of claims 83-90, characterized in that the radio 
frequency field is a near field which is a dipole that serves as a basis element to form 
a unique reconstruction. 

25 92. The method according to any one of claims 83-91, characterized in that each volume 
element is reconstructed independently in parallel with all other volume elements 
such that a scan time is no greater than a nuclear firee induction decay time. 

93. The method according to any one of claims 83-92, characterized in that said object is 
30 at least one selected &om the group consisting of a human body, an animal, and 

tissue. 

94. A method for providing a multidimensional image of inanimate and animate objects 
comprising the stq)s of: 
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magnefizmg a volume of an object to be imaged, the volume coiiq)rismg a 
plurality of voxels; 

applying a radiation field to the object to be imaged when disposed in the 
magn^c field to produce a secondary radiation field emanating fiom the object; 
S producing output signals fiom a detector in response to said secondary radiation 

at a pliirality of spatial locations outside of ibe object as a Amotion of time; 

determining a plurality of Fourier components each having the same firequency, 
an intensity and a phase angle; 

associating (he Fourier components due to each voxel of the object to be imaged 
10 by phase to form a set corresponding to the spatial variation of intensity of the 

secondary radiation due to each voxel at tiie plurality of spatial locations outside of 
the object; 

converting each set of components into a voxel location; and 
producing an image based on an intensity of the secondary radiation Gx>m each 
15 voxel, the location of each voxel, and plotting the superposition of results for each 
voxel. 

95. A method for providing a two dimensional or three-dimensional nuclear magnetic 
resonance image of an object comprising the steps of: 

20 generating an image firom a three-dimensional map of a transverse resonant radio 

frequency magnetic flux extemal to a volume of an object being scanned, the 
volume comprising a plurality of voxels, characterized in that a nuclear magnetic 
resonance signal of each voxel at any given detector gives rise to a Fourier 
component with a unique phase angle relative to the Fourier component of any other 

25 voxel of the object at a detector; 

a set of Fourier components that correspond to a nuclear magnetic resonance 
signal of a given voxel over detectors is determined; 

an intensity variation of the transverse radio frequency field over space is used to 
determine a coordinate location of each voxel; and 

30 a superposition of results for each voxel is plotted to provide a total image. 

96. A method of providing a multidimensional unage of inanimate and animate objects 
comprising the steps of: 
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aligning magnetic moments of nuclei of a volume of an object to be imaged using 
a primary field, the volume comprising a plurality of voxels; 

further aligning the magnetic moments by a radio fi:equency pulse or series of 
pulses; 

S recording free induction decay signals; 

Fburier transfonning time dependent nuclear magnetic resonance signals to give 
intensity and phase of each Fourier compon^t, characterized in that nuclear 
magnetic resonance signal of each voxel at any given detector gives rise to a Fourier 
component with a unique phase angle relative to the Fourier component of any 
10 other voxel at that detector, 

determining a matrix of Fourier components that correspond to the nuclear 
magnetic resonance signal of a given voxel over the detectors; 
determining matrices for a plurality of voxels; 

det^mining a coordinate location of each voxel from measurements of spatial 
IS variations of a transverse radio frequency field of each corresponding matrix; and 
processing and displaying the position of each voxel with a representation of the 
intensity of the nuclear magnetic resonance signal emitted by each voxel. 

97. A meOiod according to any one of claims 83-96, characterized in that a strength and 
20 duration of rotating Hj radio frequency field that is resonant with the protons of the 

magnetized volume and is oriented perpendicularly to a direction of the magnetizing 
field is applied such that final precession angle of magnetization is 90 °, - 90° , 

such that a radio frequency dipole is transverse to a primary magnetizing field and 
perpendicular to a radio fi:equency magnetic field detector. 

25 

98. The method according to any of claim 83-97, characterized in that nuclear magnetic 
resonance pulse sequences which provide signals for a 7j or 7^ image are applied. 

99. The method according to claim 98, charactmzed in that the nuclear magnetic 

30 resonance pulse sequences which provide the signals for a 2^ or 7^ image comprises 
applying a 90 pulse followed by a series of 180 pulses. 
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100. The melhod according to claim 99, characterized in that the step of flying a 90 
® pulse followed by a series of 180 pulses conq)rises the Carr-Purcell-Meiboom- 
Gill sequence. 

5 101. The method accoiding to any one of claims 96-100, characterized in that tiie step 
of determining the matrix of Fourier components that correspond to the nuclear 
magnetic lesonance signal of a given voxel over the detectors is achieved by using a 
first component having a phase angle and calculating the phase angle as a function of 
spatial position of a first detector relative to any other detector and identifying a 

10 component at each detector having the calculated phase angle. 

102. The method according to any one of clauns 96-100, characterized in that each 
matrix of components associated by phase comprises the intensity variation over the 
sample volume of the radio firequency field of bulk magnetization M of each voxel. 

15 

1 03 . The method according to claim 96, fiirther comprising reconstructing a 
multidimensional image representation according to spatial variations of the detected 
radio fi-equency magnetic field using at least one of a matrix inversion algorithm, a 
reiterative algorithm or Fourier transform algorithm. 

20 

1 04. The method according to claim 1 03, characterized in tiiat a spatial Fourier 
transform algorithm is used. 

105. The method according to claim 104, characterized in tiiat the spatial Fourier 
25 transform algorithm is performed on each set of components over the detector array 

to m^ each bulk magnetization M corresponding to a voxel to a spatial location of 
the imaged object. 

1 06. The method according to claims 1 04 or 1 05, characterized in that the spatial 
30 Fourier transform comprises the steps of : 

recording a nuclear magnetic resonance signal at discrete points in the sampling 
space at position x, y, z as a fimction of time; 

designating each point, x, y, z, as an element in matrix A ; 
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Fourier transfonning time dependent signals to give an intensity and phase of 
each component wherein &e nuclear magnetic resonance signal of each voxel at any 
given detector gives rise to a Fourier conqionent witti a unique phase angle relative 
to the Fouri^ component of any other voxel of the object at that detector 
5 determining a matrix of Fourier components that correspond to the nuclear 

magnetic resonance signal of a givoi voxel ova: the detectors; 

deteimining matrices for all of die voxels; 

detennining a coordinate location of each voxel using the measurements of the 
spatial variations of the transverse radio frequency field of a given matrix; and 
10 plotting a superposition of results for each voxel to provide a total image. 

107. The method according to claim 106, characterized in that the step of determining 
the coordinate location of each voxel using the measm-ements of the spatial variations 
of the transverse radio frequency field of a given matrix comprises the steps of 
15 a.) discrete Fourier transforming each matrix to obtain a matrix B„ ; 

b. ) multiplying each element of each matrix by a corresponding inverse value 
or reciprocal value of the Fourier transfi}rm of a system fimction which gives the 
spatial variation of flie radio frequency field emitted by each voxel as a fimction of 
distance from the voxel to each detector, evaluated at the same frequency as the 

20 element of tiie matrix 4, to form matrix ; 

c. ) generating matrix Z>„ by taldng the discrete inverse Fourier transform of 
matrix C„; 

d. ) multiplying each element of each matrix by the distance squared along the 
z-axis to which the element corresponds to generate the position of the bulk 

25 magnetization M of voxel n which conects the limitation of the sample volume to 

z>0; 

e. ) repeating steps a-d for all matrices A^\ and 

f ) plotting superposition of results for each voxel to provide a total image. 



30 108. The method according to claims 1 06 or 107, characterized in that at each 

sampled external spatial location, x, y, and z in Cartesian coordinates, a magnitude 
B' of the transverse radio frequency magnetic flux emitted by the corresponding z- 
axis aligned magnetic moment of each voxel is given by 
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where the magnetic moment corresponds to die voxel's bulk magnetization M . 



109. The method according to any one of claims 106-108, further comprising the step 
of determining a matrix of Fourier components that correspond to the nuclear 
magnetic resonance signal of a given voxel over the detectors by using a first 
component having a phase angle and calculating the phase angle as a function of 
spatial position of a first detector relative to any other detector and identifying a 
component at each detector having the calculated phase angle. 

1 10. The method according to claim 109, characterized in that each matrix of 
components associated by phase comprises the intoisity variation over the sample 
volume of the radio frequency field of the bulk magnetization M of each voxel. 



10 



IS 111. The method according to claims 106 or 107, characterized in tiiat a system 
function as a function of the detector location relative to each voxel in Cartesian 
coordinates, x, y, z, is by 



20 1 12. The method according to any one of claims 106-1 1 1, further con5)rising the step 
of correcting a point spread of each reconstructed voxel by assigning one voxel 
above a certain threshold with the bulk magnetization M and the other voxels are 
assigned a zero value. 

25 113. The method according to any one of claims 106- 1 1 2, further coii:^)rising the step 
of correcting a finite size of each detector by multiplying each nxmieric value of 
matrix by the corresponding matrix element of the inverse or reciprocal of the 
Fourier transform of a waveform representative of a detector geometry. 



30 



1 14. The method according to claim 106-1 13, characterized in that each matrix 
element of the inverse or reciprocal of the Fourier transform of a waveform 
representative of the detector geometry is given by: 
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1 

smc(^,*r,^) 

* « 

1 IS. The mefhod according to any one of claims 106, 107, or 1 1 1-114, chaiacterized 
in that each volume element is reconstructed indq>endently in parallel with all other 
5 volume elemoits such that a scan time is no greater than a nuclear fiee induction 
decay time. 

1 16. The method according to any one of claims 106, 1 07, or 1 1 1-1 1 5 wherein the 
Fourier transform of the system function is given by: 

117. The method according to any one of claims 106-1 16, characterized in that each 
matrix element for the stq) of multiplying each element of each matrix by the 
distance squared along fhe z-axis to which the element corresponds to generate the 

15 position of the bulk magnetization M of voxel n which corrects the limitation of the 
sample volume to z> 0 is givaaby: 

20 118. The method according to claim 103, characterized in that the step of 

reconstructing the multidimensional image uses a matrix inversion reconstruction 
algorithm. 

119. The method according to claim 118, characterized in that die matrix inversion 
25 reconstruction algorithm comprises the steps of : 

using the geometric system fimction to determine the spatial intensity variation of 
the transverse radio frequency field over a detector array to form a conesponding 
matrix; 

inverting the corresponding matrix; and 
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multiplying the signal over the detector array by the inverted matrix to give voxel 
sources. 



120. The mettiod according to claim 1 19, further conq>rising the steps of 
S a.) determining a coefiSdent for each voxel mathematically or by calibration 

which when multiplied by the bulk magnetization M of each voxel is tiiat voxel's 
contribution to the signal at a given detector with the corresponding unique phase at 
each detector; 

b.) repeating step a for every detectoii 
10 c.) using the coefBci^its to detemdne a matrix which when multiplied by a 

column vector of the bulk magnetization M values of the voxels, gives the signals at 
the detectors; 

d. ) inverting the matrix and storing the result in memory; 

e. ) recording signals as a function of time over the detector array; 

1 5 f.) Fourier transforming the signals as a function of time to give the Fourier 

compon^ts each having an amplitude and a unique phase; and 

g.) multiplying the components by the inverse matrix to generate a bulk 
magnetization M map which is displayed. 

20 121. The method according to claim 120, characterized in that the coefficieat for each 
voxel determined mathematically is given by 

1 V ^[^-^^^f ■^[>^"^2^f -"[^^^^fl 

^iw^Qi; 2. L i — -rn — 7T~r — ttt 

»,.o ii|-4;2* 4/2* p - n,^ J + [y /ijAi J + [z + n,Ai f J 

where K[ji^;^z] is a detector signal at Cartesian coordinate x, y, z, is a constant 
characteristic of each detector, /j, 4, and ^ are the total lengths of the detector array 
25 in the x, y, and z directions. A:,, Aj, and are spacing of elements of the detector 
array in the x, y, and z directions, and n^ y /i^, and n^ are integer indices of the 
elements of the detector array in the x, y, and z directions. 



30 122. The method according to claim 118, further comprising the stqps of 
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det^mining a matrix of Fourior components that cotrespond to the nuclear 
magnetic resonance signal of a given voxel over the detectors; 

detennining conesponding matrices for a plurality of voxels; 

detemiining a coordinate location of each voxel from measurements of spatial 
5 variations of the transverse radio frequency field of a givai matrix viierein each 

matrix of components associated by phase comprises the intensity variation over the 
sample volume of the radio frequency field of the bulk magnetization M of each 
voxel; and 

multiplying each matrix of components by the inverse matrix of system fimction 
10 to generate the bulk magnetization M msp. 

123. The method according to claim 122, further comprising the steps of 

a.) correcting a point spread of a reconstracted voxel by assigning one voxel 
above a certain threshold with the bulk magnetization M ; 
1 5 b.) assigning other voxels with a zero value; 

c. ) repeating steps a-b for other voxels; and 

d. ) plotting and displaying the superposition of the separate maps for each voxel 
which comprises the total bulk magnetization map as a nuclear magnetic resonance 
image. 

20 

124. The method according to claim 103, characterized in that said multidimensional 
image is reconstructed using a reiterative reconstruction algorithm. 

125. The method according to claim 124, characterized in that the reiterative 
25 reconstruction algorithm comprises the steps of 

a. ) using the geometric system function to determine a system of linear equations 
which gives the intensity, spatial variation, and phase of radio frequency field over a 
sample volume wherein the system of linear equations gives the signal and phase 
fix)m each voxel at each detector based on the bulk magnetization M value of each 

30 voxel and the position of the voxel relative to the detector; 

b. ) Fourier transforming the signal as a function of time to give Fourier 
components each having an amplitude and a unique phase; 

c. ) determining weighting coefiBcients based on the linear equations; 

d. ) det^mining coefficients for each voxel mathematically or by caUbration; 
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e. ) estimating bulk magnetization M for each voxel; 

f. ) calculating signals at each phase at each detecton 

g. ) comparing a result of multiplying the bulk magnetization M value of each 



voxel times its weighting coefficient and its calculated phase at a given detector to a 



5 



measured signal and phase; 

h.) making a correction to M of each voxel to give a second, or recomputed. 



estimate for M of each voxel; 

i.) computing the signal value from flus second estimate and making corrections 



10 



according to steps f-^; 

j.) repeating steps f-i until the correction for each reiteration approaches a 



predefined limit which serves to indicate that the reconstruction is within reasonable 
limits of error, 

Ic) repeating steps e-j of the procedure for other detectors; 
and 



20 where ^[ji^j^yz] is the detector signal at Cartesian coordinate x, y, z, Q is a constant 
characteristic of each detector, /,, Z^, and are flie total lengths of the detector array 
in the x, y, and z directions, kjy and are the spacing of the elements of the 
detector array m the x, y, and z directions, and n,, n^^ and n^ are the integer indices of 
the elements of the detector array in the x, y, and z directions. 

127. A method of providing a multidimensional image of inanimate and animate 
objects comprising tiie steps of : 

a.) recording secondary radiation signals corresponding to a transverse radio 



15 



1.) plotting and displaying a bulk magnetization map. 



126. The method according to claim 125, wherein the coefficient for each voxel 
determined mathematically is given by 




frequency magnetic field in response to a first radiation q)plied to a magnetized 
object over a sample volume comprising a pluraUty of voxels; 

b.) forming a matrix of signals at a plurality of spatial positions each 

corresponding to a detector element of a detector array; 
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c.) transfoiming the matrix nummc values accoiding to a Fourier transform to 
give an intensity and phase of eadi component wheiem a nuclear magnetic 
resonance signal of each voxel at any given detector gives rise to a Fourier 
component with a unique phase angle relative to the Fourier component of any other 



d) determining a matrix of Fourier conq)oneiits that correspond to the nuclear 
magnetic resonance signal of a given voxel over the detectors; 

e.) multiplying each numeric value of the transform array by a value which is an 
inverse or reciprocal of the Fourier transfonn of a waveform corresponding to an 
operation performed in step b; 

f ) multiplying each numeric value of the discrete spectrum matrix by a value 
cotresponding to an inverse Fourier transfonn of a function representative of the 
steps of exciting and detecting and providing a system corrected matrix; 

g. ) inverse-transfoiming the system corrected matrix according to a 
multidimensional discrete inverse Fourier transform providing a voxel location; 

h. ) correcting each element in the transformed matrix for any asymmetries in the 
dimensions of the sample volume; 

i. ) repeating steps c-h for other matrices, each corresponding to a distinct voxel; 
and 

j.) superimposing the results Scorn steps a-h, plotting and displaying the image. 

128. The method according to claim 127, characterized in that at each sampled 

external spatial location, x, y, and z in Cartesian coordinates, the magnitude B of the 
transverse radio frequency magnetic flux emitted by the corresponding z-axis aligned 
magnetic moment of each voxel is given by 



where the magnetic moment corresponds to the voxel's bulk magnetization M . 

129. The method according to claims 127 or 128, characterized in that the step of 
determining the matrix of Fourier components that correspond to the nuclear 
magnetic resonance signal of a given voxel over the detectors is achieved by using a 
first component having a phase angle and calculating the phase angle as a function of 



5 



voxel of the object at that detector; 
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spatial position of a fiist detector relative to any other detector and identifying die 
conqKment at each detector having the calculated phase angle. 

130. The method according to claim 129, characterized in that each nmtrix of 

S components associated by phase con^rises an intensity variation over the sample . 
volume of the radio frequency field of the bulk magnetization M of each voxel 

131. The method according to any one of claims 127-130, characterized in tbat a 
system function as a function of the detector location relative to each voxel in 

10 Cartesian coordinates, x, y, z, is given by 
^. 2z^j-x^-/ 

132. The method according to any one of claims 127-131, further comprising the step 
of correcting a point spread of each reconstructed voxel by assigning a voxel above a 

1 S certain threshold with the bulk magnetization M and other voxels are assigned a zero 
value. 

133. The method according to any one of claims 127-132, further comprising step e of 
correcting a finite size of each detector by multiplying each numeric value of the 

20 matrix fonned in step d by the corresponding matrix element of the inverse of the 
Fourier transform of a waveform representative of a detector geometry. 

134. The method according to claim 133, characterized in that each matrix element of 
an inverse of the Fourier transform of a waveform representative of the detector 

25 geometry is given by formula 

1 

sinc(^,*j„*,) 

135. The method according to any one of claims 127-135, characterized in that each 
volume element is reconstructed independently in parallel with all other volume 

30 elements such that a scan time is no greater than a nuclear free induction decay time. 
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136. The method according to any one of claims 127-13S, characterized in that the 
step f of multiplying each numeric value of the discrete spectrum matrix by a value 
corresponding to the reciprocal of the Fourier transform of a fimction rq)resentative 
of the steps of exciting and detecting and providing a system corrected matrix 
5 wherein the discrete values of the Fourier transform of the system function matrix are 
given by: 



137. The method according to any one of claims 127-136, characterized in that the 
10 step h of correcting each element in the transformed matrix for any asymmetries in 
dimensions of the sample volume wherein the discrete values of are given by: 




138. The method according to claim 103, further comprising the steps of using the 
1 5 matrix inv^ion algorithm or the reiterative algorithm in combination with the 

Fourier Transform Algorithm. 

139. The method according to claim 103, further comprising the steps of using the 
matrix inversion algorithm and the reiterative algoritiun in combination with the 

20 Fourier Transform Algorithm. 

140. The method according to claim 103, further comprising using the Fourier 
Transform Algorithm as input for the first reiteration of a reiterative algoritbixL 

25 141. The method according to claim 140, further comprising the steps of: 

using the Fourier Transform Algorithm to calculate a series of first bulk 
magnetization images for each matrix of components associated by phase; 

detecting at least one peaks or edge of each image by using a cut-off level for the 
bulk magnetization assigned to a voxel or the value of the derivative of the bulk 
30 magnetization as a function of voxel position, respectively^ 
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setting the bulk magnetization of the voxel or voxels corresponding to the peak or 
edge equal to the peak bulk magnetization and assigning the remaining voxels a zero 
value; 

calculating a signal at each detector using a system function; 
5 comparing each calculated signal to the actual signal and deteraiining a 

difference or enor; 

projecting the error over the voxels using the Fbuner Transform Algorithm, 
Matrix Inversion Algorithm, or Reiterative Algorithm; 

adding back projected errors to a preceding conected senes of bulk 
0 magnetization mq>s; and 

performing the steps reiteratively until it stabilizes. 

142. The method according to claim 96, characterized in that a J, image is produced 
by a applying at least one pulse sequence that inverts the magnetization and records 

5 the relaxation. 

143. The method according to claim 142, further comprising producing a Jj image by 
the steps of: 

switching on radio frequency receivers to follow a decay following a nuclear 
5 excitation comprising a H, pulse; 

examining a time dependait bulk magnetization in the direction of a 

primary field H^, at a time after an inverting pulse by applying another pulse 
equivalent to a rotation by 90° after waiting the time following the initial 
inversion, 

5 characterized in that the 90° pulse puts the z magnetization (f^) into the 
transverse plane for observation, and 

changing the waitmg time to allow for observation of M^(f) at different times 
during relaxation. 

) 144. The method according to claim 96, characterized in that a image is produced 
by a applying at least one pulse sequence that flips a magnetization vector into a 
transverse plane and records a transverse relaxation by producing at least one spin- 
echo. 
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145. The meOiod according to claim 144, fiirfha: comprising producing a image by 
the steps of 

flipping the magnetization initially aligned along the -axis into the transverse 
S plane by q)plying a H, pulse; 

recovering the relaxation from a FID by applying a 180° H, pulse at time 

along the yn-axis where the relaxation occurs in time 

recording flie relaxation after an additional time t^^ characterized in that the 

total time ek^sed after the 90° pulse v&2t^-t^ which is the spin-echo time and 
0 the dq)hasing due to the static inhomogeneity of Hq is canceled out, whereby the 

relaxation is due to processes that create relaxation; 

applying a plurality of 180° pulses to form a pluraUty of spin echoes; and 
determining the relaxation from the envelop of the maximum amplitude of the 

spin echoes. 

5 

146. The method according to claim 145, further comprising a pulse sequence to give 
the data for a image known as the Carr-Purcell-Meiboom-Gill sequence 
comprises applying a 90° pulse along the x^-axis followed by a series of 180° pulses 
along the y^-axis at times + Int^ where n is an integer including zero. 

Q 

147. A method of providing a multidimensional nuclear magnetic resonance image of 
an animate or inanimate object comprising the steps of: 

placing a voliune of an object having nuclei in a static magnetic field H,,, the 
volume comprising a plurality of voxels, and then subjected them to an additional 
5 rotating radio frequency field E[| , which is synchronous with their precession, such 
that the bulk magnetization M of each voxel of the object processes about Hj and 
rotates away fixmi the primary field by an angle ; 

measuring of the intensity of die radio fi^uency signal over time and space 
following r, and/or encoding pulses wherein the magnitude of M is a maximum 
) initially and decays with time by emission of the same raultipolarity radiation that it 
absorbed and by transfer of energy to the surrounding lattice and the intensity of the 
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radiation is a fimction of M and the coordinate position relative to the radio 
frequency emitting voxel; 

Fouii^ transforming the signal as a function of time at a given detector position 
to give the cony)onents each having an amplitude and a unique phase; 
S determining each set of components of the nuclear magnetic resonance signal 

over the sample volimie due to a given voxel from the phase data and the detector 
positions; 

detamining the location of each voxel through the spatial variations of the 
intensity of the transverse nuclear magnetic resonance field of the set of components 
0 associated by phase; and 

superimposing the results of the determination of the location of M of each voxel 
and the plotting and displaying the resulting image. 

148. The method according to claim 96, characterized in that the plurality of Fourier 
5 components each having the same firequency, an intensity and a phase angle each 

have a unique phase. 

149. The method according to claim 148, characterized in that the unique phase of 
each component is provided by a combination of 

0 the angle 6 suspended between each spatial location at which the secondary field 

is detected and a radial vector firom each voxel to each spatial location; and 

the angle ^ due to a separation distance r between each voxel and each spatial 
location given by a wavenumber of a detected field k times r . 

5 1 50. The method according to claim 87, characterized in that the radio firequency field 
is sampled synchronously so that the phase at any given detector may be related to 
that at any other detector. 



151. The method according to claim 87, characterized in that the radio fi-equency field 
is sampled at known times so that the phase at any given detector may be related to 
that at any other detector. 
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1 52. The method according to claim 96, characterized in tfiat the processor which 
determines a plurality of Fourier components each having the same frequency, an 
intensity and a phase angle uses a first component having a phase angle and 
calculates die phase angle as a function of spatial position of the first spatial position 

5 relative to any other spatial location and identifies the component at each spatial 
location having the calculated phase angle. 

153. The method according to claim 1 52, characterized in that the voxels that are on 
the same axis with a given detector or a plurality of detectors that align vertically 

10 with one or more voxels is determined 

1 54. The method according to claim 153, characterized in that the phase angle is a 
function of the distance r of the detector from the voxel. 

15 155. The method according to claim 154, characterized in that the components may 
be identified by the presence of at least one component of a phase given by ^ + kJSr 
where is the phase angle of the first component, Ar is the distance of the second 
detector relative to the first of the detector, and k is the wavenumber of the nuclear 
magnetic resonance radiation emanating firom the object wherein this detector(s) is 

20 defined as having a zero angle 6 . 

156. The method according to claim 155, characterized in that a phase angle ^ of a 
second component at a second detector aligned on the detection axis corresponding 
to a first component at a first detector having a phase angle is given by 
25 ^=^, +AAr 

where Ar is the distance of the second detector relative to the first of the detector 
and ^ is the wavenumber of the nuclear magaetic resonance radiation emanating 
fipom the object. 

30 157. The method according to claim 156, characterized in that, &om the components 
identified as coaxial with a detector, the phases of the components from the voxel are 
determined at the other detectors Vising the relative position of other detectors of the 
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sample volume relative to the first detector and ^, of the component at the first 
detector. 

158. The method according to claim 157, characterized in that the phase 0 due to a 
sq)aration distance between a voxel and a detector of r is given by the wavenumber 
of the radio firequency field k times r , 

where the waven\raiber k is given by 
the phase angel ^ is given by 

where Ax, Ay, and Az is the separation distance between a voxel and a detector in 
the X, and z directions, respectively, and 

in the case that a two dimensional slice is imaged in the transverse plane, die 
phase angel ^ is given by 

159. The method according to claim 158, characterized in that all voxels which have 
the same radial distance firom a given detector will have the same phase 0 . 

20 

160. The method according to claim 159, characterized in that at each case that two 
detectors have the same phase ^, since the radial distances r of two voxels to a 
given detector axe equivalent, a unique angle 0 suspended between the direction of 
the detector and the vector along r exists. 

25 

161. The method according to claim 160, characterized in that the radio frequency 
dipoles of all voxels are time synchronous. 

162. The method according to claim 96, characterized in that the nuclear magnetic 
30 resonance pulsing and detection may be synchronized. 
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163. The method accoiduig to claim 96, characterized in that at least one standard 
dipole may be used as reference phase to set the phases of the radio frequency 
dipoles of the voxels. 

5 164. The method according to claim 163, charactrnzed in that all voxels that are not 
collinear with the axis of a detector posses an angle 0 suspended between the 
direction of the detector and the radial vector from die dipole to the detector. 

1 65. The method according to claim 1 64, characterized in that since the radio 

10 frequency dipoles are time synchronous, at each point in time, the angle coiresponds 
to a phase angle 9 of the radio frequency signal from each voxel at each detector. 

166. The method according to claim 165, characterized in that in the case that a two 
dimensional slice is imaged in the transverse plane, the phase angle 0 is given by 



167. The method according to claim 166, characterized in that the total unique phase 
angle ^j. of the radio frequency signal from each voxel at each detector is given as 
the sum of ^ and d given by 



15 
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